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ABSTRACT 


Genetic algorithm (GA) is one of the widely used optimization methods. The efficiency of 
GA and its effectiveness depend to a great extent on the selection of its control parameters. 
One of the widely studied performance criteria for GAs is the on-line performance measure 
which largely depends on parameters such as population size, crossover probability, mutation 
probability, selection operator and crossover operator. 

In the present study, variation reduction of the on-line performance measure of a sim- 
ple genetic algorithm is done by the response surface methodology. The results on four 
different test functions suggest that proportionate selection with uniform crossover always 
gives consistent on-line performance of GA. For every function, the population size, crossover 
probability and mutation probability are found based on minimizing the variation in on-line 
performance measure. Since no unique values of these parameters are found to be optimal 
in all functions, it is not possible to conclude that a particular value of population size, 
crossover probability and mutation probability will always produce minimum variation in 
the on-line performance. However in each case the relative importance of these parameters 
are also found. It is observed that for a unimodal function, mutation is an important op- 
erator if used with proportionate selection and uniform crossover. For a function having 
a flat optimal basin, a smaller population with proportionate selection is the best. For a 
multimodal function, a large population with a large crossover probability is needed. This 
study suggests that before trying to solve a problem using GA, the procedure presented in 
this thesis may be applied to got consistent performance of GAs, by making the obtained 
solutions independent of initial population. This thesis has also outlined a methodology to 
obtain GA parameters which will produce consistent performance. 
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Chapter 1 
Introduction 


The performance of a genetic algorithm (GA) depends on its parameters such as population 
size, genetic operator, operator i)robability and others. Often the performance depends on 
the chosen initial population. Although most theoretical studies use Markov chain modeling 
of the GA process, they provide little guidance for the proper selection of above parameters. 
Given this, it is surprising that only few studies have been done in this area. Realizing that 
there is a need for a more thorough investigation, in the present study optimization and 
variation reduction of on-line performance of the GA is done based on methods suggested by 
Taylor[3]. But before we discuss the present work, let us review the past efforts in designing 
GAs. 

De .long K('t out to comiiare GAs with the best gradient search techni(iues. In that study he 
investigated the effect of ]>opulation size, crossover probability and mutation probability on 
on-line performance and off-line performance of GA on a number of test functions (De Jong, 
1975). The following wa.s his conclusion: 

Increasing the population size was shown to reduce the stochastic 
effects [of random sampling on a finite population] and improve 
long term performance at the expense of slower initial response. 

Increasing mutation rate was seen to improve off-line performance 
at the expense of on-line performance. Reducing the crossover rate 
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resulted in an overall improvement in performance, suggesting that 
producing a generation of completely new individual was too high 
a sampling rate [5]. 

Based on his limited study he suggested the following ranges of GA parameters for a good 
on-line performance of GA: 

Population size 50-100 

Crossover rate 0.6 

Mutation rate 0.001 

De Jong also investigated the effect of multi-point crossover on his test suite. He 
observed that a GA with multiple cross sites exhibits different schema disruption behavior, 
but no significant empirical difference was observed due to different crossover sites. 

A more robust approach was suggested by Grefenstette (1986) when he used a GA to address 
the control parameter problem by treating it as a meta-problem. That is, a meta-GA was 
used to locate parameter sets which themselves were used for GA searches on a set of tasks 
( the De Jong suite). His recommended values were [5]: 

Population size 30 

Crossover rate 0.95 

Mutation rate 0.01 

It is interesting to note that he recommended a smaller population size and much higher 
crossover and mutation rates than did De Jong. 

Goldberg (1985) did theoretical investigations to find an optimal population size for a par- 
ticular problem. Based on his simple probabilistic model of schema processing he proposed 

pop = 1.65 X 


( 1 . 1 ) 




3 


that suggests population sizes of 130, 557, 2389 and 10244 for strings of length 30, 40, 50 
and 60 respectivel}^ It is clear that this method suggests very large population for problems 
with moderately long chromosomes. It is also surprising that population size always depends 
on the string length, and not on the complexity of the problem. Realizing this, Goldberg, 
Deb and Clark(1992) devised a different population sizing equation based on signal to noise 
ratio of schema evaluation. 

Robertson’s (1988) investigation of population size on parallel machines found that GA’s 
ptu forinaiu'e nionotonically incn'ased with i)opulation size. But in serial macliities or otlier 
settings wdiere there is a fixed cost increment for each population member; we do not expect 
the ijcrformaiuai to continuously increase with the population size [5]. 

Shaffer, Caruana, Eshclman and Das [5] did a complete factorial design experiments and 
analysis of variance (ANOVA) was used to identify and quantify the influence of these control 
parameters, lie concluded the following ranges of parameters for good on-line performance: 

Population size 20-30 

Crossover rate 0.75-0.95 

Mutation rate 0.005-0.01 

He suggested two point crossovers no worse and sometimes better than one point crossover. 

Ill the all the above study the aim was to find optimum on-line performance. In none 
of the approach variation reduction of on-line performance was done due to change in the 
initial population and also no one took into account the effect of interaction between these 
parameters. 

Bagchi and Deb [8] used design of experiment approach to calibrate the GA parameters. One 
of their important conclusions was that the non-linear interaction among GA parameters 
exist in a GA simulation. Hence the parameters cannot be optimized by studying them 
one at a time. For their suite of problems, they suggested to use any one of the following 
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combination of GA parameters : 

Population size = 0.7 xsirmgleTigth, crossover probability = 0.75, mutation probability 
= 3/string length. 

Population size =1.3 x stringlength, crossover probability = 0.95, mutation probability 
= 3/string length. 

Population size =2.7 x stringlength, crossover probability = 0.5, mutation probability 
= 1/string length. 

Further for problems in which equality constraints are handled by using penalty func- 
tion, the following combination seemed to work well: 

Population size =1.3 x stringlength, crossover probability = 0.6-0. 8, mutation proba- 
bility < l/string length. 

The difficulty of the design of experiment approach is that it will give the optimum perfor- 
mance at one of the experimental settings. So one can not say that suggested combination 
of GA parameters will give global optimum performance. 

In this study, we not only investigate the effect of GA parameters on the on-line performance, 
but also find optimal settings of parameters, so that the performance of GA is less sensitive 
to the initial iiopnlation. l'\)r tlw GA to be a successful tool in science and engineering 
ap])lication, a GA must provide similar yet optimal solution with different initial poiiulations. 
We use the concept of variation reduction method of robust design (Taylor, 1991), to find 
optimal GA parameters and show the robustness of the GA by showing results on three 
functions and one engineering design problem. 

Chapter 2 describes the variation reduction technique. Chapter 3 outlines GA proce- 
dure in brief, and how the variation reduction technique is applied for designing GA. Chapter 
4 shows the results on three functions and one engineering design problem. Chapter 5 gives 
a number of conclusions of the study. 




Chapter 2 

Variation Reduction Technique 


2.1 Cause of Variation 


A system input/output model can be described by Figure 2.1. 


NOISE 



Figure 2.1; System input/output model 

The system produces certain output, be they product, invoices, sales and so forth. In 
producing these output the system requires certain inputs like materials and workforce in 
manufacturing. The system may also be affected by other inputs such as environment which 
is called noise. The system acts on the inputs to try to achieve the desired output. 

The system input/output model implies that there is some sort of relationship between 
the inputs and outputs. Given this relationship, often in the form of an equation or plot, the 
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variation transmitted to an output variable by an input variable can be determined. Figure 
2.2 shows how this might be done. 



Figure 2.2; Transmission of variation 

There are many inputs that might affect the system. Those inputs are called candidate 
input variables. Out of these candidate input variables, not every variable will actually have 
an affect on the product or process. So those inputs that can affect the output, either by 
affecting the average or by contributing to the variation in the output are called key input 
variable. 

An important part of reducing variation is the identification of those key input variables 
that make significant contribution to the variation. The variation in the output is resulted 
by two types;- 

1. Targeting the key input variable may produce variation in the output. These variables 

are called variation adjustment parameters (VAP). 

2. Variation given to the key input variables may also result in variation of output. These 

key input variables are called variation inducing parameter (VIP) . 

Optimization and variation reduction requires to find key input variables, VIPs and VAPs. 
Studying the input’s effect is one of the approaches for identifying key input variables, VIPs 
and VAPs. 
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2.2 Studying the Input’s Effect 

When studying the input’s effects, one wants to identify the key input variables, understand 
their effects on the average and use this understanding to optimize the average. However, 
one also wants to study the input’s effect in order to reduce variation. There are two different 
approaches for accomplishing this task: 

1. Transmission of variation : The variation transmitted by the inputs to the output can 

be calculated based on the relationship between the inputs and output. This relation 
of variation is used to find VAPs and VIPs. Section 2.2.2 deals with this method. 

2. Direct approach : The effect that changes of the inputs have on the variation of the 

output is observed directly. The equation of variation is found out by fitting the best 
surface through the log of standard deviation. This equation of variation is used to 
find VAPs. Section 2.2.3 deals with this method. 

The detection of key input variables is done by the screening experiment. Taylor 
(1995) has described in detail about the procedure for designing and analyzing screening 
experiments. The limitation of screening experiment is that it cannot calculate the non-linear 
effect of inputs, because it uses only two levels of input. Although screening experiment with 
some more trials including center points of the inputs can be used for calculating non-linear 
effects. The advantage of screening experiment is that large number ( upto 16 ) of inputs 

i 

can be analyzed with fewer number of trials. 

If the number of trials is less then full factorial design of experiment can be done. Here 
for 3 inputs, 27 trials experiment is proposed. This 27-triai experiment is analyzed in the 
same way as Taylor (1995) proposed, and will be used for designing GA. 

Table 2.1 provides the information required to determine where to set each input on 
each of the 27 trials. For each trial 5 to 10 runs are taken. 


2.2 Studying the Input’s EfFcct 
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A 

B 

c 

A-B 

A*C 

B*C 


Input variables corresponding 

to column 


Table 2.1: 27-trial experiment 

First three columns are used for each of the three inputs. The columns corresponding 
to the inputs contain 1, 2 and 3. 2 means to set the input to its high value, 1 means to set 
the input to its low value and 3 means to set the input to its center value. Fourth, fifth and 
sixth column shows the interaction A*B, A*C and B*C respectively. 0, 1 and 2 in these 
columns are levels of interaction. The way to calculate level of column corresponding to 
(A*B) is : ' 

Put 1 if level of parameter A and parameter B is not same i.e. A is at level 1 and B is 
at level 2 or vice versa. Put 2 if level of parameter A and parameter B is same i.e. A is at 
level 1 and B is at level 1 or A is at level 2 and B is also at level 2. Put 0 if center level 3 is 
pn'S('nt in any of paraiiu'ter A or B. 

For each trial average and standard deviation of the output is calculated. 

Analyzing the experiment 

Analyzing the experiment consists of calculating the effect of each parameter, non-liner effect 
of each parameter and the group of interactions, sorting the effects in order of importance. 
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and deciding which effects are significant. 

Effect on average of individual parameter = average of trials at level 2 — average of trials at 
level 1 

Effect on average of interaction = average of trials at level 2 — average of trials at level 1 

Non-linear Effect on average = average of trials at level 1 and 2 — average of trials at level 
3. 

Thus by analyzing the experiment, we have the effect of each parameter, interaction and 
non-linear effect of each parameter on the system’s performance. Out of these, the parameter 
wit.h small ('Ifecd.s an' neglect, ('d. The rcunaining parametius are the key input variables. 

Similar to the effect on average, effect on standard deviation of each parameter, their 
interactions and their non-linear effect arc also found. The parameter having smaller ef- 
fects are neglected. The remaining parameter which affects standard deviation are variation 
adjustment parameters. 

2.2.1 Response Surface Studies 

With all the data collected above, and with the terms having larger effect, bestfit surface is 
drawn through the average. 

The eciuatiou for the bestfit surface will be of the form 

Y = Oq -|- Oi -I- 02-^2 “b U3-X^3 T -f- 05^2 + 06-^3 T (^ 7 X 1 X 2 + U 8 W 2 W 3 

TagWiA'a (2.1) 

Where Y is the output, Xi represents the input, and Uj is a constant to be determined based 
on the data. 

In the above polynomial equation only 2^^ order terms are used. The terms with 
higher order can be included but then complexity increases. Hence for simplicity 2"*^ order 
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polynomial is used. 

Bi'stfit surface' l.lirough log of standard deviation is also found out, with the parameters 
affecting the standard deviation. 

The analysis described above shows how to find key input variable, what will be the 
terms in bestfit surface through average and what is the bestfit equation through the average. 

The analysis also shows how to find variation adjustment parameters , what will be the 
terms in bestfit surface through log of standard deviation and what will be equation through 
log of standard deviation. 

2.2.2 Variation Transniission Analysis 

The bestfit surface for the average will be of the form 
Y = g{Xr,X2,X,,....,Xn) 

Let us use the following notations: 

gi{Xi,X2, Xa, ....,Xn) = first partial with respect to Xi 

.f/2(Ai, A'2, A'a, ...., A'„) = first partial with respect to X2 


f;„(A'], A^2! A';}, ...., An) = first partial with respect to AT. 
^gi(AT, A2, AT, A„) = Second partial with respect to Xi 
52,2(AT, A2, A3, A„) = Second partial with respect to A2 


9 n,n{Xi, X2, Xa, ■■■■, Xn) = Sccoud partial with respect to An 
^i, 2(AT, AT, AT, AT) = Second partial with respect to Ai and A2 
5 i, 3(AT, A2, AT, AT) = Second partial with respect to Xi and A3 
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gn-i^ri{Xi,X 2 ,Xz., Xn) = Second partial with respect to X„_i and 
Then the squared standard deviation is approximately equal to (Taylor, 1991) 

o] - + .... + gnikjh, --Xnf + 

0.5 X gi,l{ti,t2, ••••, tn)^crf + 0.5 X 52,2(^11 ^ 2 ) -•••) in)^cr 2 T 
.... + 0.5 X 5 n,n(^l) ^2) ••••) + 5 l, 2 (^l) h, •■■■■, 'tnY <^l(^2 + 

<7l,3(ilj ^2) . — + .. + gn-l,nitl, h, • — (2-2) 

Where t\,t 2 ..., are the value of key input parameter and ai,a 2 ..., are the standard deviations 
of the key input parameter. In the above equation ai,cr2.. are known because we know 
the variation of the key input variable. The above equation is minimized for reducing the 
variation. The optimum value of ti,t 2 .. will be required at which we have to target that 
parameter (This is called the parameter design) for reducing the variation. These parameters 
are variation adjustment parameters (VAPs). 

Tolerance Design 

The purpose of tolerance design is to select tolerance for each of the inputs around the target 
selected during parameter design. As a part of the parameter design, an equation for the 
standard deviation of the output is derived. This equation expresses the standard deviation 
of the output as a function of the input’s target t(s and the input’s standard deviation a[s. 

t 

Oy = 02, ..,On) ( 2 - 3 ) 

Putting the value of t[s into this equation simplifies the equation so that Oy is simply function 
of the input’s standard deviation. 

Oy = f {0-1,02, .., On) ( 2 - 4 ) 

The transmitted contribution of an individual input an is 
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O-tl = \//(0'l,0-2,..,Cr„)2 - /(0,(72, 


crt2 = \//(o'i,o-2,-,o’n)^ -/(o-i,0,..,crn)2 (2.5) 

0-tn = A//(<7l,Cr2,..,Cr„)2 - /((Ji, Cr2, 0)^ 

These standard deviation of the inputs transmitted variation can be ranked. Those inputs 
with largest cr^s are the variation inducing parameters (VIPs). 

The tolerance on these VIPs are tightened to reduce the variation. Those key input variables 
which are not VIPs are found out from the bestfit surface for average, such that average is 
minimized. 

The procedure of the variation transmission analysis can be summarized as: 


1. 27-trial experiment is done. The analysis of experiment is done as described in previous 

section. Then response surface study is done. 

2. The bestfit equation for average is obtained from response surface studies, equation 2.1. 

3. The equation for square of standard deviation can be obtained from equation 2.2. 

4. By putting the standard deviation of each input in equation 2.2 and minimizing it we 

will get the value of parameter . This U is the value at which the corresponding 
parameter should be targeted. These parameters are called (VAPs). 

5. These t[s are put into equation 2.2, which results in equation 2.4. Equation 2.5 is used to 

calculate transmitted variation of each input. The parameter having large transmitted 
variation are VIPs. , 





2.2 Studying the Input’s Effect 


13 


2.2.3 Direct Approach 

The analysis involves analyzing the standard deviation just as one would analyze the average. 
Procedure for the direct approach is as follows: 

1. 27-trial experiment is done. The analysis of experiment is done as described in previous 
section. I’hen r('sj)onse surface study is done. 

2 The bestfit equation for the standard deviation is obtained from response surface studies. 

3 The e<ination for bestfit surface through standard deviation is minimized, and the value 

of the parameters for minimum standard deviation are found out. These parameters 
arc VAPs. The parameters has to be varied around this value of VAPs. 

4 The value of VAPs is put into equation for average. The resulting equation is minimized. 

The value of the parameter obtained after minimization is to where, we have to set it 
for optimization of average. 

The following is the limitation of the direct approach 

1. Since only fc'w samples are taken for each e.xperiment, standard deviation calculated may 

not r('pr('sent acl.ual st.aiulard deviat.ion. In that case some of the VAPs may be missed. 

2. The direct approach ident.ifios VAPs and helps to target them to minimize the variation. 

The direct approach can’t isolate the contribution of individual inputs required for 
tolerance design. Thus it can’t find VIPs. 

2.2.4 Comparison between Direct Approach and Variation Trans- 
mission Analysis 

Despite limitation in direct approach method, using this method is certainly worthwhile. As 
part of any analysis the standard deviation should be analyzed as well as the average. The 
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added information is essentially free, as no additional data are required. Any VAPs that are 
identified can result in valuable improvements. 

However one should not rely solely on direct observation. It is still worthwhile to 
perform a variation transmission analysis. Frequently, VAPs are identified that are missed 
by direct observation. Further, the variation transmission analysis can be used to perform 
tolerance design including identifying VIPs. 

Let’s turn to example to see the difference exactly. Suppose we assume input output relation 
as: 

r = 8. 125 - O.'i.V -I- 0.185A'2 - O.OOGA'^ (2.6) 

A plot of the effect of X on Y is shown in Figure 2.3. 

The variation of Y is approximately 

ay = ^(-0.4 + 0.37t^ - 0.018t2)V2 + 0.5(0.37 - (2.7) 

Further assume that a^ = 0.2. A plot of the effect of on ay is shown in Figure 2.4. The 
standard deviation of Y is 0.2001 when tx = 5, 0.2725 when tx is 7.5, and 0.3 when tx = 10. 

Y 



Figure 2.3: Effect of input X on output Y 2.4: Effect of tx on ay 

The sizes of the true linear and quadratic effects are 
True linear effect = 0.3 - 0.2001 = 0.0999 
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True quadratic effect = 0.2725 — (0.2001 + 0.3)/2 = 0.02245 


Estimation by Direct Observation 


One way of doing this would be to take ten samples when = 5 ,ten more samples when 
tx = 7.5, and a final set of ten samples when tx = 10 ■ For each group of ten samples, the 
standard deviations could be calculated and used to estimate the effect of tx on ay . One 
such data set is given in Tabic 2.2. 



Target of X 


5 

7.5 

10 


10.36 

13.40 

16.55 


10.07 

13.43 

16.91 


10.14 

13.17 

16.92 

ValiH's 

10.28 

13.28 

16.64 


9.71 

13.34 

16.67 


10.19 

13.30 

16.55 


10.05 

13.08 

16.68 


9.98 

13.13 

17.10 


10.37 

12.57 

17.04 


9.91 

12.78 

16.74 

Average 

10.11 

13.15 

16.78 

Std. Dev. 

0.208 

0.278 

0.199 


Table 2.2: Data for estimating effect of X on variation of Y 


The observed linear and (juadratic effects of adjusting tx are 

Observed linear effect = 0.199 - 0.208 = -0.009 

Observed (piadratic effect = 0.278 — (0.199 + 0.208)/2 = 0.0745 

A plot of standard deviation is shown in Figure 2.5 along with the true values. This demon- 
strates the inaccuracy involved in estimating standard deviation with small sample sizes. 

Estimation by Transmission of Variation 

The bestfit curve to the data of average in Table 2.2 is 

Y = 5.7755 + 0.6318X + 0M6X^ 


( 2 . 8 ) 
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The best fit curve is shown in Figure 2.6 along with the true curve. 




Figure 2.5: Observed standard deviations Figure 2.6: Observed averages 


Using the above equation the standard deviation of Y can be estimated as follows 

Uy = ^(0.6318 + 0.()92i^)%^ + 0.5(0.092)%^ (2.9) 

The standard deviation of X was found to be 0.1704. 

Evaluating the estimate of Oy gives 0.1875 at = 5 , 0.2274 at tx = 7.5, and 0.2674 at 

tx = 10 . 

Using these result the predicted linear and quadratic effects can be calculated as follows: 
predicted linear effect = 0.2674 - 0.1875 = 0.0799 
predicted quadratic effect = 0.2274 - (0.2674 + 0.18750)/2 = -0.00005 
The predicted effect of tx on ay is shown in Figure 2.7. 



Figure 2.7: Predicted standard deviations 

Once again there is lack of fit due to variation in the data and the fact that the true cubic 
polynomial is approximated by a quadratic polynomial. 
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By the above example we saw the limitation of both the method. In the practical case, when 
we are applying this method to reduce the variation, the bestfit surface through average 
and through standard deviation may not represent the relation between input and output 
correctly. In that case both the direct method and transmission of variation method are 


applied. 



Chapter 3 

Designing A Genetic Algorithm 


3.1 Introduction to Genetic Algorithm 

Genetic Algorithm (GA) is a search procedure based on the mechanics of natural selection 
and natural genetics. It uses the principle of survival of fittest same as nature. It abstracts 
operators from nature to form a mechanism that is suitable for a variety of search problems. 

In a simple GA, variables are coded in the form of strings containing alleles ( either 1 
or 0 ) similar to the chromosomes having genes in a biological system. In GA’s terminology, 
strings along with their decoded value are known as individuals. A set of such individuals in 
the user-set range is called a population. Like in the natural genetics, GA processes different 
population through generation. The three basic operators that modify the population to 
another are as follows: 

t 

• Reproduction 

• Crossover 

• Mutation 


The first operation employed in GA is reproduction or selection. In the population each 
individual is evaluated to find its fitness value. The fitness value is same as function value in 
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the case of maximization else in minimization it is In the selection process, individuals 
are given their presentation or copies according to their fitness value in order to form a set 
of good individuals called mating pool . 

One of the method of selection is the proportionate reproduction in which a string is 
selected for the mating pool with a probability proportional to its fitness value. 

Another method of selection is the tournament selection in which s individuals (with or 
without replacement) from the population are picked and best one is selected for the mating 
pool. In most GA applications, a binary tournament selection with s = 2 is used. 

The second operation employed in a GA is crossover. In a single point crossover, among 
the individuals in mating pool two are picked at random and a random cross site is selected 
across the string length. Alleles on one side are swapped between the individuals to form 
two different strings called the children points. For example, if si and S 2 are two individuals 
and a random cross site is selected as shown: 

si = 1 1 01 101 
s2=101|010 

After mating the two children points s[ and Sj will be 

si’ = 1 1 0 0 1 0 
s2’ = 1 0 1 1 0 1 

i 

In the uniform crossover bit by bit crossover is done with a probability 0.5. After mating 
the two children points s'l and S 2 ^^e 

sl> = 1 0 0 HI 0 1 

s2’ = 1 [T] 1 0 1 0 

Crossover is the only operator which creates new individuals by mixing the genetic informa- 
tion of parent individuals. This operation is done with a probability known as the crossover 



3.2 Application of Variation Reduction Technique for Designing GA 


20 


probability Pc- 

Mutation, the third operator, plays a secondary role in a GA. This is applied to get 
some features which cannot be achieved by crossover. Mutation is occasional alteration of 
an allele (changing 0 to 1 and vice versa) with some probability called mutation probability 
Pm- The action of mutation operator is shown below: 

10101-^ 1O101 

A flow chart of simple GA is shown in Figure 3.1. 

If proper GA parameters are not used, the performance of GA may change as the initial 
population is changed. In order to get consistent result, a variation reduction study is needed, 
which will give the approximate range of parameters for a consistent performance of GA. To 
quantify the effectiveness of different genetic algorithms, De Jong devised two performance 
measures, one to gauge convergence and other to gauge ongoing performance. He called 
these measure off-line (convergence) and on-line (ongoing) performance respectively. 

On-line performance is an average of all function evaluations upto and including current 
trial (De Jong, 1975). In a minimization problem if on-line performance is less that means 
that the performance of GA is good. 

3.2 Application of Variation Reduction Technique for 
Designing GA 

i 

The first step is to decide the performance measure. In the present study, the on-line 
performance of last 50 % generations is chosen as the performance measure. The various 
parameters that can affect on-line performance are : 
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Figure 3.1: Flow chart of simple Genetic Algorithm 
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1 String length 

2 Cardinality of alphabet i.e. binary, ternary, real coded etc. 

3 Linkage i.e. tight coding, loose coding, random coding. 

4 Population size 

5 Crossover probability 

6 Mutation probability 

7 Selection operator i.e. proportionate selection, ranking selection, tournament selection, 

genitor selection. 

8 Crossover operator i.e. single point crossover, two point crossover, uniform crossover. 

9 Mutation operator i.e. bitwise mutation, mutation clock. 

10 Termination criteria 

11 Sequence of GA operator 

If all the inputs will be considered for designing GA, the computational complexity is large. 
In the present study only a few inputs are fixed to reasonable values and others are varied. 
In the present study 

i 

• Binary GA is used. 

• String length is not taken as input parameter as affecting output, because string length 
mainly affects accuracy. 

• Tight coding is used as the linkage parameter. 

• Out of the various selection procedures tournament selection and proportionate (rowlett 
wheel) selection is chosen as input affecting on-line performance. 
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• Out of the various crossover operators single point crossover and uniform crossover is 
chosen as input affecting on-line performance. 

• Bitwise mutation is used as the mutation operator. 

• Ten thousand function evaluation is chosen as termination criteria. 

• llc'prod action then cro.ssovcr then mutation is chosen as sequence of GA operator. 

• Simple GA written in ’C’ i^rogramming language is used. 

Following are the parameters which is considered as candidate input variable: 

1 Population size (PS) in the range 0.25 x stringlength to 2.5 x stringlength 

2 Crossover probability (CP) in the range 0.5 to 1.0 

3 Mutation probability (MP) in the range 0.001 to 0.1 

4 Selection operator (SP) 

5 Crossover operator (CO) 


Out of the above j)aranieter key input variable, VAPs anti VIPs arc to bo found. Since 
selection procedure and crossover operator are discrete variable, handling this is found to be 
difficult when bestfit surface through average and standard deviation is found. So for each 
of the test function study, four cases are made and in each case these parameters is fixed. 

Case 1 Proportionate selection, single point crossover is used. 

Case 2 Proportionate selection, uniform crossover is used. 

Case 3 Tournament selection, single point crossover is used. 
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Case 4 Tournament selection, uniform crossover is used. 

I'\)r ('iu'li of t.lu' cnsc's population size, crossover probability and mutation probability are 
used as candidate input variable. 

In the present case following test functions is studied for optimization and variation reduction 
of on-line performance. 

1 DcJong FI : 

Minimize f{x) = J2i=i 

-5.12 < Xi < 5.12 

This is a three variable uuimodal function having a minima at (0,0,0)^. The GA 
will find the minimum point easily. The purpose of taking this function is to reduce 
the variation in the on-line performance of one of the simplest problem. That is why 
variation reduction is first tried on this function. 

2 DeJong F2 : 

Minimize f{x) = 100(a:f — + (1 ~ 

-2.048 < Xi < 2.048 

This is a two-dimensional minimization problem with the global minimum at (1, 1)'^. 
In this function, the slope towards the optimum point is very small and the population 
gets converged to a point along the ridge. The idea behind using this function is to 
apply variation reduction on slightly difficult function. 

3 Rastrigin’s function: 

Minimize f{x) = 200 -I- ~~ 10cos(2nxi)) 

-5.12 < Xi < 5.12 

This function has 20 variables and contains 2^° local minimum points, of which only 
one is the global minimum point. The global minimum point is at Xi = 0 for all 
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i = 1,2, ...,20. At this point the function value is zero. The other local minima 
occur when the cosine term is one. This function is a difficult optimization problem, 
because there are 2^° different local optimal points (surrounding the global optimum) 
with function values very close to that at the global optimal point. This function is 
used to show the variation reduction on one of difficult optimization problem. 

4 welded beam problem: 

For showing the variation reduction on one of engineering design problem, welded beam 
problem is selected. 

Minimize cost = \.\x'\x 2 + .0483;3.'e,i(L + x^) 
subjected to 

(;j:) = > 0 

52 (.x) = Ssy -r{x)>Q 
g,{x) = Pc{x) -F>Q 
5,1 (;r) = x.i -Xy >0 

The welded beam is shown in Figure 3.2. 

The parameter L is the length of the beam not welded and this is kept fixed to 14 
inches. The cost of fabrication includes material cost and the welding labour cost. The 
first constraint limits the bending stress {cr{x)) anywhere in the welded beam to the 
maximum allowable strength {Syt) of the beam material. The second limits the shear 
stress (r(.'r)) in the weld to the allowable shear stress (Ssy) of the material. The third 
constraint restricts the buckling of the plate [Pc{x)) due to applied load. The fourth 
constraint makes sure that the weld size (/i) is less than the thickness (6) of the beam. 
The details of different stresses and loads are given by: 

Bending stress, a{x) = 

Shear stress is given in terms of primary stress (r') and secondary stress (r") as: 
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Figure 3.2: Diagram of welded beam problem 


r(.x) = + 2r'r"cos(0) + (r")^ 

where 


\/2xiX2 ’ 

t( _ MR 
T — jr , 

M = F(L 4- f ) , 

J = 1.414xiX2(f| + (^)") 

The parameter M is the moment of F about the center of gravity and J is the polar 
moment of inertia of the weld group and cosO = ||. 


Buckling load is given by : 


ip ( r, A 4.Q13\/K/a7 l ^.fML 

I — j^2 2L V a 


where moment of inertia I = , a — ^Gxzx\ , 
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£■ = 30 X and shear modulus G = 12 x lO'^psz. The allowable shear strength 

taken is Syt = 30000 psi and maximum shear stress is Ssy = 13600psi . The load F is 
6000 lb. 



Chapter 4 
Results 


4.1 Test Function 1: De Jong Function FI 

4.1.1 Initial Capability Study 

In order to emphasize the need for a variation reduction study a capability study is first 
performed. In the capability study, twenty different trials are taken. For each trial all five 
variables are selected at random and ten runs are performed each time with a different initial 
population. The on-line performance measured for last 50 % of the generation are shown in 
Table 4.1. The average and standard deviation of ail 20 trials are also shown in the table. 



Run number { 

average 

mt 

j— 



pummi] 

HIBEIIil 






1 

30,99 

31.22 

32.57 

31.52 

31.91 

31.23 

31.13 

31.79 

31.69 

31.47 

31.55 

0.46 

2 

19.39 


17.80 

18.88 

19.51 

18.14 

21.50 

17.95 

19.87 

20.66 

19.18 

1.26 

3 

4.94 


00.59 

5.85 




00.55 



01.30 

2.17 

4 

23.47 

21.77 

23.43 

20.86 

21.40 

21.08 

21.99 

22.79 

24.12 

19.95 

22.09 

1.33 

6 

20.24 

11.47 

02.93 

01.93 


03.14 

21.62 

21.27 

23.50 

17.14 

12.93 

8.78 

6 

32.01 

33.72 

33.32 

32.45 

33.26 

33.45 

34.93 

36.09 

32.34 

33.15 

33.47 

1.23 

7 

24.64 

23.75 

24.25 

23.75 

24.31 

23.89 

24.15 

23.82 

23.52 

24.08 

24.01 

0.33 

8 

14.57 

15.39 

15.95 

14.38 

21.27 

21.39 

22.19 

20.61 

16.14 

20.29 

18.22 

3.17 

9 

66.84 

66.82 

67.18 

67.03 

65.85 

65.14 

66.62 

66.24 

66.02 

66.22 

66.30 

0.76* 

10 

09,70 

09.00 

00.02 

70.01 

00.00 

00.50 

70. 1 0 

70.24 

00.34 

00.75 

00.84 

0.27 

ll 

43.30 

42.76 

40.08 

42.78 

41.27 

41.4.3 

41.38 

41.30 

41.88 

40.72 

41.76 j 

0.90 

12 


02.06 

63.09 

00.73 

60.92 

62.20 

62.04 

01.14 

02.16 


61.90 1 

0.76 

13 

71.71 

71.64 

71.83 

72.21 

72.00 

71.94 

71.91 

71.98 

71.67 

72.17 

71,91 

0.20 

14 

20.55 

20.33 

23.14 

20.38 

20.92 

21.82 

20.47 

20.95 

20.32 

19.59 

20.85 

0.99 

15 

37.46 

36.39 

36.35 

37.52 

37.97 

37.34 

38.16 

37.33 

37.26 

35.80 

37.16 

0.75 

16 

53.81 

50.44 

50.41 

52.19 

50.20 

53.20 

52.47 

51.73 

51.59 

51.66 

51.77 

1.20 

17 

22.31 

21.10 

22.97 

24.95 

22.26 

23.67 

22.21 

22.48 

23.65 

24.57 

23.02 

1.19 

18 

10.92 

17.49 

10.41 

06.38 

05.46 

06.89 


06.84 

11.53 

■IKIMM 

08.92 

3.93 

19 

34.44 

33.43 

35.09 

34.65 

34.45 

36.81 

34.73 


33.86 

37.59 

34.81 

1.42 

20 

37.28 

35.60 


36.09 

32.83 

34.31 

35.58 

i 35.08 

36.38 

33.97 

34.96 

1.55 


Table 4.1: Initial capability study for function FI 
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The standard deviation of each trial is the measure of variation within the trial due to change 
in initial population. This table is used to construct the average and standard deviation 


control charts as shown in Figure 4.1 and Figure 4.2. 



Figure 4.1: Control chart for average on-line Figure 4.2: Control chart for standard devi- 
performance before variation reduction for ation of on-line performance before variation 
function FI reduction for function FI 


If Xi be the average of every trial and Si be the standard deviation of every trial, then the 
grand average is given as follows: 


t^20 v. 


( 4 . 1 ) 


and the standard deviation is given as follows: 
5 = 




20 


( 4 . 2 ) 


Accepting acceptable trials within ±3cr from mean, the upper control limit of the average is 
3S 


UCL^ = X + 


n 


( 4 . 3 ) 


and tlie lower control limit of the average is 

LCL^ = X-^=. (4-4) 

r/n ■ . 

Similarly, the upper control limit of the standard deviation is given as 


UCL, = B^S 


( 4 . 5 ) 
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and the lower control limit of the standard deviation is given by 


LCLs = Br,S 


(4.6) 


where Dq = 1.CG9 and — 0.2759 for number of 10 runs (Taylor, 1991). 
For the data in Table 4.1, we compute the above quantity; 



= 34.297 

s 

= 2.483 

UCLx 

= 36.65 

LCLx 

= 31.94 

UCLs 

= 4.145 

LCLs 

= 0.6852 


The average control chart is a plot of the averages of all twenty trials. The control limit 
represents the extent that the group average should vary if the input variables are changed. 
The standard deviation control chart is a plot of the standard deviation of all twenty trials. 
Figure 4.1 shows the plot of average of on-line performance of all 20 trials. The figure clearly 
shows that points of the trials are not within acceptable region. Similarly Figure 4.2 shows 
the standard deviation of all 20 trials, of which 5 trials are outside acceptable limits. 

Thus, the GA considered above is not a capable process. In order to make the GA 
search process cai)able, we need to use the variation reduction study as described in the 
previous chapter. 

Among the five parameters, the selection operator and the crossover operator are con- 
sidered only of two types each. Thus, we consider four different cases each with a particular 
combination of selection and crossover operator and the effect of three parameters — popu- 
lation size, crossover probability and mutation probability are studied. 
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4.1.2 CASE 1 ( Proportionate Selection and Single Point Crossover) 


There are only three parameters to study in this case. Thus, 27 trial experiment is performed 
as described in chapter 2. For each trial, with proportionate selection as selection operator 
and the single point crossover as the crossover operator and by choosing the combination 
of three parameter from Table 4.2, GA runs are performed. For each trial, ten runs are 
performed by changing the initial population. The average and the standard deviation of 
these 10 runs are calculated for each trial and are tabulated in Table 4.2. 


T'rial 

no. 

candidate input 
variables 

Run number 

on-line performance 



MP 

r 


3 



er 


g- 


nr 

average 

8td. dev. 

1 

112 

0.50 

0.100 

24.9 

25.6 

24.9 

25.4 

24.8 

24.9 

24.0 

24.9 

25.3 

26.5 

25.02 

0.46 

2 

112 

1.00 

0.001 

13.3 

13.1 

13.1 

8.3 

17.6 

15.1 

13.4 

14.0 

10.3 

11.8 

13.00 

2.54 

3 

11 

0.50 

0.100 

15.2 

13.6 

12.9 

12.4 

13.9 

14.3 

13.8 

12.4 

12,6 

12.6 

13.36 

0.94 

4 

112 

0.50 

0.001 

14.2 

15.8 

14.9 

7.4 

11,3 

17.0 

15.5 

17.7 

16.0 

21.0 

15.08 

3.67 

5 

' 11 

1 .00 

0.001 

2.0 

l.l 

1.4 

1.0 

1.3 

1.3 

0.8 

0,1 

0.0 

0.0 

0.90 

0.68 

0 

1 1 

l.OO 

0.100 

13.9 

17.5 

12.8 

13.0 

13.5 

15.8 

12.7 

12.6 

12.4 

12.6 

13.68 

1.67 

7 

112 

1.00 

0.100 

25.0 

24.6 

24.7 

24,8 

24.5 

24.0 

25.0 

24.8 

24.1 

25.3 

24.68 

0.42 

8 

11 

0.50 

0.001 

21.0 

17.7 

18.1 

0.4 

0.0 

0.0 

O.I 

0.0 

O.I 

0.1 

5.75 

9.15 

9 

61 

0.75 

0.051 

20,5 

23.4 

19.9 

20.6 

21.0 

21.3 

22,2 

20.8 

20.1 

20.4 

21.02 

1.05 

10 

61 

0.50 

O.OOl 

2.4 

4.8 

3.1 

1.5 

0.7 

1.2 

2.2 

5.5 

1.5 

10.4 

3.31 

2.93 

11 

Cl 

0.50 

0.100 

22.8 

22.7 

23.1 

23.7 

23.3 

23.1 

22.6 

22.7 

23.6 

23.4 

23.08 

0.40 

12 

61 

1.00 

0.001 

4.4 

1.5 

1.6 

2.0 

7.7 

0.2 

5.3 

3.4 

4.9 

0.9 

3.18 

2.37 

13 

61 

1.00 

0.100 

23.7 

23.6 

23.5 

23.4 

23.5 

22.9 

23.0 

23.2 

23.3 

23.0 

23.31 

0.27 

14 

11 

0.75 

0.001 

8.2 

0.3 

0.1 

0.1 

0.1 

1.2 

1.1 

0.9 

0.7 

0.4 

1.30 

2.47 

15 

11 

0.75 

0.100 

13.7 

14.0 

18.8 

15.6 

15.8 

13.2 

15.9 

13.4 

15.4 

18.6 

15.43 

1.99 

16 

112 

0.75 

0.001 

10.7 

19.0 1 

11.9 

14.9 

23.9 

18.0 

14.3 

15.2 

22.2 

10.6 

16.06 

4.62 

17 

112 

0.75 

0.100 ! 

24.9 

25.1 i 

25,1 

25.2 

25.3 

24.7 

24.4 

24.7 

25.1 

24.5 

24.90 

0.30 

18 

11 

0.50 

0.051 

13.7 

20.3 

19.2 

18.9 

19.8 

18.9 

19.0 

19.4 

18.5 

18.4 

18.60 

1.82 

19 

11 

1.00 

0.051 i 

9.4 

7.9 

7.6 

8.2 

7.3 

7.6 1 

12.0 

10.0 

11.6 

10.2 

9.20 

1.72 

20 

112 

0.50 i 

0.051 

23.2 

22.8 

23.1 

23.3 

24.4 

24.5 i 

22.9 

24,3 

23.9 

23.7 

23.63 

0.65 

21 

112 

1.00 1 

0.051 

23.1 

24.2 

22.3 

24.8 

24.0 

24.5 ! 

23.2 

23.7 

24.6 

25.0 

23.95 

0.85 

22 

61 

0.75 ! 

0.001 

5.1 

0.4 

0.3 

3.4 

4.4 

1.4 ! 

2.6 

5.7 

1.8 

0.3 

2.53 

2.03 

23 

61 i 

0.75 

0.100 

23.1 

23.4 

22.7 

22.9 

23.0 

23.5 

22.6 j 

23.2 

23.2 

23.6 I 

23.21 

0.42 

24 

1 1 

0.75 

0.051 

10. 1 

14.3 

0 8 

10,0 

14.0 

15.1 1 

10.2 

8.5 

7.4 

8.8 i 

10.01 

2.67 

2 ^ 

112 

0,75 

0.051 

24.0 

25,3 

2:1.0 

21.0 : 

24 3 

23.;! 

24.1 

24.8 

24.7 

23.5 j 

24.17 

0.75 

20 

61 

0.50 

0.051 

21.2 

20.7 

22.0 

21.1 

20.8 

23.4 j 

18.0 1 

21.8 

22.2 

22.0 j 

21,50 

1.28 

27 

01 

1.00 

0.051 

21.2 

21.6 

19.3 

20.7 

22.6 

21.3 

21.9 I 

20.5 

22.4 

20.9 

21.25 

0.97 


Table 4.2: Data obtained from 27 trials on DeJong FI function case 1 

The effect of each parameter and their interaction with each other is calculated. The 
quadratic effect of each parameter is also calculated by the same way as described in chapter 
2 . 

The significant effects of parameters, their interactions and quadratic effect are shown 
in Table 4.3 for the average. The table shows that PS, CP and MP affect the average on-line 
performance. The significant effects are linear in PS, CP, MP and PS*CP, and quadratic in 
MP. These linear effects of PS, CP and MP represent the amount that average of on-line 
performance of the trial changes as a result of adjusting the parameter from their low value 
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to high value. 


Design parameter 

Effect 

MP 

mutation prob. 

13.951 

PS 

population size 

11.26 

MP*MP 

mutation prob. (quad) 

5.593 

PS*CP 

population size * crossover prob. 

1.9734 

CP 

crossover prob. ; 

-1.796 


'lablo 4.3; Effects table for average case 1 


The effect that design parameters have on standard deviation is shown in Table 4.4 . Since 
MP and CP*MP have significant effect on standard deviation means MP and CP are VAPs. 


Design parameter 

Effect 

MP 

CP*MP 

mutation prob. 

crossover prob.* mutation prob. 

2.62 

1.7864 


Table 4.4: Effects table for standard deviation case 1 


Response surface study 

We fit the best surface on the average of on-line performance from data presented in Table 
4.2. The resulting equation is : 

On-line performance = 6.8046 -f 0.059PS' -I- 226. lAfP — 7.7CP — 873. 

+0.069PP X CP (4.7) 

This equation can be used to [)crform the variation transmission analysis. As part of 
variation transmis.sion analysis the tolerances for the design parameters can be established 
for a capable GA. 

The best fit surface through the standard deviation of on-line performance presented 
in Table 4.2, is also found for doing analysis by direct observation. The surface through 
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logatithin (due to statistical leasoii) of standard deviation is piissed. The resulting equation 
is: 


log((7) = 1.09 - 17.87MP + 4.66CP x MP (4.8) 

Analysis by direct observation 

We see that MP and CP are VAPs as their effects are significant on the standard deviation 
of the on-line performance. The minimization of the equation 4.8 gives 

a = 0.G28 at CP = 0.5 and MP = 0.1. 

These values of CP and MP can be put into equation 4.7 of on-line performance and 
on-line performance is minimized, keeping PS as the only variable. Since, in this case there 
is only one variable, we can minimize the above equation by hand calculation. But in case, 
there are more variables, GA is used for minimization. 

Minimum on-line performance = 17.8568 at PS = 11 

The direct observation can’t find VIPs. Based on these analysis, we conclude the 
following: 

• Set CP = 0.5 and MP = 0.1 in order to minimize variation. 

• Set PS = 11 to reduce the average. 

Analysis by Variation transmission 

Variation transmission analysis requires an equation such as one obtained for the on-line 
performance. It assumes that the average of the process parameters are tps, tcp-, Imp and 
that these parameters vary around their averages with standard deviations aps,<Jcp,crMP 
respectively. Based on these assumptions, the standard deviation of on-line performance can 
be determined as follows: 

= (0.059 + Q.Omcpfals + (-7-7 + 0M9tpsfcrlp + 
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(226.1 - 2 X 873.63tA/p)V^^P + 0.5 x (2 x 873.63) + 0.069V^5 Ct|p (4.9) 

VVe see that the terms tp 5 , tcp, t^p are present in the above equation of standard 
deviation. This means that PS, CP and MP all are VAPs. In the above equation, the 
following values of standard deviation in the parameters are assumed. 

(T i>s ~ 5, ^cp ~ 0.05, (^MP ~ 0.005. 

Putting these values in the equation 4.9 and minimizing the equation, yields; 
a = 0.53 at PS = 105, CP = 0.5, MP = 0.1 

The outcomes of both methods are not the same because the analysis by direct obser- 
vation shows that only CP and MP are VAPs, whereas analysis by variation transmission 
shows that PS, CP and MP are VAPs. 

The standard deviation from bestfit surface will be 0.628, at PS = 105, CP =0.5 and 
MP=0.1. Both methods showed that at CP = 0.5 and MP =0.1, the standard deviation 
will be minimum. Keeping PS = 11, in the equation 4.7 will minimize the average of on-line 
performauce. 

So with mean at PS = 11, CP = 0.5 and MP = 0.1 and standard deviation aps = 
5, Gcp = 0.05, and a^p = 0.005, final capability study is done to verify the expected 
improvements. 

4.1.2. 1 Final Capability Study for Case 1 

The oil-line performance measured for last 50 % of the generation are shown in Table 4.o, 
for 20 trials. The control chart for average and standard deviation is also given in Figure 4.3 
and Figure 4.4 . The process seems to be capable except two trials falling slightly outside 
LCL and UCL of average. We now try the next case. 
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Trial 

number 


Run number | 

average 

std . 

dev . 

1 

2 

3 


5 



8 

9 

10 

1 

16.01 

15.15 

14.51 

16.76 

15.53 

17.41 

15.89 

18.01 

18.13 

14 . 5 ^”"' 

16.19 

1.34 

2 

17.77 

15.70 

22.40 

21.03 

16.31 

20.24 

16.61 

17.06 

17.82 

22.11 

18.71 

2.51 

3 

20.14 

19.44 

19.87 

14.25 

18.39 

14.65 

13.56 

17.72 

16.13 

20.32 

17.45 

2.61 

4 

16.73 

14.92 

16.97 

14.08 

21.68 

17.16 

17.78 

16.24 

14.62 

19.84 

17.00 

2.35 

5 

13.88 

17.67 

16.88 

19.60 

21.06 

21.29 

20.74 

18.10 

16.01 

‘ 15.34 

18.06 

2.57 

6 

18.90 

18.58 

19.55 

19.21 

17.93 

18.22 

18.27 

20.48 

18.53 

18.50 

18.82 

0.76 

7 

19.77 

21.11 

21.28 

17.45 

20.71 

17.27 

16.69 

20.34 

20.80 

20.73 

19.61 

1.77 

8 

14.82 

14.29 

19.86 

16.84 

18.82 

20.35 

12.97 

18.46 

16.20 

14.76 

16.74 

2.54 

9 

22.22 

19.40 

19.31 

20,17 

20.19 

19.15 

20.06 

19.88 

22.62 

19.53 

20.26 

1.20 

10 

18,30 

14.03 

17.05 

20,25 

13.99 

19.57 

18.01 

14.82 

17.63 

17.14 

17,08 

2.18 

1 1 

18,99 

19.05 

19.71 

18.85 

20.08 

19,04 

20.21 

20.89 

17.95 

18.52 

19.34 

0.87 

12 

17.69 

19.50 

21.73 

17.99 

16,93 

17.73 

17.88 

17.24 

17.12 

17.15 

18.11 

1.48 

13 

14.46 

12.34 

13.06 

12.49 ! 

12.83 

18.31 

17.19 

14.20 

12.18 ; 

14.80 

14.19 

2.10 

14 

13.84 

16.35 

16.83 

21.38 

13.43 1 

20.57 

16.56 

15.46 

17.74 i 

15.45 i 

10.76 

2.58 

15 

15.74 

15.81 

17.45 

19.70 i 

13.67 1 

21.21 

20.74 

17.77 

21.63 

15.69 ! 

17.94 

2.75 

16 

19.36 

22.33 

19.49 

19.22 I 

18.34 i 

20.24 

19.75 

17.97 

19.41 

17.76 ; 

19,39 

1.31 

17 

17.00 

20.35 

20.97 ' 

17.30 

21.02 

18.65 

18.88 

20.89 

19.13 

20.34 

19.45 : 

1.50 

18 

15.69 1 

13.97 

17.15 i 

19.67 

22.29 1 

21.25 

13.20 

18.57 

15.19 

13.49 

17.05 * 

3.27 

19 

19.84 1 

22.52 

21.37 

20.07 

22.72 

21.61 

20.18 

20.00 

19.47 

18.68 

20.65 

1.34 

20 

18.72 

15.07 

18.73 

20.11 

13.82 

14.15 

18.04 

14.85 

13.83 

16.87 

16.42 

2.36 


Tabic 4,5: 


Final capability study for case 1 function FI 




Figure 4.3: Control chart for average of on- 
line perforniance after variation reduction 
for case 1 on function FI 


Figure 4.4: Control chart for starulard devi- 
ation of on-line performance after variation 
reduction for case 1 on function FI 
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4.1.3 CASE 2 (Proportionate Selection and Uniform Crossover) 

The same way as in case 1, 27 trials of experiment are done with proportionate selection as 
the selection operator and uniform crossover as the crossover operator. For each trial, the 
combination of variables are taken from Table 4.6 and ten runs are taken for each trial by 
changing the initial population. The average and standard deviation are calculated for each 
trial. The data are shown in Table 4.6. 


TViftl 

no . 

candidnte input 
variables 

Rtitt number 

on-line performance 

PS 

UP 

WP 

1 

T 






g - 


- nr - 



1 

112 

0.50 

0.100 

24.7 

25.0 

25.5 

25.1 

26.0 

24.7 

25.0 

24.8 

HEIB 

mfzmm 

25.01 

0.43 

2 

112 

1.00 

0.001 

6.9 

5.3 

21.8 

17.7 

20.0 

16.0 

14.4 

19.5 

19.5 

12.4 

15.34 

5.63 

3 

11 

0.50 

0.100 

14.4 

16.6 

17.8 

17.8 

15.0 

12.9 

13.3 

13.8 

14.1 

13.2 

14.90 

1.87 

A 

112 

0.50 

0.001 

10.0 

18.9 

15.4 

16.9 

17.5 

18.1 

22.9 

18.3 

15.1 

15.6 

16.87 

3.31 

5 

11 

1.00 

0.001 

1.4 

0.0 

1.0 

0.8 

1.7 

0.1 

0.0 

0.1 

0.0 

0.0 

0.52 

0.66 

0 

11 

1.00 

0.100 

15.4 

16.2 

17.5 

14.3 

14.2 

15.3 

14.6 

16.3 

13.8 

14.4 

15.21 

1.18 

7 

112 

1.00 

0.100 

25.3 

24.8 

25.6 

25.5 

25.4 

25.9 

25.7 

25.4 

25.3 

25.2 

25.39 

0.29 

a 

11 

0.50 

0.001 

12.9 

6.1 

7.9 

7.4 

8.1 

7.1 

7.1 

0.9 

7.1 

8.0 

7.86 

1.89 

9 

61 

0.75 

0.051 

23.8 

24.2 

22.9 

23.0 

23.9 

23.4 

23.3 

24.4 

23.5 

22.7 

23.51 

0.56 

10 

01 

0.50 

0.001 

9.1 

10.4 

1.6 

2.3 

11.1 

2.0 

12.0 

12.9 

7.0 

0.9 

0.92 

4.79 

1 1 

61 

0.50 

0.100 

24.1 

23.7 

23.9 

23.6 

24.3 

24.4 

23.7 

23.4 

24.1 

24.3 

23.93 

0.35 

12 

61 

1.00 

0.001 

8.6 

2.4 

2.5 

0.3 

1.2 

3.2 

13.2 

2.8 

1.4 

11.3 

4.68 

4.59 

13 

61 

1.00 

0.100 

24.4 

24.2 

24.0 

24.0 

24.7 

24.3 

24.4 

24.5 

24.5 

24.1 

24.33 

0.22 

14 

11 

0.75 

0.001 

0.0 

0.0 

0.4 

0.2 

1.7 

0.5 

0.4 

0.8 

0.7 

0.8 

0.56 

0.49 

15 

11 

0.75 

0.100 

15.7 

15.8 

13.5 

13.9 

14.3 

13.3 

15.3 

14.5 

14.2 

15.7 

14.63 

0.93 

16 

112 

0.75 i 

0.001 

24.1 ' 

25.7 

22.0 ' 

19.0 

19.0 

18.6 

8.6 

8.8 

12.3 

15.0 

17.31 

5.99 

17 

112 : 

0.75 

0.100 

25.0 1 

25.1 

25.0 

24.7 ! 

25.7 

25.8 

25.6 

24.7 

25.1 

24.7 

25.12 

0.42 

18 

11 

0.50 

0.051 

16.4 

14.0 

18.9 

16.9 i 

16.6 

9.2 

8.2 

8.8 

8.4 

0.0 

12.64 

4.31 

19 

11 

1.00 

0.051 

11.3 ' 

10.8 

18.2 

21.2 

23.0 

22.1 

22.0 

21.2 ' 

22.3 

23.3 

19.54 

4.69 

20 

112 

0.50 

0.051 

25.9 

25.8 

25.2 

24.5 

25.1 

24.9 

24.7 

24.6 

24.6 

23.9 

24.92 

0.60 

21 

112 

1.00 

0.051 

25.3 

24.8 1 

25.3 

24.0 

24.9 1 

24.4 

24.8 

25.4 

25.1 

25.2 

24.93 

0.46 

22 

61 

0.75 

0.001 

3.1 

3.6 

7.7 

0.8 

3.9 

0.4 

10.2 

10.7 

15.2 

12.0 

6.75 

6.10 

23 

61 

0.75 

0.100 

23.8 

24.3 

23.9 

24.0 

24.7 

24.1 

23.8 

23.8 

24.7 

23.5 

24.04 

0.41 

24 

11 

0.75 

0.051 

10.0 

11.6 

23.5 

19.7 

20.4 

19.7 ' 

20.0 

20.8 

21.0 

19.9 

18.66 

4.29 

25 1 

112 

0.75 

0.051 

24.8 

25.3 

24.9 

25.1 

24.9 

24.6 

24.4 

25.3 

24.2 

25.1 

24.87 

0.38 

26 

61 

0.50 

0.051 

23.9 

23.2 

21.4 

23.0 

23.9 

23.2 

22.6 

21.7 

23.0 

23.3 

22.93 

0.83 

27 1 

61 

1.00 

0.051 

23.1 

24.1 

22.9 

23.5 

23.7 

23.8 

22.7 

23.4 

24.3 

22.3 

23.37 

0.64 


Table 4.6: Data obtained from 27 trials on DeJong FI function case 2 

The effect of each parameter and their interaction with each other is calculated. The 
quadratic effect of each parameter is also calculated. The significant effects of parame- 
ters, their interaction and quadratic effects are shown in Table 4.7 for the average. PS and 
MP arc significant parameters affecting the average. In addition, the interaction PS*MP, 
CP*MP and the (luadratic effect of MP arc also found to be significant. 

The effect that design parameters have on standard deviation is shown in Table 4.8 . Since 
MP and PS*MP have significant effects on the variation, this means MP and PS are VAPs. 
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Design parameter 

Effect 

MP 

mutation prob. 

12.86 

PS 

population size 

10.58 

MP*MP 

mutation prob. (quad) 

6.739 

CP*MP 

crossover prob. * mutation prob. 

2.03 

PS*MP 

population size * mutation prob. 

-1.63 


Table 4.7; Effects table for average case 2 


Design parameter 

Effect 

MP 

l’S*MP 

mutation prob. 

population size * inutaiion prob. 

-2.92 

-2.45 


Table 4.8: Effects table for standard deviation CJisc 2 


Response surface study 

Based on 27 trials data in Table 4.6, the best surface through the average is fitted.. The 
resulting equation is : 

on-line performance = 2.4936 + 0.1159PS' + 21A.9MP — 1313.87AfP^ 

+77.494C'P X MP - .2206P5' x MP (4.10) 

The best fit surface through standard deviation is also found using the data of Table 
4.6. The resulting equation is ; 

log(a) = 0.968 - 6.1429MP- 0.1597P5 X MP (4-11) 

Analysis by direct observation 

Wc see that MP and PS are VAPs as its effect are significant on the standard deviation of 
the on-line performance. The minimization of the equation 4.11 for the standard deviation 

gives 


a = 0.238 at PS = 112 and MP = 0.1 
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These values of PS and MP can be put into equation 4.10 of average of on-line perfor- 
mance and on-line performance is minimized. 

Minimum on-line performance = 25.229 at CP = 0.5 

Based on these analysis we concluded the following: 

• Set PS = 112 and MP = 0.1 in order to minimize variation. 

• Set CP = 0.5 to reduce the average. 

Analysis by Variation transmission 

The equation for standard deviation of on-line performance can be found by equation 4.10 
of average of on-line performance. 

= {.im- .2206tMpf<yj,s + {77.mtMp)^(ylp 

-|-(214.9 - 2 X 1313.87iAfp + 77.494«cp - .2206ip5 ) V^p + 

0.5 X (2 X 1113.87)Vl^p -b (77.494)VcpaMP + 0.2206V|,5a^p (4.12) 

Putting aps = 5,acp = 0.05, aMP = 0.005 in the equation 4.12 and minimizing the 
equation yields; 

(7 = 0.5928 at PS = 105, CP = 0.54, MP = 0.0818. 

The standard deviation from eciuation 4.11 will be 0.4 at PS = 105, CP = 0.54 and 
MP = 0.0818. Comparing the standard deviation found from equation 4.11 from both the 
method, shows that the the results from the variation transmission method is not bad. By the 
equation obtained from standard deviation of on-line performance, the transmitted variation 
can be calculated by the method explained in chapter 2. The transmitted variation are : 

CtPS ~ 0.489, Ctcp ~ 0.3175, ctMP ~ 0.1066, 

We see that the transmitted variation by PS and CP is more, so there is need to reduce 
the variation given to PS and CP as input. Hence PS and CP are VIPs. So with mean at 
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PS — 105, CP — 0.54 and MP = 0.0818 and standard deviation ops = 5, ccp — 0.05, and 
<^MF — 0.005 final capability study is done to verify the expected improvements. 

4. 1.3.1 Final Capability Study for Case 2 

The data are given in Table 4.9. The control chart for average and standard deviation is 
also given in Figure 4.5 and Figure 4.G . The proce.ss seems to be capable as all points are 
within UCL and LCL of average and standard deviation. Also the average and the standard 
deviation for each trial is reduced and is consistent. We now try the next case. 


Trial 

number 

Run number j 

average 

stci. 

dev. 

1 

2 



5 

6 

7 

8 



1 

25.05 

24.81 

24.35 

25.03 

25.00 

24.93 

24.84 

24 rod 

25.70 

25.24 

24.99 

0.34 

2 

24.09 

24.44 

24.91 

25.05 

24.00 

24.09 

25.16 

25.09 

25.24 

24.95 

24.85 

0.35 

3 

24.04 

24.30 

25.43 

24.92 

25.31 

23.91 

24.00 

24.87 

24.89 

23.63 

24.60 

0.00 

4 

24.38 

24.00 

25.20 

25.24 

25.03 

24.52 

24.64 

24.65 

2.5.03 

25.16 

24.80 

0.48 

5 

25.52 

24.70 

25.04 

24.42 

24.01 

23.73 

23.52 

24.48 

24.58 

24.49 

24.52 

0.57 

6 

24.42 

24.54 

25.12 

24.69 

24.27 

25.48 

24.07 

25.15 

24.03 

24.81 

24.66 

0.49 

7 

23.86 

25.02 

24.61 

25.21 

24.79 

24.51 

24.99 

24.52 

24.93 

25.41 

24.79 

0.44 

8 

24.72 

24.89 

24.17 

24.76 

24.72 

24.76 

24.05 

24.78 

24,80 

24.82 

24.65 

0.29 

9 

25.01 

24.86 

25.23 

25.12 

25.71 

24.76 

25.25 

24.85 

25.37 

24.33 

25.05 

0.38 

10 

24.23 

24.99 

24.51 

24.96 

24.69 

25.25 

25.59 

24.88 

24.81 

24.62 

24.85 

0.38 

11 1 

25.41 

25.36 

24.99 

25.40 

25.33 

25.22 

24.82 

25.27 

25.00 

24.82 

25.16 

0.23 

12 

24.86 

24.80 

25.48 

25.46 1 

25.30 

24.82 

25.63 

25.19 

25.05 

. 25.55 

25.21 

0.32 

13 

25.52 

25.20 

24.90 

24.73 ; 

24.33 

24.33 

25.62 

25.06 

24.88 

24.51 

24.91 

0.45 

14 

24.77 

24.76 

24.70 

24.70 

24.87 

23.95 

24.05 

24.67 

25.02 

24.33 

24.58 

0.35 

15 

24.94 

24.39 

24.68 

24.78 

24.87 

25.14 1 

24.79 

25.45 

25.12 

24.40 

24.86 i 

0.33 

Id 

24.02 

24.93 

25.25 

24.87 

24.96 

25.05 

23.92 

24.53 

24.00 

25.08 

24.84 

0.47 

17 

24.95 

24.82 

24.39 

24.80 

25.40 

24.04 

24.87 

24.42 1 

24.75 

25.03 

24.75 

0.38 

18 

24.00 

23.81 

24.45 

24.54 

25.28 

24.75 

24.19 

24.46 

25.14 

25.24 

24.65 

0.47 

19 

25.65 

24.87 

24.75 

24.51 

24.96 

24.44 

24.53 

24.89 

24.83 

24.66 

24.81 

0.35 

20 

24.48 

25.82 

25.20 

25.23 

24.99 

25.33 

24.98 

25.18 

24.58 

24.70 

25.05 

0.40 


Table 4.9: Final capability study for case 2 function FI 




Figure 4.5; Control chart for average of on- Figure 4.6: Control chart for standard devi- 
line performance after variation reduction ation of on-line performance after variation 
case 2 on function FI reduction case 2 on function FI 





4.1 Test Function 1: De Jong Function FI 


40 


4.1.4 CASE 3 (Tournament Selection and Single Point Crossover) 

As in case 2, 27 trials of experiment are done with tournament selection as selection operator 
and single point crossover as crossover operator. For each trial the variable combination is 
taken from Table 4.10. For each trial, ten runs are taken by changing the initial population. 
The average and standard deviation are calculated for each trial. The data are shown in 
Table 4. 10. 


'Trial 

no. 

cai 

(lid ate 
variabli 

nput 


Run nutnber 

On-line performance 


Df> 


Ml’ 

1 

2 

3 

4 


r~' 

^ 

jpi 

jy-n 

^ 

average 

■ td. dev, 

1 

1 12 

6.50 

0,100 

38.4 

37.8 

37.5 

37.6 

38.1 

38.8 

38.2 

38.3 

37.7 

38.1 

38.05 

(TTr” 

2 

1 12 

1.00 

0,001 

77.5 

77.3 

77.6 

76.1 

77.5 

77.2 

77.7 

77.7 

76.1 

77.1 

77.20 

0.61 

3 

1 1 

0.60 

0.100 

25,2 

22.6 

23.2 

24.6 



24.7 

22.7 

22.8 

25.5 

23,87 

1.08 

4 

112 

0.50 

0.001 

77.6 

77.0 

76.9 

76.2 

76.5 

77.1 

74.0 

77.6 

75.9 

76.9 

76.55 

1.07 

5 

1 1 

1.00 

0.001 


48.4 

55.2 

68.9 

53.0 

65.9 

54.3 

50.6 

49.1 

68.7 

58.32 

8.75 

0 

1 1 

1.00 

0.100 

25.5 

24.6 

24.5 

24.4 

22.6 

23.6 

24.3 

25.3 

23.2 

28.1 

24.62 

L53 

7 

112 

1.00 

0.100 

37.4 

36.9 

37.5 

38.2 

38.4 

37.2 

37.7 

38.5 

37.3 

38.3 

37.74 

0.67 

R 

11 

0.50 

0.001 

69.3 

57.4 

58.3 

G0.O 

56.8 

58.8 

50.5 

51.9 

68.9 

69.5 

61.03 

7.48 

9 

01 

0.75 

0.051 

51.5 

51.3 

53.0 

51.2 

51.0 

51.2 

52.6 

54.4 

53.0 

54.1 

62.33 

1.27 

10 

61 

0.50 

0.001 

76.6 

75.7 

76.8 

76.8 

75.9 

76.2 

76.3 

70.6 

76.8 

76.5 

76.41 

0.39 

1 1 

01 

0.60 

0.100 

30.3 

37.1 

35.7 

36.6 

35.4 

38.3 

35.7 

38.4 

37.4 

38.2 

36.91 

1.16 

12 

01 

1.00 

0.001 

76.6 

76.5 

76.9 

76.8 

76.0 

76.8 

70.4 

76.6 

76.5 

76.6 

76.47 

0.43 

13 

01 

1.00 

0.100 

38.1 

37.1 

34.6 

37.4 

36.0 

38.5 

36.1 

37.2 

38,0 

36.5 

36.93 

1.17 

H 

11 

0.75 

0.001 

59,4 

68.5 

57.6 

68.8 

49.9 

69.0 

57.3 

49.3 

69.2 

51.7 

59.06 

7.71 

15 

11 

0.75 

0.100 

29.2 

30.7 

29,2 

29.3 

33.3 

29.8 

29.4 

28.8 

29.1 

30.6 

29.96 

1.35 

16 

112 

0.75 

0.001 

77.7 

77.1 

77.7 

76.7 

77.4 

77.6 

77.2 

77.6 

77.2 

77.7 

77-40 

0.32 

17 

112 

0.75 

0.100 

38.6 

38.4 

38.2 

38.5 

38.1 

37.8 

37.6 

39.3 

37.7 

38.0 

38.22 

0.51 

18 

n 

0.50 

0.051 

29.8 

30.4 

29.4 

28.7 

28.3 

29.3 

28.1 

28,1 

47,4 

38.0 

31.76 

6.24 

19 

11 

1.00 

0.051 

29.4 

29.6 

28.4 

29.2 

28.3 

28.6 

29.5 

31.1 

32,6 

30.7 

29.73 

1.37 

20 

112 

0.50 

0.051 

64.0 

53.2 

52.1 

54,7 

54.0 

53.7 

54.0 

63.4 

53.3 

62.5 

53.48 

0.79 

21 

112 

1.00 

0.061 

64.5 

54.5 

62.8 

54.3 

53.8 

52,3 

63.1 

54.0 

62.3 

53.3 

53.49 

0.86 

22 

Gl 

0.75 

0.001 

76.6 

76.8 

77.0 

76.0 

75.9 

76.8 

76.6 

76.6 

76.1 

76.6 

76.39 

0.67 

23 

61 

0.75 

0.100 

35.9 

37.1 

39.6 

37.0 

36,6 

37.7 

35.9 

38.2 

36.6 

36.9 

37.16 

l.H 

24 

11 

0.75 

0.051 

30.1 

29.2 

28.1 

27.6 

30.7 

29,5 

27,0 

30.9 

29.8 

29.8 

29.27 

1.31 

25 

112 


0.051 

53.0 

52.7 

54.1 

54.6 



54.6 

54.7 

54.5 

53.6 

63.75 

0.91 

26 

61 


0.051 

52.6 

63.0 

53.1 

53.2 

50.1 

52.8 

52.8 

54.1 

52.3 

62,0 



27 

01 


0.051 

54.2 

54.4 

64.0 

53.6 

52.3 

53.3 

53.2 

60.4 

53.4 

53.2 

53.20 

1.15 1 


Table 4.10: Data obtained from 27 trials on DeJong FI function case 3 


The significant effects of parameters, their interaction and quadratic effect is shown in Table 
4.11 for the average, and in Table 4.12 for the standard deviation. 


Design parameter 

Effect 

MP 

mutation prob. 

-37.265 

PS 

{lopvilation size 

17.58 

PS*PS 

population size (quad) 

7.96 

MP*MP 

mutation jirob. ((jiiad) 

-6.83 

PS*MP 

population size * mutation prob. 

-2.86 


Table 4.11: Effects table for average case 3 
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Dc.sign parameter 

Effect 

PS 

PS*MP 

MP 

PS*PS 

population size 

population size * mutation prob. 
mutation prob. 
population size (quad.) 

-3.42 

3.24 

-2.05 

-1.45 


Table 4.12: Effects table for standard deviation case 3 


Response surface study 

Based on 27 tiials data in Table 4.10, we fit the best surface through the average. The 
re.sulting equation is : 


on-line performance = 52.534 -4- 0.5977P5 - C37.577MF - 0.00315P5'^ 

+2982.668MF^ - 0.6992P5 x MP (4.13) 


The best fit surface through standard deviation is also found using the data in Table 
4.10. The resulting equation is: 

log(c7) = 2.217 - 0.048P5 - 12.75MP -t- 0.0002P5^ -b 0.1429PS' x MP (4.14) 

Analysis by direct observation 

We see that MP and PS are VAPs as its effect are significant on standard deviation of on-line 
performance. The minimization of the equation 4.14 for standard deviation gives 

a = 0.52 , at P5 = 112 and MP = 0.001 

These value of PS and MP can be put into equation 4.13 of on-line performance then, 
on-line performance = 79.24 

We see that on-line performance doesn’t depend on CP. So, CP is not the key input 
variable. We can keep CP at any point, it will not affect either average or standard deviation. 
Tlie r('sults from this analysis can be concluded jus: 

• Set PS = 112 and MP = 0.001 in order to minimize variation. 
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Analysis by Variation transmission 


The ('(luiitiou for st.aiulard 


deviation of on-line performance can be found by equation 4.13. 


cr^ - (0.5977 - 2 X 0.00315ips - O.G992tMpfaps + (-637.577 

+2 X 2982.668^Arp - 0.6992ips)V^^p -I- 0.5 x (2 x .00315)Vp5 

+0.5 X (2 X 2982.6C8)Vl,p + {0. 6992)^ alscr^p (4.15) 


Putting aps — 5,crcp — 0.05, (jAfp — 0.005, in the equation 4.15 and minimizing it 

gives 

a = 0.51 at PS = 82, MP = 0.1. 

Putting those values of PS and MP in the equation 4.14 of standard deviation of on- 
line performance, we get a = 0.69. The standard deviation obtained from direct observation 
method was 0.52. 

By the equation 4.15, transmitted variation can be calculated by the method explained 
in chai)tcr 2. The transmitted variation are : 

aips = 0.1, atMP = 0.5, 

We see that transmitted variation by MP is more so there is need to reduce the variation 
given to MP as input. Hence MP is a VIP also. 

So with mean at PS = 82, CP = random value between 0.5 to 1, and MP = 0.1 and 
standard deviation aps = 5 and cmp = -005, the final capability study is done to verify the 
expected improvements. 

4. 1.4,1 Filial Capability Study for Case 3 

The data are given in Table 4.13. The control chart for average and standard deviation is 
also given in Figure 4.7 and Figure 4.8 . The process seems to be capable, as all points are 
within the UCL and LCL of average and standard deviation. Also the average and standard 
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of all trials are reduced a.rd are cousistent. We ucw try the next case. 



Table 4.13: Final capability study for case 3 function Fl 




Trial Number 


Figure 4.7: Control chart for average of on- Figure 4 8- Control oUart r . a j j • 

etc a^Xn “‘7 performance IZ “ariato 

reduction case 3 on function Fl 


4.1.5 CASE 4 (Tournament Selection and Uniform Crossover) 

27 trials of cxiicriincnt arc done with tonrnaincnt selection as selection operator and uniform 
crossover as crossover operator. The data are shown in Table 4.14. 

The significant effects of parameters, their interaction and quadratic effects arc shown 
in Table 4.15 for average, and in Table 4.16 for standard deviation. 
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'lYial 

no. 

■i 

tlidnie 

variable 




Run number 

on-line performance 



OI* 

■KUM 

1 

2 

3 


5“ 

^ 

^ 

8 



^ 

average 

»td. dev. 

1 

1 

0.50 

U.iOU 

34.7 

33.9 

33.9 

33.0 

33.8 

33.3 

34.2 

34.4 

33.9 

34.2 

33.93 

0,49 

2 

112 

1.00 

0.001 

77.8 

77.7 

77.5 

77.7 

77.7 

77.8 

77.7 

77.7 

77.7 

77.7 

77.70 

0.08 

3 

1 1 

0.50 

0.100 

22.1 

22.1 

23.2 

23.2 

21.9 

23.7 

23.8 

26.5 

21.7 

23.6 

23.17 

1.41 

4 

112 

0.60 

0.001 

77.6 

77.7 

77.5 

77.8 

77.7 

77.7 

77.7 

77.5 

■a 

74.7 

77.36 

0.93 

6 

11 

1.00 

0.001 

68.1 

50.2 

52.2 

56.1 

59.1 

58,2 

55.9 

52.2 

iii 

59.2 

58.00 

6.28 

6 

11 

1.00 

0.100 

25.4 

26.1 

26.5 

27.4 

31.2 

22.4 

27,1 

27.9 

W-iWM 

27.7 

26.78 

2.19 

7 

112 

1.00 

0.100 

31.8 

31.1 

31.2 

30.3 

31.5 

32.1 

31.3 

31.2 

31.3 

31.9 

31.37 

0.51 

8 

11 

0.50 

0.001 

69.1 

49.2 

59.0 

68.6 

58.5 

50.8 

57.7 

62.6 

48.7 

54.9 

56.89 

7.31 

9 

61 

0.75 

0.051 

53.2 

52.8 

50.6 

52.4 

51.1 

51.6 

50.0 

51.9 

51.7 

52.3 

51.76 

0.99 

10 

61 

0.50 

0.001 

76.4 

76.6 

77.1 

76.7 

71.8 

76.5 

76.4 

76.4 

76.8 

76.7 

76.13 

1.53 

11 

61 

0.50 

0.100 

33.4 

34.4 

31.1 

33.2 

32.1 

32.4 


33.1 

32.6 

34.1 

33.02 

1.02 

12 

61 

1.00 

0.001 

76.6 

76.8 

77.0 

76.8 

76.0 

76.9 

76.5 

76.5 

76.8 

76.0 

76.57 

0.34 

13 

61 

1.00 

0.100 

31.3 

29.4 

33.6 

29.0 

30.3 

30.5 

30.1 

31.0 

31.8 

31.6 

30.83 

1.29 

14 

1 1 

0.75 

0.001 

08.7 

54.9 

58.4 

59.9 

50.3 

53.8 

54.5 

61.4 

56.1 

48.7 

56.50 

5.84 

15 

U 

0,75 

0.100 

25.7 

22,9 

23.0 

26.5 

v3Q.7 

25.8 

27.6 

24.5 

28.6 

27.3 

26.20 

2,44 

10 

112 

0.75 

0.001 

77.7 

77.7 

77.7 

77.8 

77.7 

77.8 

77.7 

77.8 

77.7 

77.7 

77.73 

0.03 

17 

112 

0.75 

0.100 

31.5 

31.9 

31.1 

32.6 

33.9 

31.7 

32.3 

32.9 

32.6 

31.6 

32.21 

0.83 

18 

11 

0.50 

0.051 

28.4 

36.1 

35.1 

37.4 

38.6 

30.3 

36.9 

28.6 

29.5 

36.1 

34.31 

3.89 

19 

11 

1.00 

0.051 

25.8 

28.0 

35,9 

35.1 

38.0 

35.2 

29.6 

28.4 

36.7 

46.7 

33.95 

6,18 

20 

112 

0.60 

0.051 

55.5 

52.8 

53.8 

53.3 

54.6 

52.7 

■tKW 

53.0 

63.9 

53.8 

53.65 

0.80 

21 

112 

1.00 

0.051 

54.4 

55.4 

54.4 

52.0 

53.4 

53.5 

65.0 

54.1 

54.6 

64,3 

54.11 

0.95 

22 

61 

0.75 

0.00 1 

70.6 

76,8 

76.8 

76.8 

70.2 

70.9 I 

76.5 

um 

76.8 

76.6 

76.05 

2.07 

23 

01 

0.75 

0.100 

31.7 

30.7 

32.7 

30.0 

31.9 

31.5 ; 

31.8 

31.6 

31.8 

31.9 : 

31.57 : 

0.74 

24 

11 

0.75 

0.051 

28.8 

24.6 

26.5 

26.8 

38.7 

34.8 

35.9 

36.0 

28.7 

36.1 

31.69 

5.10 

26 

112 

0.75 

0.051 

55.0 

55.1 

54.1 

53.3 

55.1 

53.5 

53.4 

55.1 

54.6 

64.3 

54.34 

0.72 

26 

61 

0.50 

0.051 

53.6 

51.8 

53.4 

52.8 

49.6 

53.1 

62.9 

60.5 

51.5 

62.3 1 

52.15 

1.31 

27 

61 

1.00 

0.051 

49.8 

48.7 

54.4 

50.5 

50.7 

53.1 

50.8 

51.7 

60.3 

52.3 ! 

51.23 

1.67 


Table 4.14; Data obtained from 27 trials on DeJong FI function case 4 


Design paraiiieter 

Effect 

MP 

nmtatioii prob. 

-40.428 

PS 

population size 

16.087 

PS*MP 

population size * mutation prob. 

-6.67 

PS*PS 

population size (quad) 

6.589 

MP*MP 

mutation prob. (quad) 

-3.76 

PS*CP 

population size * crossover prob. 

-1.01 


Table 4.15: Effects table for average case 4 


Design parameter 

Effect 

PS 

population size 

-3.91 

PR*MP 

j)OI)ulation size * nnit.T.tion prob. 

2.3G2 

MP 

mutation prob. 

-1.49 

PS*PS 

population size (quad.) 

-1.33 


Table 4.16: Effects table for standard deviation case 4 
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Response surface study 


Based on 27 trials data in Table 4.14, 
resulting equation is : 


wo fit the best surface througii the average. 


Tlie 


on-line performance = 50.29 -b 0.5557P5 - 338.81MP - 0 . 0026 ^ 5 ^ + j p 2 

-1-32PSxMP-QM19PSxCP (4.16) 

The best fit surface equation through standard deviation is: 

log(or) = 2.3788 - 0.0389PS' - 13.047MP -b 0.00009Pp2 + o.ieGP^ x MP (4.17) 
Analysis by direct observation 

Analysis by direct observation is done in same way as in case 3, i.e. the equation 4.17 is 
niinirnized. These values of PS and MP obtained for minimized standard deviation is put 
into equation 4.16. Then resulting equation in terms of only CP, is minimized for minimum 
on-line performance. The result from this analysis can be concluded as: 

• Set PS = 112 and MP = 0.001 in order to minimize variation. 

• Set CP = 1 to reduce the on-line performance. 


Analysis by Variation transmission 

The equation for standard deviation of on-line performance is obtained from equation 4.16. 
Putting aps = ^,crcr — 0.05, ctmp = 0.005, in the resulting equation and minimizing it 
gives: 

cr = 2.04 at PS = 60, CP = 1.0, and MP = 0.1. 

Putting these values of PS and MP in the equation 4.17, we get cr = 1.064. The value 
of standard deviation obtained from direct observation analysis is 0.43. This means that 
result from variation transmission method is not good. This shows that direct observation 
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result is better for ininiinuni standard deviation. Since we are assuming that the result of 
variation transmission analysis is not good, hence there is no need to calculate transmitted 
variation of PS ,CP and MP. 

So with mean at PS = 112, CP = 1, and MP = 0.001 and standard deviation <Jp 5 = 5, 
(Jcp — O.Oo and ctmp = 0.005 final capability study is done to verify the expected iinprove- 
nuMit. 


4. 1.5.1 Final Capability Study for Case 4 

The data are given in Table 4.17. The control chart for average and standard deviation is 
also given in Figure 4.9 and Figure 4.10. The process seems to be capable except one point is 
falling slightly outside UCL of standard deviation. We see that average on-line performance 
is higher compared to case 1, 2 and 3. 


HusBBSEIi 

Hun number j 

average 

mrntm 

. 1. 

2 

3 

4 



7^ 


r- 

10 

1 

77.20 

77.14 

77.21 

77.26 

77.24 

77.25 

77.17 

77.22 

77.34 

77.26 

77.24 

KMX 

2 

77.92 

76.33 

77,87 

77.92 

77.92 

77.93 

77.92 

77.92 

77.94 

77.93 

77.76 

0.50 

3 

77.08 

77.02 

77.02 

77.33 

77.28 

■ia 

77.91 

77.33 

77.59 

77.67 

77.57 

0.20 

4 

77.00 

77.58 

77.09 

77.62 

77.22 


77.61 

77.65 

77.75 

77.68 

77.60 

0.15 

fi 

77.03 

77.59 

77.20 

77.70 

77.27 


77.91 

77.33 

77.18 

77.67 

77.62 

0.24 

0 

70.84 

70.86 

76.01 

76.85 

76.40 

76.82 

77.01 

76.76 

miMm 

76.79 

76,77 

0.16 

7 

76.79 

77.28 

77.18 

77.22 

76.98 

77.19 

77.14 

77.22 

77.26 

77.21 

77.15 

0.15 

8 

77.65 

77.02 

77.62 

77.70 

77.28 

77.64 

77.91 

77.33 

77.59 

77.67 

77.60 

0.18 

9 

75.84 

75.88 

75.50 

75.99 

75.57 

75.92 

75.84 

75.98 

75.78 

75.95 

75.83 

0.17 

10 

77.92 

77.92 

77.90 

77.92 

77.71 

77.92 

mkMm 

77.92 

77.92 


77.90 

0.06 

11 


75,74 

75.69 

75.75 

75.93 

75.83 

75.86 

75.88 

75.54 

75.86 

75.80 

0.12 

12 

76.16 

76.24 

76.31 

76.27 

76.37 

76.30 

76.36 

76.21 

76.20 

76.37 

76.28 

0.08 

13 

77.62 

77.54 

77.48 

77.54 

77.49 

mmm 

77.53 

77.48 

77.54 

77.56 

77.54 

0.05 

14 

77.92 

70.33 

77.91 

77.92 

77.92 

77.93 

77.92 

77.92 

77.94 

77.93 

77.76 


15 

77.03 

77.62 

77.26 

77.70 

77.07 

77.66 

77.91 

77.33 

77.18 

77.67 

77.50 

0.27 1 

10 

70.82 

76.81 

76.57 

76.80 

76.39 

76.77 

76.97 

76.71 

76.68 

76.74 

76.73 


17 

70.71 

77,25 

77.17 

77.19 

76.95 

77.18 

77.12 

77.15 

77.25 

77.18 

77.12 

0.16 

18 

77.68 

77.02 

77.60 

77.70 

77.28 

77.64 

77.91 

77.33 

77.18 

77.67 

77.56 

0.22 

19 

75.76 

75.80 

75.76 

75.94 

75.50 

75.82 

75.73 

76.91 

75.68 

75.86 

75.78 

0.13 

20 

77.92 

77,92 

77.90 

77.92 

77.92 

77.92 

77.91 

77.88 

77.92 

77.93 

77.91 i 

■mif 


Table 4.17: Final capability study for case 4 function FI 


4.1.6 Discussion of Result of Function FI 

In the Table 4.18 the average and the standard deviation of the final capability study is 
shown for each case. The value of PS, CP and MP are given at which we have to set it with 
standard deviation of 5, 0.05 and 0.005 respectively for variation reduction. 
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Figure 4.9: Control chart for average of on- Figure 4.10; Control chart for standard devi- 
line perlorinance after variation reduction ation of on-line performance after variation 
case 4 on function FI reduction case 4 on function FI 



PS 

CP 

MP 

average 

std. deviation 

CASE 1 

Proportionate selection 
single point crossover 

11 

0.5 

0.1 

17.96 

2.086 

CASE 2 

l)roportionate selection 
uniform crossover 

105 

0.54 

0.0818 

24.83 

0.414 

CASE 3 

tournament selection 
single point crossover 

82 

(0.5 to 1) 

0.1 

37.55 

0.7211 

CASE 4 

tournament selection 
uniform crossover 

112 

1.0 

0.001 

77.14 

0.218 


Table 4.18: Result of function FI 
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We see that if oui aim is to minimize only on-line performance then we will choose 
pr()i)oitionat(^ selet.tion as selection oi)erator and single point crossover as crossover operator 
(case 1). But if our aim is to minimize only variation in on-line performance then tournament 


selection with uniform crossover (case 4) will be chosen. But in this case on-line performance 
is very high. So we will choose the case in which variation in on-line performance is less 
and also on-line peifoimance is less. So for function Fl proportionate selection with uniform 
crossover (case 2) will give better result. 


In every case PS, CP and MP which arc to, be targeted, vary. Since we found that 
proportionate selection with the uniform crossover is better, we will set PS at 105, CP at 
0.51 and MP at 0.081 as mea.n with standard deviation of 5 for PS, 0,05 for CP and 0.005 
for MP. We now see the variation reduction process applied on De Jong Function 2. 


4.2 Test Function 2: De Jong Function F2 


4.2.1 Initial Capability Study 


To check the need for variation reduction a capability study is done. The on-line performance 
measured for last 50 % of the generation are shown in Table 4.19. The average and standard 
deviation of each trial are also shown in Table 4.19. 


These data arc used to construct the average and standard deviation control charts shown in 
Figure 4.11 and Figure 4.12. We see that lot of points are going outside the control limits. So 
there is need to do variation reduction on this function. For reducing the variation four cases 
as described in previous chapter is selected and for each case variation reduction technique 
is applied. 
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TrmI 

number 


“T 


Run number 


T 


'T 


"ir 


average 


“itlT 

dev. 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 
1 1 
12 

13 

14 

15 

16 

17 

18 

19 

20 


1254.7 

220.7 
11.7 

168.7 
25.6 

1493.9 

273.4 
96.2 

3273.0 

3201.4 

1806.0 

3083.2 
3500.6 

224.4 

1381.5 
2430.0 

163.9 

77.5 

1563.2 

1665.9 


1275.9 

186.4 
50.0 

163.9 
89.9 

1491.3 

251.9 

100.0 

3132.2 

3312.4 

1812.5 

3011.6 
3514.4 

233.4 

1358.3 

2406.6 

156.4 
170.3 

1576.1 

1542.0 


1366.7 

203.6 
3.7 

151.9 

40.3 

1421.7 
273.4 
112.3 

3249.2 
3284.6 

1862.4 

2980.5 

3481.6 

230.6 

1328.2 

2355.2 

202.2 

88.9 

1590.4 

1483.0 


1331.4 
160.0 

10.8 

164.4 

90.1 

1489.6 

264.0 

112.1 

3263.6 
3311.8 
1831.3 

3076.2 
3496.1 

237.3 

1392.3 

2417.5 
204.2 

30.0 

1543.6 
1632.5 


1251.0 

187.4 
2.4 

177.4 
144.8 

1474.7 

296.0 

108.1 

3254.2 
3310.6 

1759.5 
3066.9 

3538.2 
238.1 

1393.6 

2344.1 
152.7 

81.2 

1545.2 
1626.1 


1253.5 

176.3 
1.7 

180.4 
69.6 

1463.1 
258.6 

91.0 

3279.2 

3300. 1 

1777.5 

3037.6 
3466.0 

279.3 

1350.5 

2485.8 

173.5 
62.9 

1688.7 

1455.2 


1175.9 

188.5 

7.1 

169.8 

104.7 

1410.2 
265.3 

116.9 

3208.6 

3306.8 

1750.3 

3056.6 

3514.6 

236.8 

1372.2 
2376.5 

201.0 

36.4 

1536.3 

1685.9 


1240.9 

187.8 
5.2 

172.4 

90.2 

1465.7 

258.8 

116.7 
3306.6 

3322.0 

1838.0 

3110.9 
3507.2 

228.6 

1280.4 

2499.5 

171.7 
78.7 

1525.1 
1493.0 


1305.7 
189.4 

4.8 

158.6 

96.2 

1481.3 

285.3 

105.0 
3247.1 

3.3 U . 5 

1785.6 

2980.6 

3484.5 

257.3 

1217.8 

2389.0 

108.1 
69.6 

1554.0 

1600.6 


1124.1 

195.6 
2.7 

176.7 

135.4 
1518,8 

291.8 

100,1 

3257.3 

3274.3 
1856.5 

2987.2 

3467.4 

248.5 
1328.1 
2404.7 

167.8 
40.1 

1525.4 
1688.0 


1258.0 
189.6 

10.0 

168.4 
88.7 

1471.0 

271.8 

105.9 

3253.2 
3302.6 

1808.0 
3040.0 

3497.9 

241.4 

1340.3 

2411.9 
176.2 

74.2 

1554.8 

n n 


70.6 

15.9 

14.5 
8.9 

37.0 

33.2 

15.0 
8.9 

45.9 

14.8 

30.2 

46.8 

23.1 

10.4 

55.4 

51.1 

19.2 

39.0 

24.1 

80.6 


Table 4.19: Initial capability study for function F2 




Trial Number 


Figure 4.11: Control chart for average on- Figure 4.12: Control chart for standard devi- 
line performance before variation reduction ation of on-line performance before variation 
on function F2 reduction on function F2 
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4.2.2 CASE 1 (Proportionate Selection and Single Point Crossover) 


27 trials of exj)eiiniciit is done with proportionate selection as selection operator and single 
point crossover as crossover operator. For each trial the combination of variable is taken 
Iroin fable 4.20 and ten runs are perlornred by changing the initial populatiorr. The average 
atrd the standard deviation are calculated for each trial. The data are shown in Table 4.20. 


TrinP 

no . 

cni 

(Jlclni e 
vnriabl < 

npnt 

'S 



Rttn nuinhor 

on-line performance 


PS 

CP 

M P 

1 

2 

3 

4 



7~ 



9 

10 

1KSS3£I^33II 


1 

iuo 

0.50 

0.100 

284.9 

298.1 

300.0 

302.6 


297.9 

288.5 

308.7 

319.0 

315.4 

Milil'iMiill M 


2 

100 

1.00 

0.001 

58.3 

79.5 

49.1 

42.3 

52.8 

38.4 

62.4 

14.9 

11.6 

64.2 

47.34 

21.4.3 

5 

10 

0.50 

0,100 

205.7 

135.7 

194,8 

141.3 

146.3 

156.6 

143.2 

137.9 

160.3 

202.1 

161.30 

28.0.3 

4 

100 

' 0.50 

0.001 

20.0 

2.H.7 

48.7 

10.6 

121.9 

40.2 

20.8 

77.8 

25.5 

.30.5 

43.04 

.33.0.3 

5 

10 

1 .00 

O.OOl 

l.O 

11. 6 

2.6 

0.6 

6.9 

0 1 

0.1 

1.8 

2.4 

5.2 

3.23 

.3.69 

i ] 

10 

1.00 

0.100 

197.2 

239.2 

168.3 

141.9 

147.9 

140.1 

145.6 

301.1 

141-0 

238.7 

185.10 

56.66 

7 

100 

1 .00 

0.100 

274.3 

306.1 

303.6 

327.7 

328.6 

297.7 

292.9 

317.8 

290.1 

286.6 

v302.55 

17.88 

8 

10 

0.50 

0.001 

94.8 

2.4 

6.4 

0.6 

223.1 

76.9 

1.4 


3.6 

1.0 

41.08 

72.92 

9 

56 

0.75 

0.051 

201.8 

196.3 

191.5 

19G.4 

199.5 

208.6 

228.5 

198.3 

192.1 

197.0 

201.00 

10.84 

10 

65 

0.50 

0.001 

14.7 

22.0 

34.1 

24.8 

81.7 

34.0 

37.9 

21.3 

32.8 

69.2 

37.30 

21.56 

11 

55 

0.50 

0.100 

209.G 

254.9 

263.2 

268.3 

252.9 

263.4 

252.8 

258.0 

281.6 

271.0 

263.55 

9.29 

12 

55 

1.00 

0.001 

1.3 

0.6 

3.8 

3.1 

80.6 

2.3 

2.9 

1.5 

65.2 

15.6 

17.69 

29.64 

13 

55 

l.OO 

0.100 

280.3 

249.2 

283.3 

298.8 

248.0 

2G4.6 

205.2 

279.4 

282.3 

288.7 

273.98 

16.70 

14 i 

10 

0.75 

0.001 

7.3 

83.0 

0.5 

2.6 

1.7 

5.9 

16.9 

0.5 

3.9 

24.2 

14.66 

25.23 

15 

10 

0.75 

0,100 

140.4 

157.6 

174.3 

148.5 

134.7 

141.7 

236.2 

168.6 

161.8 

140.8 

160.45 

29.71 

10 

100 

0.75 

0.001 

73.9 

29.9 

61.8 

33.0 

88.2 

48.3 

47.1 

48.7 

53.6 

54.1 

53.87 

17.49 

17 1 

100 

0.75 

0.100 

298.4 

293.5 

306.9 

284.1 

278.9 

316.6 

■cltWCM 

301.4 

281.6 

301.2 

296.31 

11.89 

18 

10 i 

0.50 

0.051 

94.6 

197.4 

83.3 

97.4 

112.4 

71.3 

168.4 

218.9 

79.8 

124.8 

124.83 

52.11 

19 

10 

1.00 

0.051 

140.7 

117.1 

105.0 

122.3 

75.9 

123.2 

96.0 

133.8 

113.5 

116.6 

114.42 

18.65 

20 

100 

0.50 

0.051 

240.3 

209.5 

228.9 

208.5 

240.6 

224.9 

229.3 

194,1 

253.9 

231.7 

226.19 

17.80 

21 

100 

1.00 

0.051 

269.4 

250.9 

239.3 

233.1 

243.0 

221.3 

225.8 

240.2 

230.5 

245.5 

239.91 

13.79 

22 

55 

0.75 

0.001 

30.6 

54.9 

2.4 

2.5 

49.6 

10.0 

1.5 

11.0 

14.3 

4.4 

18.13 

19.97 

23 

55 

0.75 

0.100 

268.5 

264.3 

284.7 

263.1 

280.4 

260.0 

287.8 

251.6 

255.8 

277.7 

269.40 

12.55 

24 

10 

0.75 

0.051 

78.8 

75.3 

141.3 

176.6 

133.5 

85.2 

101.5 

97.8 

82.3 

198.3 

117.05 

43.51 

25 

100 

0.75 

0.051 

211.8 

226.3 


248.4 

235.6 

214.6 

236.0 


249.2 

243.7 

233.58 

17.02 

26 

55 

0.50 

0.051 

177.9 

160.0 

216.0 

196.6 

183.2 




187.6 

177.9 1 

195.21 

22.21 

27 

55 

1.00 

0.051 

199.3 

192.3 


211.2 

208.6 

EIBI 

212.0 

w£yi 

220.7 

210.6 

206.33 

7.99 


Table 4.20: Data obtained frorrr 27 trials on DeJorrg F2 furretion case 1 


The signilica.nt effects of parameters, their irrtcractiorr atrd the quadratic effects are 
shown in Table 4.21 for average and in Table 4.22 for standard deviation. 


Design parameter 

Effect 

MP 

niutation prob. 

215.11 

PS 

population size 

91.35 

PS*MP 

population size * mutation prob. 

51.07 

MP*MP 

mutation prob. (quad) 

45.92 

PS*PS 

population size (quad) 

16.596 

CP*MP 

cro.ssovcr prob.* mutation prob. 

15.03 

ps*cr 

population size * crossover prob. 

7.25 


Table 4.21: Effects table for average case 1 
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Design parameter 

Effect 

CP*MP 

crossover prob.* mutation prob. 

19.1 

PS 

population size 

18.68 

PS*CP 

population size*crossover prob. 

10.72 

PS*PS 

population size (quad.) 

-10.63 

CP 

crossover prob. 

-9.167 

PS*MP 

population size* mutation prob. 

-7.124 

MP 

mutation prob. 

-5.58 

CP*CP 

crossover prob. (quad.) 

-4.38 


Inblc 4.22: li/ffccts table for standard deviation case 


Response surface study 


Based on 27 trials data in I able 4.20, we fit the best surface through the average. The 


resulting eejuation is : 


on-line performance = 4 + 1A7PS + 1784.56MP - O.OOSIPP^ - 7332.14MP^ 

-0.2589PP X CP + 549.91CP x MP 

+12.58P5 X MP (4.18) 


The best fit surface through standard deviation is: 


log((T) = 16.37 - 0.044P5 - 29.32CP - 58.4MP -t- 0.000197P5'2 + IS.OSICP^ 

-t-().029PP X CP + 82.4CP x MP - 0.1P5 x MP (4.19) 


Analysis by direct observation 


The result from this analysis can be concluded as: 


• Set PS = 90, CP = 0.64 and MP = 0.0991 in order to minimize variation. 


Analysis by Variation transmission 


Minimization of the equation for standard deviation obtained from equation 4.18 gives 

central uBRAm 

f I r., (CANPUii 
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a - 10.43 at PS = 82, CP = 0.5, MP = 0.1 

The standard deviation from equation 4.19 will be 5.98, at PS = 82, CP =0.5 and 
MP=0.1, while standard deviation from direct observation method was 4.773. This means 
that result from analysis of variation transmission method is not bad. The transmitted 
variation of each parameter w^as found to be: 

(^trs = 6-318, (Tfcp = l.C94,crjA/p = 8.13 

Ihis shows that PS and MP arc VIPs. So there is need to reduce the variation of 
PS and MP. So with mean at PS = 82, CP = 0.5 and MP = 0.1 and standard deviation 
O' PS — 6 , Ocp = 0.05 , and omp = 0.005 final capability study is done to verify the expected 
improvement. 

4.2. 2.1 Final Capability Study for Case 1 


The data are shown in Table 4.23. The control chart for average and standard deviation is 
also given in Figure 4.13 and Figure 4.14. The process is capable as all points are within 
control limits except one point which is slightly out of UCL limit in control chart for standard 
deviation. We now try the next case. 


Trial 

number 

Run number | 

average 








umim 




1 

299.1 

295.2 

301.1 

290.9 

294.2 

309.9 

286.7 

287.9 

324.5 

305,9 

299.6 

11.5 

2 

305.9 

304.9 

278.9 

mklmm 

289.5 

279.1 

279.0 

277.3 

291.0 

297.4 

287.9 

11.5 

3 

296.2 

296.9 

315.1 

272.3 

256.5 

272.4 

283.4 

266.1 

263.1 

299.8 

282.2 

19.1 

4 

297.5 

285,2 

281.4 

280.8 

301,3 

295.8 

280,6 

267.7 

284.6 

296.9 

287.2 

10.4 
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'lable 4.23: Final capability study for cfise 1 function F2 
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Trial Number 


Figure 4.13: Control chart for average of on- Figure 4.14: Control chart for standard devi- 
line perlonna.nco alter variation reduction ation of on-line i)crforinance after variation 
lor case 1 on lunction F2 reduction for case 1 on function F2 


4.2.3 CASE 2 (Proportionate Selection and Uniform Crossover) 

27 trial experiment data are shown in Table 4.24. 
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Table 4.24: Data obtained from 27 trials on DeJong F2 function case 2 

The significant effects of parameters, their interaction and the quadratic effects are 
shown in Table 4.25 for average and in Table 4.26 for standard deviation. 
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Design parameter 

Effect 

MP 

mutation prob. 

207.4 

PS 

population size 

96.6 

PS*MP 

population size * mutation prob. 

78.32 

MP*MP 

mutation prob. ((iiiacl) 

46.66 

PS+CP 

I)opulation size * crossover prob. 

16.39 

PS*PS 

I)oi)nIation size (quad) 

12.19 

CP*MP 

cro.Hsovcr prob. * mutation j)rob. 

10.29 


Table 4.25: Effects table for average case 2 


Table 4.26: Effects table for standard deviation case 2 


Design parameter 

Effect 

MP 

mutation prob. 

-64.73 

PS 

population size 

-61.44 

PS*MP 

population size * mutation prob. 

59.7 

PS*PS 

population size (quad) 

-36.77 

MP*MP 

mutation prob. (quad) 

-26.127 

PS*CP 

population size * crossover prob. 

14.17 

CP*MP 

crossover prob. * mutation prob. 

10.38 

CP 

crossover prob. 

-9.88 







4.2 Test Function 2: De Jong Function F2 


55 


Response surface study 

Based on 27 trials data in Table 4.24, the best surface through the average is fitted. The 
resulting equation is : 

on-line performance = 51.015 + 0.8585P5 + 1639.9MP - O.OOfiPP^ - 8034.29Mp2 

-0.082P5 X CP -I- 281.23CP x MP 

-bl8.725P5 X MP (4.20) 


The best lit surface through the standard deviation is: 

log(fT) = ().71 - OXmPS - 1.G19CP - 30.08A/P + 0.00039P5^ + 22.237 MP'^ 

-hO.OlPC X CP -I- 15.546CP x MP + 0.065P5 x MP (4.21) 

Analysis by direct observation 

• Set PS = 55 , CP = 0.5 and MP = 0.1 in order to minimize variation. 

Analysis by Variation transmission 

Miniinization of the cciuation for standard deviation obtained from equation 4.20 gives 

a = 9.70 at PS = 10, CP = 0.5, MP = 0.02 

The standard deviation from equation 4.21 will be 133.38 at PS = 10, CP = 0.5 and 
MP = 0.02, while value of standard deviation obtained from direct observation method is 
7.917. This means the result from Variation method method is not good. 

So with mean at PS = 55, CP = 0.5 and MP = 0.1 and standard deviation <jps = 
5, acp = 0.05, and cta/p = 0.005 final capability study is done to verify the expected 


iniproveiiient. 
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4.2.3. 1 Final Capability Study for Case 2 

The data are given in Table 4.27. The control chart for average and standard deviation is 
also given in Figure 4.15 and Figure 4.16. The process is capable as all points fall within the 
control limits. Also the standard deviation is less compared to initial capability study. We 
now try the next case. 


'IVinl 

iwuubcr 





Him uutuhr>r “ 

Rvernga 

M t (1 . 

1 

2 

8 

4 

5 

6 

7 

5 

g 

10 

1 

2 

3 

4 

5 

a 

1 

u 

10 

1 1 

12 

13 

M 

15 

10 

17 

la 

19 

20 

300.0 
276.9 

291.8 

304.2 

302.2 

300.9 

300.5 

251.5 

314.2 

333.8 

292.4 

314.5 

298.4 

291.0 

291.5 

300.9 

300.7 
254.4 

330.9 

303.7 

27 4.8 

293.3 
314.9 
281.1 

287.5 
3 U 1.4 

295.2 

208.7 

311.5 

287.4 

303.0 

282.1 

287.0 

290.7 

288.2 

299.0 
307.2 

250.4 

320.5 

297.1 

279.5 

310.5 

314.5 

304.6 

327.0 

319.0 

295.6 

271.4 

299.6 

208.5 
280.2 

287.1 

285.5 

300.4 
278.0 

307.4 

281.2 

273.4 

313.6 

300.4 

283.1 

282.7 
278.9 

319.8 

293.4 

341.8 

298.0 

205.5 

317.2 

278.3 

322.9 

300.9 

297.1 

287.2 

272.2 

296.6 

300.2 

259.2 

302.8 

292.9 

291.2 

285.2 

346.2 

319.6 

302.5 

310.8 
■ 340.4 

304.5 

317.0 

302.4 

311.1 

302.8 

297.7 

285.4 

270.1 
293.0 

308.9 

263.3 

311.8 

294.4 

312.1 

272.9 

301.1 

296.9 

274.8 

307.2 

340.0 

200.9 

325.7 

324.3 

314.4 

293.0 

295.1 
281.6 

274.7 

297.8 
311,3 

279.1 

309.5 

310.2 

294.8 

280.6 

275.7 
302.3 

284.6 

318.0 

300.0 

282.8 

309.2 

290.0 
319.9 

300.0 

320.7 

298.3 

269.2 

308.2 

273.8 

280.9 

332.7 

299.8 

297.8 

305.3 

288.8 
312.9 

283.1 

308.1 

273.7 

295.2 
315.0 

285.2 

357.8 

286.4 

288.7 

318.7 

266.3 

276.4 

292.7 

269.7 

294.9 

263.9 

299.8 

303.0 

284.2 

280.8 

275.4 

292.7 

287.7 

250.5 

318.8 

281.1 

333.2 

288.7 

311.5 

293.2 

289.1 

299.9 

208.9 

276.1 

313.5 

292.7 i 

282.1 

269.7 

297.9 
282.2 
265.2 

300.4 

277.4 

274.1 

307.2 

200.2 

320.0 

303.1 

312.0 
200.3 

298.7 

305.0 

286.9 

270.0 

321.2 

270.5 

291.5 
288.0 

299.4 

300.4 

289.6 

311.6 

301.6 

272.5 

313.6 

294.7 
315,5 
295,9 

300.0 
290. T 
280.4 

299.0 

296.8 
268.2 

315.1 

293.2 

11.5 

14.3 
21.2 

15.0 

17.7 

13.6 

23.9 

17.6 
7.2 

20.6 

21.5 

10.0 
13.0 

14.9 

10.8 

9.5 
12.8 
9.5 
11.7 

15.4 


Table 4.27: Final capability study for case 2 function F2 




Figure 4.15: Control chart for average of on- Figure 4.16: Control chart for standard devi- 
line performance after variation reduction ation of on-line performance after variation 
case 2 on function F2 reduction case 2 on function F2 


4.2.4 CASE 3 (Tournament Selection and Single Point Crossover) 

27 trials experiment data are shown in table 4.28. 
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1242.12 
1225.78 
2562.89 
2435.11 
2474.01 

35.12 

49 83 
48.56 
74.81 

718.02 
28.59 
39.35 

859.41 

53.91 
286.21 

48.34 

251.57 

52.42 

582.02 

37.88 
60.02 

47.91 
80.15 

93.89 
35.55 
53.26 

221.68 

69.85 

72.44 

84.32 

51.87 

66.72 


Table 4.28: Data obtained from 27 trials on DeJong F2 function case 3 

The significant effects of parameters, their interaction and the quadratic effects are 
shown in Table 4.29 for average and in Table 4.30 for standard deviation. 


Design parameter 

Effect 

MP 

mutation prob. 

-2281.87 

PS 

population size 

1070.61 

PS*PS 

population size (quad) 

359.87 

MP*MP 

mutation prob. (quad) 

-153.184 

PS*MP 

population size * mutation prob. 

-120.83 

CP*MP 

crossover prob. * mutation prob. 

-89.17 

PS*CP 

population size * crossover prob. 

-71.006 

CP 

crossover prob. 

64.8 

CP*CP 

crossover prob. (quad.) 

31.16 


Table 4.29: Effects table for average case 3 
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Design parameter 

Effect 

PS*MP 

population size mutation prob. 

330.35 

MP 

mutation prob. 

-299.5 

PS 

population size 

-266.75 

MP*MP 

mutation prob. (quad) 

-129.29 

PS*PS 

population size (quad.) 

-44.22 

CP*MP 

crossover prob. * mutation prob. 

31.55 

PS*CP 

population size * crossover prob. 

24.09 

CP*CP 

crossover prob. (quad.) 

-22.98 

CP 

crossover prob. 

-18.48 


Tal)lc 4.30: 


Effects table for standard deviation case 3 


Response surface study 

Based on 27 (, rials data in fable 4.28, the best surface through the average is fitted. The 
resulting cciuation is : 

on-line pcrfonuancc = 1084.62 + 35.83P5 4- 2886.36C7P - 19469.5MP - 0.178P5^ 

- 1516.32CP2 -F 22060.7MP2 - 3.65P5 x CP 
-5342.28CP X MP - 30.91P5 x MP (4.22) 


The best, fit surface t.hrough standard deviation is: 


log((7) = 3.9C8 - 0.017P5 + 7.19CP - 36.78MP - 0.000071P52 - 4.097C'p2 

-I-102.2CMP''* 4- 0.0022P5 x CP - 5.29CP x MP -F 0.264P5 x MP(4.23) 


Analysis by direct observation 

The results froui this analysis can be concluded as: 

• Set PS = 10 , CP = 1 and MP = 0.1 in order to minimize variation. 


Analysis by Variation transmission 
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Minimization of the equation for standard deviation of on-line performance obtained by 
cciuation 4.22 gives: 

O' = 105.13 at PS = 82, CP = 0.57 and MP = 0.1. 

Putting the value PS — 82, CP = 0.57 and MP = 0.1 in the equation 4.23 we get 
a — 51.2. The value of standard deviation obtained by direct observation method is 29.34. 
This moans that tin; result from variation transmission method is not good. 

So with mean at PS = 10, CP = 1, and MP = 0.1 and standard deviation aps = 5, 
f^cp — 0.05 and (Ta//> = 0.005 final capability study is done to verify the expected improve- 
ment. 

4. 2.4.1 Final Capability Study for Case 3 


The data arc given in Table 4.31. The control chart for average and standard deviation is 
also given in Figure 4.17 and Figure 4.18. The process in this case is not capable since some 
points in control chart are marginally out of the control limits. This means that tournament 
selection with single point cross over for this function does not gives capable process. We 
now try for next case. 
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Ta,ble 4.31: Final capalulity study for case 3 function F2 














4.2 'rest" l‘^nic-Uoii 2: Do. Jong Function F2 


GO 




Trial Number 


Figure 4.17: Cont rol cliart for average of on- Figure 4.18: Control chart for standard devi- 
liiio, performance after variation reduction ation of on-line performance after variation 
Cixse 3 on function F2 reduction case 3 on function F2 


4.2.5 CASE 4 (Tournament Selection and Uniform Crossover) 


27 trials expcriinent data are shown in Table 4.32. 


I8BI 

enn 

cli(inlQ input 
varinhlcfl 

Run number 


HI 

mtim 

CP 

MP 

1 

2 

3 

4 


g- 


g- 


nr 



r 

100 

0.60 

O.lOO 

1304.2 

“TTJTT” 

1395.9 

1323.3 

1365.2 

1385.8 

1369.7 

1451J 

1420.6 

1307,7 

1381.02 

45 (50 

2 

100 

1.00 

0.001 

3844.4 

3860.2 

3846.4 

3844.1 

3853.9 

3844.0 

3855.0 

3849.8 

3841.3 

3844.0 

3847.93 

5.41 

3 

10 

0.60 

0.100 

002,5 

507.7 

502.7 

570.5 

620.3 

528.4 

000.6 

599.0 

611.7 

018.9 

594.29 

27 15 

4 

100 

0.50 

0,001 

3842.1 

3841 8 

3840.3 

3838.0 

3654.8 

3835.1 

3788.4 

3847.2 

3840.3 

3840 9 

3837. 80 

18 28 

5 

10 

1.00 

0.001 

3332.3 

2002.5 

3103.4 

3340.3 

1752.0 

3440.7 

1979,2 

2693.6 

3421.3 

2027.8 

2846.20 

009 70 

0 

10 

l.OO 

0.100 

538.7 

555.7 

620.7 

020.3 

505.0 

532.9 

497.2 

691.0 

580.8 

574.0 

572 30 

40 04 

T 

100 

l.OO 

0.100 

1372.0 

1323.0 

1328.4 

1340.5 

1417.0 

1401.2 

1387.0 

1372 2 

1438.7 

1401.7 

1384.90 

45 07 

A 

10 

0.50 

0.001 

3432,0 

3183.3 

3008.1 

1808.0 

2265.0 

2410.8 

3470,7 

1680.8 

3207.5 

3325.7 

2783.90 

714 17 

0 

5fl 

0.75 

0.051 

2020.7 

2300.0 

2418.0 

231)7.1 

2375.9 

2305.1 

2348 8 

2422.4 

2520.3 

2374.4 

2428.15 

85 HU 

10 

65 

0.50 

0.001 

3803.4 

3054.4 

3735.0 

3073.0 

3684. 0 

3730.4 

3752.5 

3714.4 

3154.7 

3721.0 

3003.04 

18.1 68 

U 

65 

0.50 

0,100 

1380,5 

1200.4 

1217.2 

1065.0 

1363. 0 

1336,5 

1329.0 

1284.3 

1133.3 

1400.9 

1272.07 

112 04 

12 

65 

1.00 

0.001 

3704.3 

3732.7 

3004.5 

3594.8 

3743.1 

3686.3 

3747.0 

3730.7 

3741.0 

3733.4 

3710,78 

65 23 

13 

85 

1 00 

0.100 

1177.0 

1102.3 

1244.0 

1107.2 

1228.7 

1185.0 

1144 7 

1253.2 

1200.6 

1180.4 

1 108 09 

35 00 

H 

10 

0 75 

0.001 

3320 4 

3175.7 

3482.3 

3178.9 

2373.8 

3153.1 

.3058.1 

3188.0 

3370,0 

2091.1 

3128,87 

302 4 1 

in 

10 

0.75 

0.100 

675 0 

508.2 

053.0 

661.1 

804.8 

504,7 

580 a 

607.1 

000 8 

01.3.0 

587.42 

30 80 

1(1 

100 

0.75 

0 001 

3852,0 

3848 5 

3830.3 

3841,0 

3861.1 

3846 8 

38.53 7 

3854.0 

3840.3 

38.38.7 

3817.10 

5 82 

17 

UK) 

0,75 

0.1 00 

13(104 

1352 4 

1434.0 

1330.3 

1330.7 

1,308.2 

1302.6 

1337.8 

1208 2 

1418 0 

1,3.58 74 

4(1 10 

IB 

10 

0.50 

0.061 

IIBO.O 

1215.1 

1182.7 

1210.0 

1330.4 

1285.6 

1175.1 

1210.3 

1104.5 

1208.1 

1210.02 

54.09 

10 

10 

1.00 

0 051 

1 100.2 

1227 8 

1357.1 

1107.0 

1448.8 

1119.3 

1201.3 

1251.5 

1245.5 

1220.9 

1246.00 

06,75 

20 

100 

0.60 

0 051 

2590.0 

2035.0 

2004.0 

3620.3 

2070.0 

2037.5 

2032.4 

2632.3 

2048.4 

2042,0 

2637.51 

22 35 

21 

100 

1.00 

0.051 

2030.3 

2002.4 

2054.0 

2031.0 

2720.8 

2025.5 

2020.4 

2578.0 

2029. T 

2670.4 

2042.02 

37,31 

22 

55 

0.76 

0 001 

3014.3 

3732.8 

3737.7 

3748.3 

3160.1 

3735.2 

3730 4 

3748.0 

3748.0 

3735.7 

3000 80 

183 53 

23 

65 

0.75 

0 100 

1207.2 

1130.8 

U79.0 

1194.3 

1064.3 

1171.6 

1184.2 

1159.7 

1170.1 

1180.8 

1170.81 

51 01 

24 

10 

0,75 

0.061 

1443.7 

U08.8 

1212.3 

1333.1 

1281.1 

1234.8 

1101.0 

1279.6 

1090.1 

1384.9 

1256.04 

11283 

2S 

100 

0.76 

0,061 

2003.1 

2585.7 

2666.0 

2052.4 

2040.7 

2054.9 

2038.7 

2661.0 

2548.8 

2033.1 

2618.39 

30 00 

20 

56 

0.60 

0.061 

2637.2 

2460.0 

2604,0 

2000.4 

2619.0 

2491.0 

2541.4 

2023.8 

2527.3 

2570.7 

2547.15 

5 1 no 

27 

65 

1.00 

0.051 

23204 

2342.7 

2465.1 

2370.6 

2420.1 

2458.1 

2382.5 

2446.1 

2388.2 

2308.0 

2300.05 

45.75 


Table 4.32: Data obtained from 27 trials on DeJong F2 function case 4 


The signilituuit elfects of parameters, their interaction and the cjuadiatic ellects aie 


shown in Table 4.33 for average and in Table 4.34 for standard deviation. 







4.2 Test Function 2: De Jong Function F2 




Uesign parameter 

Effect 

MP 

PS 

PS*PS 

MP*MP 

PS*MP 

CP*MP 

CP*CP 

mutation prob. 
population size 
population size (quad) 
mutation prob. (quad) 
population size * mutation prob. 
crossover prob. * mutation prob. 
crossover prob. (quad.) 

-2424.54 

1035.92 

352.3 

-159.48 

-67.55 

-36.2 

18.94 


Inblc 4.33. Effects table for average case 4 


Design parameter 

Effect 

PS*MP 

population size * mutation prob. 

272.6 

PS 

population size 

-191.42 

MP 

mutation prob. 

-182.69 

MP*MP 

mutation prob. (quad.) 

-78.86 

PS*PS 

population size (quad.) 

-36.015 

CP*MP 

crossover prob. * mutation prob. 

30.455 

CP 

crossover prob. 

-28.6 

GP*CP 

crossover prob. (quad.) 

-26.9 


Table 4.34: Effects table for standard deviation case 4 
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Response surface study 


Unsi'd oil 27 trials data, in 'rabh' 4.32, tho best 
resulting equation is : 


surlace through the 


average is fitted. 


Tfie 


on-line performance = 2571.32 + 31.054P5 - 24666.67MP - Q.17PS^ - 134.43C'P^ 

-20733. 15MP2 -f 3625.19C'P x MP - 10.3P5 x MP (4.24) 


Tho best fit surfaeo through standard deviation is: 

log(cr) = -5.2 - 0.016P5 + 30.92CP - 23.89AfP - 0.0002P52 - 19.24C'P^ 

+298.578A/P'^ - 4G.1CP x MP -f- 0.4338P5' x MP (4.25) 

Analysis by direct observation 

• Set PS = 100 , CP = 0.5 and MP = 0.001 in order to minimize variation. 

Analysis by Variation transmission 

Minimization of the eejuation for standard deviation of on-line performance obtained by 
equation 4.24 gives: 

a = 110.9 at PS = 87, CP = 1.0, and MP = 0.001. 

Putting the value PS = 87, CP =1.0 and MP = 0.1, in the equation 4.25 we get 
cr = 31.28. The value of standard deviation obtined by direct observation method is 5.62 . 
This means that result from variation transmission method is not good. 

So with mean at PS = 100, CP = 0.5, and MP = 0.001 and standard deviation 
aps = 5 ,(Tc 7 > = 0.05 and ctmp = -005 final capability study is done to verify the expected 


improvement. 
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4.2.5. 1 Final Capability Study for Case 4 


Tho (lata aie given in lal)lc 4.35. The control chart for average and standard deviation 
is also given in hignie 4.19 and figure 4.20. We see that one point in standard deviation 
control chart is out of the control limit. Also the on-line performance obtained in capability 
study is much more compared to case 1 and case 2. So this means that with the tournament 
s(d(!ction and unilorrn crossover it is not possible to tirake the process capable. 
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Tabhi 4.35: Final capability study for case 4 function F2 




Trial Number 


Figure 4.19: Contr ol cha.rt for average of on- Figure 4.20: Coirtrol chart for standard devi- 
line performance.: after variation reduction ation of on-lirre performance after variation 
case 4 on function F2 reduction case 4 on function F2 
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4.2.6 Discussion of Result of Function F2 


111 the 'Inblc 4.36 tlu^ average aiui I, he standard deviation of tiie final capability study are 


shown for each case. Ihe value of PS, CP and MP are given at which we have to fix it with 


standard deviation of 5, 0.05 and 0.005 respectively, for variation reduction. 


Table 4.36: Result of function F2 



PS 

CP 

MP 

average 

std. deviation 

CASE 1 

Proi>ortionate selection 
single point crossover 

82 

0.5 

0.1 

295.3 

15.38 

CASE 2 

proi)ortionate selection 
uniform crossover 

55 

0.5 

0.1 

295.88 

15.23 

CASE 3 

tournament selection 
single point crossover 

10 

1 

0.1 

624.44 

1 

91.98 

CASE 4 

tcmmaiuent selection 
uniform crossover 

100 

0.5 

0.001 

3801.69 

30.08 


We see that proportionate selection is always better for minimization of both on-line 
performance and variation in on-line performance. We see that if proportionate selection is 
the selection operator, both single point crossover and uniform crossover gives almost same 
value of average and st.andard deviation of capability study. So any of the crossover operator 
can be chosen with proi)ortionate selection as selection operator. Also we have seen that in 
case 3 and case 4, when tournament selection is used it is not possible to make the process 
capable. So for this function proportionate selection is recommended. 

In every case PS, CP and MP which has to be targeted, varies. Since we found that 
proportionate selection with any of the crossover operator is better, we can set PS at 55, CP 
at 0.5 and MP at 0.1 as mean with standard deviation of 5 for PS, 0.05 for CP and 0.005 
for MP with uniform crossover as crossover operator. Otherwise we can set PS at 82, CP at 
0.5 and MP at 0.1 as mean with standard deviation of 5 for PS, 0.05 for CP and 0.005 for 
MP with single point crossover as crossover operator. 
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4.3 Test Function 3 : Rastrigin’s Function 

4.3.1 Initial Capability Study 

To check the need tor vaiiation reduction capability study is done. The on-line performance 
measured for last oO A) of the geneiation are shown in Table 4.3T. The average and the 
standard deviation of each trial are also shown in Table 4.37. 
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Table 4.37; Initial capability study for Rastrigin’s function 


These data are used to construct the average and standard deviation control charts shown 
ill Figure 4.21 and Figure 4.22. 

We can see from Figure 4.21 and Figure 4.22 that lot of points are going outside the control 
limit. So there is need to reduce the variation. For reducing the variation four cases as 
described in previous chapter is selected and for each case variation reduction technique is 
applied. 

4.3.2 CASE 1 (Proportionate Selection and Single Point Crossover) 

27 trial oxixuimonl, (la.ta are sliown in Table 4.38. 

The vSignilicaut olfocts of parameters, their interaction and the quadratic effects are 
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Figure 4.21: Control chart for average on- Figure 4.22: Control chart for standard devi- 
Imc porforn.ana. I.doro var.ptio.. ro,luc.tio„ ation of performance before variation 

o„ Haslngm s lunol.on reduction on Itetrigin’s function 
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Table 4.38: Data obtained from 27 trials on Rastrigin’s function case 1 
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sliowii ill Tfibic 4.39 ioi a,veia,ge and in Table 4.40 for standard deviation. 


Design paraiucter 

Effect 

MP 

mutation prob. 

184.37 

MP*MP 

mutation prob. (quad) 

71.992 

PS 

population size 

68.08 

PS*MP 

population size * mutation prob. 

-27.55 

PS*CP 

population size * crossover prob. 

7.8 

CP*MP 

crossover prob.* mutation prob. 

-3.32 


Table 4.39: Eflects table for average case 1 



Dc'sign para.nictor 

Effect 

MP 

mutation prob. 

-5.21 

I’S 

population size 

-4.23 

PS*CP 

population size*crossover prob. 

-3.04 

MP*MP 

mutation prob. (quad.) 

3.02 

PS*PS 

population size (quad.) 

-2.441 

CP*CP 

crossover prob. (quad.) 

2.14 

CP 

crossover prob. 

-2.02 

CP*MP 

crossover prob.* mutation prob. 

1.71 


Table 4.40: Effects table for standard deviation case 1 


Response surface study 

Based on 27 trials data in Table 4.38, best surface through the average is fitted. The resulting 
equation is : 


on-line performance = 130.16 + O.dOSSP^T 3036. 84MP 10881. 16MP 

-0.01P5 X CP + 200.76CP x MP - 1.897P5 x MP (4.26) 

The best fit surface through standard deviation is: 

log(,r) = 4.920 - omnrs + l.aiCP - 2 .IIMP + O.OOOOOSPS* - 1.446CP’‘ 

-Sl.lSAfP" - O.OOOOSSPS' X CP + 7.622CP x MP (4.27) 
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Analysis by direct observation 

The result from this analysis can be concluded as: 


Set PS = 250, CP = 0.5 and MP = 0 1 in nrHpr to • • • • ■ 

mi - u.i in order to minimize variation. 


Analysis by Variation transmission 

Mmimizatioii of equation for stadard deviation, obtained from equation 4.26 of on-line per- 

fonnaiice givi's: 

a = 2.83 at PS = 250, CP = 0.5, and MP = 0.1. 

I his means that result Irom analysis of variation transmission method and by direct 
observation method lead to same value of PS, CP and MP. Now the transmitted variation 
of each paramettn- was found and given as: 

atps = 1.003, fficp = 0.84 and atMP = 2.46 

This shows that more control on variation of MP is required than PS and CP. So there 
is need to reduce the variation of MP i.c. MP is a VIP. 

So with mean at PS = 250, CP = 0.5 and MP = 0.1 and standard deviation aps — 5, 
Ocp ~ .05, and cyi^ip = 0.()05 liual capability study is done to verify the expected improve- 
ment. 


4.3.2. 1 Final Capability Study for Case 1 

The data are given in Table 4.41. The control chart for average and the standard deviation 
is also given in Figure 4.23 and Figure 4.24. We see that process is capable as all points are 
within the control limil.s. Now we try next case. 

4.3.3 CASE 2 (Proportionate Selection and Uniform Crossover) 


27 trial data, ai(' shown in 'ra.l)le 4.42. 
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Figure 4.23: Control chart for average of on- Figure 4.24: Control chart for standard devi- 
line performance after variation reduction ation of on-line performance after variation 
for case 1 on Rastrigin’s function reduction for case 1 on Rastrigin s function 
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369.7 

381.8 

298.1 

325.9 

376.9 

381.6 

280.6 

379.9 

318.3 
379.9 

364.4 
376.7 


381.6 

359.1 

341.0 

354.4 

147.5 

344.2 

383.1 

152.3 

372.3 

287.6 

376.1 

300.6 

379.2 

150.8 

329.6 

360.9 

382.5 

304.3 

320.5 

376.9 

382.7 

332.3 

379.2 

309.2 

379.4 

370.7 

376.4 


r 

381.8 
367.7 

313.3 

347.9 
146.1 

354.0 

385.3 

142.3 

378.6 

294.7 

376.7 

322.1 

378.9 

148.1 

318.0 

357.6 
382.5 

273.1 

318.7 

374.1 

381.3 

295.0 

377.9 

200.1 

379.7 
357.1 

375.9 


”W 


381.9 

367.1 

338.4 

358.9 
148.3 

348.7 

382.5 

150.8 

372.2 

308.9 

378.6 

306.0 

380.9 

140.0 

326.9 

359.0 

382.8 

298.8 

325.3 

374.6 

380.4 
310.2 

378.4 

307.0 

380.1 

368.1 

377.1 


on-line perlorroance" 


average 


377.78 

362.84 

333.33 
344.08 
143.42 
348.21 

379.96 
147.63 
369.25 
286.03 

372.20 

302.91 
376.12 
155.82 

331.96 

355.34 

378.21 
290.45 
321.94 

371.91 
378.06 
302.65 
375.28 
305.71 

376.92 
362.90 
372.86 


»td, dev. 


6.66 

6.98 
9.3 T 

15.76 

4.48 
4.81 
6.13 

2.99 
6.24 

15.29 

5.48 

10.72 
5.58 

16.73 
8.67 
9.63 
6.29 

10.12 

10.68 

5.48 

5.76 

15.48 

4.80 

13.74 
6.02 
6.10 
6.07 


lablc 4.42: Data obtained Iroin 27 trials on Rastrigin’s function case 2 

The significant effects of parameters, their interaction and the quadratic effects are 
shown in Table 4,43 for average and in Table 4.44 for standard deviation. 


Design parameter 

Effect 

PS 

population size 

105.1 

MP 

mutation prob. 

96.82 

PS*MP 

population size * mutation prob. 

-82.15 

MP*MP 

mutation prob. (quad) 

34.62 

PS*PS 

population size ((piad) 

29.98 

CP 

crossover prob. 

-11.01 


Tabh' 4.43: Flft'c.ts tjd)le for av('rage case 2 


Response surface study 

Based on 27 trials data in Table 4.42, the best surface through the average is fitted. The 
resulting equation is : 

on-line perfonnauce = 118.2 + 1.49P5' + 23.25CP ■+ 2259.16MP - 0.00239PP 

-2918.11MP2 - 7.235P5 x MP (4.28) 
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Design paraniotor 

Effect 

MP 

mutation prob. 

-4.57 

PS*MP 

population size * mutation prob. 

-2.15 

CP*CP 

crossover prob. (quad.) 

2.1 

CP 

crossover prob. 

-1.67 

PS 

population size 

-1.54 

CP*MP 

crossover prob. * mutation nrob. 

1.31 


4.44. fcjfiects table for standard deviation case 2 


The Ix'st lit surfiux^ 


(.h rough standard deviation is: 


log(fT) 


-0.4438 + O.OOlli’S' + 7X)CP -f 7.569MP - 4.53C'P^ - 9.55CP x MP 
-().{)33G/\9 X MP (A n 


Analysis by direct observation 


• Set PS = 250 , CP = 1 and MP = 0.1 in order to minimize variation. 

Analysis by Variation transmission 

Minimizing the cciuation for standard deviation obtined from equation 4.28 gives: 

a == 1.45 at P6' = 250, CP = 0.G8 and MP = 0.055. 

'llie. standard d('viation from (Kunvtion 4.29 will bo 8.114 at PS = 250, CP = 0.68 and 
MP = 0.055. The stajidard deviation obtained from equation 4.29 by direct observation- 
method is 3.61. Tins means that the result from Variation transmission method method is 
not good. 

So with mean at PS = 250, CP = 1 and MP = 0.1 and standard deviation aps = 5, 
<ycp = .05, and (Tmp .005 linal capability study is done to verify the expected improvement. 
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4.3.3. 1 Final Capability Study for case 2 


IS 


Tlic aro givan in Tal.lo .1.45, The control cliart for average and standard deviation i 
also given in I" igiirc 4.25 and I' igiire 4.26. We see that all points are within the control limits. 

Hence the process is capable. Now we try for next case. 


TVuvl 

number 

1 

2 



^ 

Hun n 
5 

umber 
5 

7 

5 


TO 

average 

ltd. 

2 

3 

4 

5 
(1 

7 

S 

1) 

10 

1 1 

12 

13 

14 

15 

HI 

IT 

la 

H) 

20 

370.0 

372.2 

370.0 

300.0 

aoa.o 
300, a 
aoa.a 
aoa.a 

370.2 

370.1 
307.0 
300.5 
3f50.1 

300.3 

300.4 

307.5 

307.3 

370.3 
30a.5 

371.2 

371.1 

370.2 

371.2 

371.4 

300.0 

300.5 

371.0 

300.8 

370.0 

370.8 

370.8 

308.0 

300.5 

300.0 

300.3 
300.7 

300.0 

371.0 

374.7 

374.4 
375. a 

375.8 

375.0 

375.1 

373.8 

375.0 

373.5 

374.2 
374.2 

370.0 

375.1 

370.2 
374.4 

373.3 

373.4 

375.8 

375.4 

382.8 

381.5 
380.0 

380.0 

380.5 

370.0 

380.0 

381.7 

382.0 

381.8 

380.3 

382.8 

383.3 

381.1 

380.1 

382.0 

381.1 

380.4 

383.5 
380.4 

383. f 

383.2 

382.3 

381.3 
383.6 

383.1 

384.1 

382.3 
383.5 

383.0 

385.1 
383.0 

383.0 

385.3 

382.1 
382.9 

383.8 

381.9 

382.2 
383.0 

383.3 

383.7 

382.4 

381.7 

382.5 
383.9 

384.2 

384.0 

383.2 

385.0 

386.3 

384.1 
383.9 

384.7 

384.2 

383.4 
v383.3 

383.4 

384.4 

383.7 

384.1 

383.1 

382.4 

383.6 

382.5 

383.7 
382.9 
'382.9 

382.9 

383.2 

384.7 

385.4 

383.7 

382.3 

382.6 

382.5 
384.2 

381.9 

382.9 

384.6 

383.6 

384.5 
382.9 

382.5 

384.0 

383.2 

384.5 

383.3 

385.0 

384.3 

385.0 
382.9 

381.1 

382.4 

381.8 

384.2 

383.3 
382.1 

381.3 

382.9 

383.0 

383.3 

382.4 
384.7 

384.9 

386.0 

383.4 

382.2 

382.6 

3.82.9 j 

382.5 ! 

383.7 i 

383.1 

383.6 

383.5 

384.9 

383.3 
384.0 

380.6 

384.4 

384.2 

382.9 

381.9 

382.9 

383.2 

383.6 

384.7 

384.3 
383.6 

383.6 

383.7 

383.3 

385.4 

384.4 

380.0 

383.7 

382.7 
381,6 

384.8 

382.1 

380.2 
379.9 

379.2 

379.2 
379.7 

379.7 
379.9 

379.3 
379.9 

379.8 

380.2 

379.9 

380.0 

379.7 

379.0 

379.8 

379.2 

378.6 

379.6 

379.7 

5.54 

5.48 

4.70 
6.40 

6.70 
6.04 
6.03 
6.13 
5.87 
6.08 
6.37 
6.44 
6.81 
0.32 
5.50 
6.15 
0,50 
6.07 
5.61 
5.79 


Tal)l('. 4.45: Final capability study 


for case 2 Rastrigin’s function 




Figure 4.25: Control chart for average of on- Figure 4.26: Control chart for standard devi- 
line performance after variation reduction ation of on-line performance after variation 
case 2 on R.astrigin’s function reduction case 2 on Rastrigin’s function 


4.3.4 CASE 3 (Tournament Selection and Single Point Crossover) 

27 trial experinxmt data are shown in Table 4.46. 
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"Trial 

no. 

cai 

"nTT" 

Tilda to 
variabl* 

pnr"'”"’ 

npiifc 

_ "-j— 

""250 

d.nd 

O.lOO 

2 

250 

l.OO 

0.001 

3 

25 

0.50 

0.100 

4 

250 

0.50 

0.001 

6 

25 

1.00 

0.001 

6 

25 

1 .00 

0.100 

7 

250 

1 .00 

0.100 

B 

25 

0.50 

0.001 

9 

137 

0.75 

0.061 

10 

137 

0.60 

0.001 

11 

137 

0.50 

0.100 

12 

137 

1.00 

0.001 

13 

137 

1.00 

0.100 

14 

25 

0.75 

0.001 

15 

25 

0.75 

0.100 

16 

250 

0.75 

0,001 

17 

250 

0.75 

0.100 

18 

25 

0.50 

0.051 

19 

25 

1 .00 

0,051 

20 

250 

0.50 

0.051 

21 

250 

1.00 

! ().()5l 

22 

137 

0.75 

0,001 

23 

137 

0.75 

0.100 

24 

25 

0.75 

0,051 

2 IS 

250 

0.75 

0,051 

26 

137 

0.50 

o.oni 

27 

137 

1.00 

0,051 


3 71). a 

073.4 
077.9 

383.4 

404.4 

000.2 

443.3 

733.0 
401.7 

761.1 

402.4 

644.2 

381.0 

098. 2 

405.3 

413.6 

400.6 

439.1 

443.7 

740.3 
402.6 

400.3 
44.6,8 

444.5 

442.2 


3 

Too. 9“ 

706.0 

382.1 

093.6 

053.7 

381.0 

405.0 
039.6 

446.8 

740.5 

401.8 

759.0 

400.0 
658.4 

378.3 

703.0 

404.4 

417.4 

411.6 

448.6 

444.6 

748.1 

400.1 
405.3 
447.0 

430.5 

442.7 


T 

408.9 

720.4 

387.1 

092.7 

677.7 

391.1 

412.6 

663.1 

443.9 

705.2 

407.5 

762.8 
405.4 

657.7 

385.3 

712.8 

411.1 

421.1 

420.1 

450.9 

451.4 

770.3 
404.7 

419.3 

452.2 

447.3 
448.1 


410.8 

719.5 

387.3 

726.3 

660.9 

388.8 

417.8 

654.3 

454.7 

746.9 

411.2 

757.8 

410.9 

663.6 

389.4 

735.2 

417.2 

423.6 

421.7 

462.4 

460.8 

767.0 

411.0 

419.8 

455.1 

451.4 

455.4 


Run nuniljnf 


418.2 

745.6 

389.2 

721.0 

695.2 

391.2 

419.4 

662.6 

460.5 

769.5 

421.0 

783.5 

415.8 

648.4 

389.3 

716.9 

417.9 

425.0 

439.9 

460.9 
464.2 

783.6 

416.1 

423.9 

460.1 

457.5 
458.0 


r 

417.7 

741.2 

385.8 

706.2 

599.1 

388.5 

416.6 

694.4 

456.9 

742.5 

416.3 

783.2 

416.1 

631.6 

387.7 
721.6 
419.5 

419.2 

418.3 

458.0 

460.8 

756.1 

414.8 

422.2 

401.9 

458.4 

453.3 


418.2 

736.5 

390.8 

716.2 

609.0 

391.5 

420.6 

678.9 

455.9 

756.5 

413.9 

775.2 

413.5 

623.7 

384.4 

736.9 

419.5 

420.5 

419.1 

459.7 

463.6 

766.6 
412.4 

421.7 

461.8 

455.8 
455.7 


TT 


Tio.o 

739.2 

388.2 

711.2 

633.0 

385.0 

418.5 

606.3 

460.5 

756.2 

413.7 

789.1 

416.8 

647.8 

387.6 

714.4 

417.0 

419.4 

420.0 

460.8 

460.5 

754.7 

418.3 

416.7 

400.8 

459.4 

455.4 


418.6 

742.6 

394.3 

722.2 
673.9 
391.5 

421.4 

676.1 

459.8 

766.3 

414.5 

763.9 

418.4 

650.0 

392.0 

739.2 

417.4 

424.2 

424.8 

460.8 

458.7 

775.8 

414.2 

423.4 
463.7 

457.4 

461.2 


— nr 

421.0 

733.4 

385.8 

708.2 
621.6 

388.5 

416.3 

616.7 

454.0 

758.4 

415.1 

755.7 

416.5 
602.3 

388.5 

730.5 

416.9 

416.7 

419.2 

460.8 
461.0 
769.7 

416.2 

418.6 

458.3 

451.4 

456.6 


on-line performance 


average 

414.93 
730.45 

386.98 
707.09 

650.22 
388.17 
415.32 

636.22 
453.63 
747.48 
411.67 
769.14 
411.58 
642.77 
380.35 
720.86 
414,62 
420.05 
420.12 
456.19 

456.93 
762.83 

410.92 
417.72 
456.75 

451.98 

452.93 


. dev. 
5.92 
13.23 

4.24 
16.41 
32.77 

3,54 

6.97 

33,70 

6.66 

18.80 

6.24 
12.40 

6.64 
18.68 

4.10 

14.33 

5.65 
3.63 
8.68 
7.60 
7.56 

12.25 

6.29 

6.66 
6.29 
7.38 
6.47 


'l'a,bki 4.46: Data obtained from 27 trials on Rastrigin’s function ca^e 3 


The signilicant efie(!t„s of the parameters, their interaction and the quadratic effects are 
shown in Table 4.47 for av(u-age and in Table 4.48 for standard deviation. 


Design parameter 

Effect 

MP 

mutation prob. 

-302.83 

MP*MP 

mutation prob. ((piad) 

-113.003 

I’S 

population size 

47.28 

PS*PS 

population size (quad) 

34.64 

PS*MP 

population size * mutation prob. 

-24.46 

cn>*MP 

crossover prob. * nmtation prob. 

-9.75 

CP 

cros.sover prob. 

7.029 


I'able 4.47; Elfccts table for average case 3 


Response surface study 

Based on 27 trials data in 'ra,ble 4.4G, the best surface through the average is fitted. The 
resulting equation is : 

on-line performance = 543. C + 1.13P5' + 85.36CP — 3558. 68MP — 0.00274P5 

+2()9G2.92A/P'-' - 1447.79CP xMP- 3.35P5 x MP (4.30) 
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Table 4.48: Effects table for standard deviation case 3 



besign parameter 

Effect 

MP 

mutation prob. 

-13.78 

PS*MP 

population size * mutation prob. 

7.797 

MP*MP 

mutation prob. (quad) 

-5.52 

PS 

population size 

-3.67 

CP*CP 

crossover prob. (quad.) 

-2.18 

PS*PS 

population size (quad.) 

-1.82 

CP*MP 

crossover prob. * mutation prob. 

1.71 


The best lit surfaces through standard deviation is: 


log(o) 


0.1)02-12/ ’,b’ 2;!.!)G/l//' - O.OGiG’/''-^ 4- A[A2MP'^ 


+2.mCP X MP + 0.043P5 x MP 


(4.31) 


Analysis by direct observation 


• Set PS = 25 , CP = 0.5 and MP = 0.1 in order to minimize variation. 

Analysis by Variation transmission 

Minimization of e(iua.tion for standard devaition obtained from equation 4.30 gives: 

O' = 3.05 at PS = 98, CP = 0.5 and MP = 0.1. 

Putting the value PS = 98, CP = 0.5 and MP = 0.1 in the equation 4.31 of standard 
deviation of on-liiui performance we get a - 4.29. The value of standard deviation obtained 
by direct observation method is 3.72. This means that the result from variation transmission 
method is good. 

Tlie transmitted va.riation of the parameter is: 

(^i.ps “ ~ 4.25, (7tMP ~ 1-69 


We see that traiismitted variation by CP is largest, hence CP is a VIP. So with mean 
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at PS = 98, CP = 0-5, and MP = 0.1 and standard deviation cps = 5 and ctmp = -005 final 
(capability study is done to verify the cx})ected imi)roveincnt. 

4.3.4. 1 Final Capability Study for Case 3 

The data are given in Table 4.49. The control chart for average and standard deviation is 
also given in Figure 4.27 and Figure 4.28 . We see that process is capable as all points fall 
within the control limits. Now we try for next case. 


'IVinl 

number 

Run number | 

avemge 

"TuT” 

dflV. 

1 

2 

3 

4 

L_ S 


7 

g 

9 

lo 

1 

400.5 

400,1 

407.7 

411.4 

410.9 

ll3.5 

413.8 



Ui.i 



2 

400.0 

402.9 

406.4 

414.2 

410.2 

415.4 

413.5 

414.0 

417.4 

412.3 

411.2 

5.90 

3 

396.3 

402.5 

407.5 

412.5 

412.7 

410.5 

417.9 

412.4 

412,9 

411.3 

409.7 

6.15 

4 

401.6 

398.5 

404.3 

412.5 

412.0 

410.7 

413.1 

411.0 

410.1 

4U.4 

408.5 

5.14 

5 

398.6 

403.0 

410.9 

411.4 

416,3 

410.8 

413.0 

417.6 

410.8 

410.5 

410.2 

5.66 

6 

396.3 

399.7 

404.5 

411.9 

411.2 

410.9 

415.6 

412.1 

412.5 

411.9 

408.7 

6.31 

7 

401.7 

400.7 

406.6 

412.5 

414.7 

410.7 

421.5 

413.3 

409.7 

408.9 

410.0 

6.15 

8 

395.1 

396.3 

400.7 

410.1 

408.1 

411.2 

408.6 

412.9 

413.0 

412.4 

406.8 

6.88 

9 

402.2 

402.8 

409.4 

411.5 

416.6 

413.3 

413.9 

415.7 

414.4 

414.0 

411.4 

5.10 

10 

397.1 

401.8 

403.6 

408.6 

406.2 

412.1 

412.7 

411.1 

411.4 

413.1 

407.8 

5.43 

11 

397.3 

400.8 

403.0 

410.9 

413.3 

412.9 

414.1 

410.2 

411.6 

411.2 

408.5 

5.90 

12 

396.9 

398.9 

404.7 

411.1 

414.4 

414.6 

412.8 

412.0 

416.1 

411.4 

409-3 

6.75 

13 

401.2 

400.7 

408.0 

414.5 

413.5 

414.9 

415.1 

416.2 

414.8 

414.4 

4U.3 

5.92 

14 

399.7 

400.2 

404.8 

410.6 

413,9 

416.5 

413-7 

412.0 

418.7 

411.4 

410.1 

6.50 

15 

395.0 

399.4 

401.0 

410.8 

411.5 

412.6 1 

408.8 

412.3 

412.4 

410.4 

407.4 

6.44 

16 

399.2 

404.1 

403.4 

409.6 

410.0 

413.2 

409.9 

411.9 

412.8 

413.9 

408.8 

4.90 

17 i 

399.4 

401.3 

408.0 

410.5 

413.5 

414.9 ! 

411.2 

418.5 

413.7 

411.7 

410.3 

5.94 

18 ' 

397.6 

398.9 

401.6 

407.8 

411.0 

410.2 

409.7 

408.8 

409.0 

411.5 

406.6 

5.20 

19 

402.1 

402.6 

406.3 

411.2 

413.8 

417.1 

414.3 

413.9 

413.5 

415.3 

411.0 

6.38 

20 

400.5 

400.6 

409.4 

413.8 

415.0 

416.8 

416-8 

417.0 

416.3 

416.9 

412.3 

6.62 


Table 4.49: Final capability study for case 3 Rastrigin’s function 




Figure 4.27: Control chart for average of on- Figure 4.28: Control chart for standard devi- 
line performance after variation reduction ation of on-line performance after variation 
case 3 on Rastrigin’s function reduction case 3 on Rastrigin’s function 
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4.3.5 CASE 4 (Tournament Selection and Uniform Crossover) 

27 trials experiment data are shown in Table 4.50. 


imQii 

candidate input 
variables 

Run number 


mn 


UF 




3 


5 - 


7^ 


. A 




■■Bl: 


0.50 

0.100 

390.1 

391.7 

395.4 

“403.2 

401.4 

405.0 

404.2 


K 1 tlCTi 

404.0 

400.17 

5.57 i 

2 


1.00 


BtaiKS 

467.5 

452,6 

516.0 

486.9 

■.ItKPtM 

506.1 


K : sfiiK 

499.9 

486.89 


3 


0.50 

0.100 

372.7 

369.7 

372,4 

377.1 

376.7 

377.3 

375.3 


K * 

376.9 

375.62 

3.17 

4 

250 

0.50 

0.001 

489.1 

506.1 

507.2 

484.9 

520.6 


556.8 

487.2 

M ^ mIm 

511.3 

505.19 

22.74 

5 

25 

1.00 

0.001 

G01.9 

581.8 

026.6 

663.5 

650.9 

636.2 

635.0 

641.2 

636.0 

617.2 

629.03 

23.72 

6 

25 

1.00 

0.100 

366.0 

367.4 

371.2 

370,7 

371.2 

368.3 

368.0 

369.3 

374.3 

372.0 

369.83 

2.50 

7 

250 

1.00 

0.100 

380.5 

380.6 

385.3 

393.1 

394.1 

392.9 

393.9 

393.1 

393.3 

393.8 

390.06 

5.64 

8 

25 

0.50 

0.001 

543.0 

605.3 

643.7 


633.9 


627.8 


678.0 

627.0 

629.93 

36.19 

9 

137 

0.75 


388.7 

391.7 

397.2 

402.4 

406.5 

404.2 

407.5 

405.2 

406.2 

407.4 

401.69 

6.82 

10 

137 

0.50 



658.2 

673.7 


706.6 

742.8 

mMtwm 

718.7 

665.2 

702.2 

685.30 

32.18 

11 

137 

0.50 

0.100 

386.2 

389.0 

390.7 

398.5 

399.8 

399.1 

399.4 

400.2 

399.6 

401.2 

396.37 

5.49 

12 

137 

1.00 

0.001 

641.8 

654.7 

662.0 

688.0 

696.0 

692.5 

657.5 

683.5 

682.6 

664.5 

672.32 

18.49 

13 

137 

1.00 


377.8 

376.8 

383.9 

391.6 

390.9 

393.1 

391.2 

391.1 

391.4 

391.6 

387.94 

6.12 

14 

25 

0.75 

0.001 

657.9 

648.9 

655.1 

643.5 

651.9 

615.0 

611.9 

625.4 

646.3 

610.3 

636.62 

18.89 

15 

25 

0.75 

0.100 

369. 0 

367.9 

373.6 

374.3 

372.7 

371.3 

373.4 

373.9 

381.6 

371.1 

372.88 

3.72 

16 

250 

0.75 


491.1 

462.6 

471.7 i 

474.4 

480.5 

507.7 

491.0 

489.0 

499.3 

532.1 

489.95 

20.00 

17 

250 

0.75 

0.100 

386.3 

382.8 

388.2 1 

398.7 

397.9 

398.0 

398.0 

397.7 

398.9 

396.0 

394.25 

6.04 

18 

25 

0.50 

0.051 

383.4 

388.1 

396. C 

391.2 

397.5 

394.5 

398.9 

396.3 

399.9 

394.4 

394.08 

5.15 

19 

25 

1.00 

0.051 

373.0 

380.9 

386.5 

389.9 

384.1 

380.9 

387.6 

388.5 

394.2 

384.2 

385.58 

5.69 

20 i 

250 

0.50 

0.051 

404.0 

405.1 

405.1 i 

415.5 

414,4 

418.9 

419.9 

415.2 

416.0 

417.9 

413.27 

6.99 

21 i 

250 

l.OO 

0.051 

380.3 

385.5 

393.0 

399.5 

401.2 

402.3 

399 6 

401.6 

401.9 

398.1 

396.90 

6.39 

22 

137 

0,75 

0 001 

nnn o 

0(12 5 

676 4 ^ 

6n6,4 

675 .3 

669 7 

667 8 

6M7.7 

666 6 

669 8 

672,28 

12 28 

23 

137 . 

0.75 

0.1 00 

3.M3 0 

381.3 

3H7,5 

394.2 

393.0 

391 6 

304 0 

395.0 

307.1 

393.0 

301 81 

8,44 

24 

25 

0.75 

0.051 

383.1 

381.0 

387.7 

393.0 

393.5 

387.6 

389.6 

390.4 

400.5 

386.5 

389.29 

5.57 

25 

250 

0.75 

0.051 

392.9 

392.7 

399.1 i 

408.0 

406.7 

406.4 

409.1 

408.0 

408.8 

407.4 

403.90 

6.48 

26 

137 

0.50 

0.051 

403.0 

400.0 

407.3 ‘ 

412.1 

414.8 

415.8 

414.1 

414.1 

416.1 

411-4 

410.85 

5.61 

27 

137 

mm 

0.051 

387.0 

387.2 

392.4 

396.2 

398.4 

399.7 

401.9 

402.8 i 

402.1 

398.9 

396.65 

5.91 


Table 4.50: Data obtained from 


27 trials on Rastrigin’s function case 4 


The significant 
shown in Table 


effects of the parameters, their interactions and the quadratic effects are 
4.51 for average and in Table 4.52 for standard deviation. 


Design parameter 

Effect 

MP 

mutation prob. 

-214.32 

MP*MP 

mutation prob. (quad) 

-94.54 

PS*MP 

population size * mutation prob. 

79.95 

PS*PS 

population size ((piad) 

42.576 

PS 

population size 

-33.58 

CP 

crossover prob. 

-10.62 

PS*CP 

population size * crossover prob. 

-4.933 


Table 4.51: Effects table for average case 4 


Response surface study 

Based on 27 trials data in Table 4.50, the best surface through the average is fitted. The 
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Table 4.52: Effects tabic for standard deviation case 4 


Design parameter 

Effect 

MP 

mutation prob. 

-17.98 

MP*MP 

mutation prob. (quad) 

-7.887 

CP*MP 

crossover prob. * mutation prob. 

4.656 

PS*MP 

population size * mutation prob. 

3.82 

CP 

crossover prob. 

-2.96 

CP*CP 

crossover prob.(qiiad.) 

-2.61 

PS*CP 

population size * crossover prob. 

1.887 


vesultiug C'xiuatiou is : 

on-line perfonnancM' = G2G.28 -t- 0.474C/’5 - IS.37 MP - 4424.5MP - 0.003344 PS''' 

+14031. 33MP2 - 0.041PS x CP + 6.45PS x MP (4.32) 

The best fit surface through standard deviation is: 

log(a) = 6.127 - 7.64CP-38.98MP + 4.383CP2 + 53.92MP2 
+18.36CP X MP + 0.03PS x MP 

Analysis by direct observation 

• Set PS = 43 , CP = 0.66 and MP = 0.1 in order to minimize variation. 

Analysis by Variation transmission 

Minimization of the equation for standard deviation obtained from equation 4.32 gives: 

(T = 2.47 at PS = 204, CP = 0.54, and MP = 0.1 

Putting the value PS = 204, CP =0.54 and MP = 0.1 in the equation 4.33 of standard 
deviation of on-line performance we get a = 4.27.The value of standard deviation obtained 
by direct observation method is 2.6. This means that result from variation transmission 
method is not good. 


- 0.0007P5 X CP 
(4.33) 
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So with mean at PS — 43, CP — 0.66, and MP = 0.1 and standard deviation ops — 5) 
ffcP ~ — d-()05 final capability study is done to verify the expected improve- 


ment. 

4,3. 5.1 Final Capability Study for Case 4 

The data are given in Table 4.53. The control chart for average and standard deviation is 
also given in Figure 4.29 and Figure 4.30 . We see that process is capable as all points are 
within the control limits except one point in Figure 4.30, which is slightly above the UCL. 


IVial 

number 

Run number j 

average 

std. 

dev. 

j 

2 

3 

4 

5 


1 

1 

g 

— ns — 

1 

372.8 

374.2 

380.4 

379.3 

382.2 

382.0 

383.2 

380.2 

382.3 

380.0 

379.8 

3.54 

2 

370.4 

374.8 

378,9 

378.1 

380.9 

382.5 

383.7 

380.4 

378.8 

378.5 

379.3 

2.70 

3 

371.9 

307.4 

370.4 

382.4 

393.2 

381.8 

385.1 

394.0 

371.8 

386.9 

381.1 

0.20 

A 

309.2 

371.1 

380.0 

377.0 

373.5 

371.7 

378.4 

373.4 

378.8 

376.0 

374.9 

3,66 

5 

370.3 

308.2 

308.1 

379.4 

378.2 

374.3 

377.2 

381.3 

370.8 

376.7 

375.0 

4.65 

0 

371.0 

372.6 

370.2 

379.4 

378.0 

379.6 

378.4 

376.7 

383.0 

377.2 

377.2 

3.48 

7 

371.9 

370.3 

381.7 

388.7 

388,6 

384.8 

387.9 

372.8 

387.8 

380.0 

382.0 

7.49 

8 

374.2 

373.4 

381.9 

380.8 

381.2 

371.2 

387.9 

367.5 

386.2 

386.8 

379.0 

7.42 

n 

370.0 

374.8 

370.0 

380.2 

380.7 

380.0 

381.2 

382.0 

379.8 

380.6 

379.3 

4.45 

10 

372,9 

374.1 

377.2 

378.8 

381.4 

381.1 

377.5 

379.0 

382.1 

379.1 

378.3 

3.03 

11 

374.0 

375.9 

377.2 

382.1 

384.1 

383.6 

384.5 

385.1 

383.7 

382.9 

381.3 

4.00 

12 

374.4 

376.3 

380. 7 

381.3 ! 

384.2 

382.5 

380.8 

384.0 

383.2 

382.3 1 

381.0 

3.22 

13 

373.5 

373.7 

378.2 

381.3 

384.9 

382.5 

381.4 

379.5 

330.1 

381.3 

379.6 

3.64 

14 

373.7 

373.0 

378.9 

382.1 

382.6 

378.6 

381.1 

380.3 

378.3 

379.7 

378.8 

3.21 

16 

373.3 

307.3 

370.6 

370.1 

376.6 

369.2 

383.6 

370.7 

385.6 

384.3 

376.3 

6.46 

10 

371.9 

372.3 

370.0 

378.0 

376.3 

378.5 

380.8 

378.7 

377.2 

377.0 

376.7 

2.76 

17 

371.8 

372.0 

370.9 

379.2 

381.2 

376.4 

378.0 

377.8 

383.5 

376.2 

377.4 

3.53 

18 

373.3 

305.3 

374.2 

378.4 

379.5 

375.0 

371.3 

377.6 

385.0 

383.5 

376.3 

5.82 

19 

370.0 

375.0 

376.5 

379.3 

379.4 

384.1 

378.5 

380.4 

379.0 

380.3 

378.3 

3.63 

20 

372.9 

374.2 

378.4 

380.0 

379.8 

379.5 

380.1 

381.4 

379.2 

379.8 

378.5 

2.75 


Table 4.53: Final capability study for case 4 Rastrigin’s function 



Figure 4.29: Control chart for average of on- 
line perfonnanc.e alter variation reduction 
case 4 on Rastrigin’s function 


Figure 4.30: Control chart for standard devi- 
ation of on-line performance after variation 
reduction case 4 on Rastrigin’s function 
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4.3.6 Discussion of Resuit of Rastrigin’s Function 

In the Table 4.54 the average and the standard deviation of the final capability are shown for 
each case. The value of PS, CP and MP are given at which we have to fix it with standard 
deviation of 5, 0.05 and 0.005 respectively. 

Table 4.54: Result of Rastrigin’s function 



PS 

CP 



std. deviation 

CASE 1 

Proportionate selection 
single point crossover 

250 

0.5 

0.1 

371.49 

6.057 

CASE 2 

proportionate selection 
tmiform cro.ssover 

250 

1 

0.1 

379.01 

5.89 

CASE 3 

tournament selection 
single point crossover 

im 

0.5 

0.1 

409.2 

5.845 

CASE 4 

tournament selection 
uniform crossover 

■ 

0.G6 

0.1 

378.54 

4.78 


We see that if our aim is to minimize only on-line performance then we will choose 
proportionate selection as selection operator and single point crossover as crossover operator 
(case 1). But if our aim is to minimize only variation in on-line performance then tournament 
selection with uniform crossover (case 4) will be chosen, but we have seen that one point is 
outside the control limit, hence this case can not be recommended. The next case in which 
standard deviation is less is case 3, but in this case average on-line performance is high. The 
next higher standard deviation is case 2, so for Rastrigin’s function proportionate selection 
with uniform crossover are recommended. 

In every case PS, CP and MP which are to be targeted, vary. Since we found that 
proportionate selection with uniform crossover is better, we will set PS at 250, CP at 1.0 
and MP at 0.1 as mean with standard deviation of 5 for PS, 0.05 for CP and 0.005 for MP. 

To show that the process capability can not be achieved at the set up of parameters 
corresponding to minimum value of on-line performance, we are taking one set up from initial 
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capability study, corresponding to minimum on-line performance. The minimum value is 
achieved in initial capability stiuly for proportionate selection, uniform crossover, PS = 61, 
CP = 0.9969, and MP = 0.0022. 

Now with mean at PS = 61, CP = 0.9969, and MP = 0.0022 and standard deviation 
aps = 5, acp = 0.05 and <7mp = 0.005 capability study is done. The data are shown in 
Table 4.55, and the control chart for average is shown in Figure 4.31 and control chart for 
standard deviation is shown in Figure 4.32. We see that lot of points are going outside the 
control limits, hence the process is not capable. 


T ’ ria ! 

niiiuber 




I 

2 — 

3 


5 

6 



9 

ns 

1 

192 . Ti 

170.0 

162.7 

177. .5 

182.1 

194.8 


174.0 

101.3 

174.0 

■■IrTiXillli 

11.02 

2 

134.1 

128.0 

130.2 

128.6 


129.0 

140.3 

134.0 

130.0 

125.8 

131.7 

4.45 

3 

139.0 

123.3 

151.9 

135.1 

152.5 

141.8 

130.0 

104.7 

127.4 

143.6 

141.5 

12.39 

4 

138.0 

138.5 

124.7 

157.4 

130.7 

150.3 

190.0 

141.4 

130.4 

142.5 

146.3 

MlfUfMI 

5 

123.5 

140.8 

142.7 

154.1 

145.2 

154.3 

140.8 

175.3 

150.4 

152.7 

148.0 

13.30 

6 

206.4 

198.2 

208.9 

218.2 

185.1 


230.2 

191.6 

214.2 

212.4 

207.0 

13.13 

7 

158.3 

167.5 

166.6 

210.5 

192.1 

161.7 

178.2 

189.1 

222,3 

169.6 

180.6 


8 

139.0 

133.7 

130.7 

132.4 

166.9 

170.2 

145.1 

155.7 

175.9 

164.5 

150.4 

16.70 

9 

279.8 

240.1 

251.4 

293.8 


276.6 

248.0 

273.5 

281.4 

270.2 

267.5 

16,99 

10 

126.7 

125.4 

135.7 

129.4 

137.9 

136.9 

140.6 

134.6 

134.1 

131.2 

133.3 

4.95 

11 

275.4 

278.7 

290.4 


322.8 

337.2 

295.1 

336.0 

288.4 

316.5 


22.94 

12 

277.0 

292.8 

297.5 

302.3 

286.6 

282.2 

285.0 

268.6 

275.2 

268.7 i 

283.6 

11.56 

13 

183.1 

178.1 

145.9 

152.0 

148.9 

152.8 

151.3 

165.4 

158.6 

148.1 

158.4 

13.02 

14 

155.6 

156.7 

210.0 

219.4 

229.4 

212.4 

230.0 

204.3 

224.0 

178.3 

202.0 

28,42 

15 

153.6 

162.7 

142.7 

121.6 

145.3 

149.2 

143.8 

183.5 

160.2 

157.3 

161.0 

15.84 

16 

195.7 

219.7 

215.4 

216.2 

217.3 

185.4 

218.1 

219.4 

197.0 

210.9 

209.5 

12.22 

17 

178.5 

174.0 

180.2 

194.6 

197.7 

163.1 

182.8 

198.1 

176.7 

214.1 

186.0 

14,91 

18 

134.0 

133.7 

131.3 

133.0 

157.3 

153.3 

146.5 

154.0 

160.6 

187.1 

148.1 

16.99 

19 

270.8 

254.3 

268.2 

244.1 

286.2 

268.2 

280.4 

271.0 

249.7 

270.6 

266.4 

13.22 

20 

121.4 

123.2 

125.5 

127.4 

126.1 

124.6 


123.6 i 

125.6 

134.7 

126,2 

3.77 


Table 4.55: Capability study for checking the process capability at one of minimum point in 
initial capability study 
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Figure 4.31: Control chart for average of on- 
line performance 


Figure 4.32: Control chart for 
viation of on-line performance 


standard de- 


4.4 Test Function 4: Welded Beam Problem 

4.4.1 Initial Capability Study 


To check the need for variation reduction capability study is done. The on-line performance 
ureasured for lent 50 % of the generation are shown in Table 4.56. The average and standard 


deviation of each group are also shown in Table 4.56. 


IVIal 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 
1 1 
12 

13 

14 

15 

16 

17 

18 

19 

20 


I 

380.1 

26.2 

9.1 

21.2 

23.4 

426.7 

32.4 

18.6 

2179.2 

1176.2 

600.0 

1354.7 

1745.1 

35.4 

259.7 

1044.8 

15.9 

9.1 

512.8 

392.1 


2 “ 

376.1 

24.5 

9.6 

17.2 

20.6 

408.1 

31.9 

26.3 

2366.7 

782.7 

470.0 

1876.0 

2143.4 

34.8 

264.6 

1035.7 

17.0 

12.5 

487.8 

401.2 


349.3 

22.1 

9.8 

16.6 

20.9 

440.5 

31.8 

19.4 

2194.6 

948,8 

510.2 

2043.6 

1806.0 

35.6 

266.2 

1008.3 

17.4 

9.3 

440.6 

460.0 


333.0 

23.7 

5.3 

18.7 

10.0 

437.4 

31.6 

20.5 

1972.8 

1202.1 

481.9 

1179.0 

2187.8 

34.2 

253.7 

1112.0 

15.0 

8.4 

494.2 

372.2 


Run n umber 


374.0 

22.9 

5.1 

16.9 

14.1 

444.5 

30.7 

21.1 

1416.1 

1176.3 

556.2 

1839.0 

1718.1 

36.5 

272.1 

981.1 

15.6 

10.8 

458.1 

419.3 


358.0 

24.3 

10.8 

17.9 
22-3 

397.8 

29.9 

22.4 

2382.5 

1037.9 

433.2 

1900.5 

2276.7 

36.5 

275.0 
X050.0 

16.9 

II. 1 

477.7 

423.4 


364.6 

23.2 

7.0 

23.4 

17.6 
496-3 

30.6 

17.1 

2359.7 

857.7 

498.8 

1908.1 

2047.2 

35.7 

269.9 

976.8 

16.9 

7.4 

500.8 
366.4 


300.2 
22.8 

3.3 

18.7 
9.2 

441.3 

32.3 
21.2 

1798.3 

957.0 
474.8 

1657.7 

1760.4 

33.7 
254.7 

1042.5 

16.4 
9.9 

487.3 

374.1 


22.6 

3.6 

17.9 

9.0 

432.1 

31.8 

22.6 

1761.8 

987.7 

465.0 

1707.3 

1981.2 

34.6 

241.5 

945.5 

15.0 

11.1 

424.9 

430.1 


30.4 

3.9 

16.8 

14.3 

459.1 

33.0 

22.6 

2036.6 

1039.6 

608.3 
1943-5 

1814.5 

35.9 

280.9 

917.9 
17-2 

7.8 

466.3 

441.1 


average 


363.5 

24.2 
0.7 

18.5 

16.1 

438.4 

31.6 

21.2 

2046.6 

1016.5 

510.5 

1741.0 
1948-0 

35.3 

263.8 

1011.5 

10.4 

9.7 

475.1 

408.0 


Table 4.56: Initial capability study for welded beam problem 


-liKT 

dev. 


31.62 

2.49 

2.87 

2.17 

5.58 

27.02 

0.93 

2.64 
310.23 
139.42 

59.27 

270.31 

206.97 

0.95 

11.09 

57.23 

0.82 

1.64 
27.57 
31.81 


instruct tlie average and standard deviation control charts shown 

These data were used to con 

in Figure 4.33 and Figure 4.34. 
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Figure 4.33; Control chart for average on- Figure 4.34: Control chart for standard dcvi- 
line pcrfonnancc before variation reduction ation of on-line performance before variation 
on welded beam problem reduction on welded beam problem 


VVe see from Figure 
So there is need to 
in previous chapter 


4.33 and FiguK! 4.34 that lot of points are going outside the control limit, 
reduce the variation. For reducing the variation four cases as described 
is selected and for each case variation reduction technique is applied. 


4.4.2 CASE 1 (Proportionate Selection and Single Point Crossover) 

27 trials experiment data are shown in Table 4.57. 



can 

didatc input 
variables 

kiin number 

mB&Bm 


CP 

MP 

1 

2 

3^ 


5 - 


7^ 

g -' 

g- 




1 

100 

0.50 

0.100 

39.5 

37.9 

40.0 

38.3 

37,3 

38.0 

37.2 

37.1 

39.7 

39.6 

■KSEOi 

1.16 

2 

100 

1.00 

0.001 

7.5 

6.3 

7.8 

8.0 

6.9 

7.6 

8.1 

7.3 

7.5 

6.3 

7.34 

0.65 

3 

10 

0.50 

0.100 

28.6 

37.2 

32.0 

34.8 

25.5 

31.6 

26.4 

32.0 

27.1 

42.1 

31.72 

5.22 

4 

100 

0.60 

0.001 

4.3 

7.2 

8.8 

8.4 

4.6 

6.7 

10.0 

10.2 

7,5 

4.7 

7.22 

2.17 

5 

10 

1.00 

0.001 

13.9 

11.3 

13.1 

7.1 

6.4 

7.8 

5.3 

11.0 

11.5 

8.0 

9.54 

2.99 

6 

10 

1.00 

0.100 

34.7 

33.1 

29.4 

33.8 

40.0 

44.2 

28.8 

28.4 

29.0 

29.6 

33.11 

5.35 

7 

100 

1.00 

0.100 

39.8 

40.5 

41.4 

40.9 

40.2 

40.9 

39.9 

41.5 

41.9 

41.6 

40.86 

0.72 

8 

10 

0.50 

0.001 

11.6 

18.7 

12.1 

14.5 

8.7 

4.9 

7.5 

9.0 

7.3 

9.5 

10.38 

3.98 

9 

55 

0.75 

0.051 

22.2 

24.6 

23.1 

23.7 

22.9 

24.2 

23.7 

22.1 

24.0 

25.0 

23.65 

0.98 

10 

55 

0.50 

0.001 

4.3 

4.3 

4.0 

4.1 

4.5 

4.4 

4.1 

4.6 

5.0 

4.4 

4.38 

0.28 

11 

55 

0.60 

0.100 

35.3 

34.8 

38.1 

37.1 

36.6 

36.0 

35.9 

33.7 

36.4 

34.0 

35.85 

1.30 

12 

65 

1.00 

0.001 

4.0 

3.8 

3,1 

2.8 

3.2 

3.1 

2.7 

3.2 

3.1 

3.1 

3.21 

0.39 

13 

55 

1.00 

0.100 

37.7 

38,3 

35.6 

39.6 

38.7 

40.5 

38.3 

37.2 

38.2 

39.4 

38.36 

1.36 

14 

10 

0.75 

0.001 

14.6 

9.1 

11.8 

12.1 

8.1 

6.9 

8.5 

7.3 

10.2 

8.3 

9.66 

2.44 

15 

10 

0.75 

0.100 

37.9 

35.6 

38.2 

32.8 

31.3 

29.0 

32.8 

27.6 

29.5 

30.0 

32.47 

3.71 

16 

100 

0.76 i 

0.001 

7.7 

6.3 

8.0 

8.9 

, 6.2 

6,9 

8.2 

8.7 

9.8 

8.8 

7.96 

1.17 

17 

100 

0.75 

0.100 

38.2 

40,6 

41.2 

39.1 

40.1 

40.4 

38.8 

39.0 

39.4 

40.0 

39.68 

0-93 

18 

10 

0.50 

0.051 

15.3 

28.2 

23.3 

16.9 

24.3 

17.1 

33.8 

17.8 

18.8 

35.1 

23.06 

7.22 

19 

10 

1.00 

0.051 

23.9 

33.6 

17.2 

29.6 

22.2 

18.5 

33.2 

26.0 

25.3 

28.1 

25.77 

5.57 

20 

100 

0.50 

0.051 

25.2 

24.0 

24.0 

24.1 

23.8 

25.0 

24.4 

24.5 

23.6 

25.0 

24.36 

0.64 

21 i 

100 

1.00 

0.051 

25.1 

27.0 

25.8 

29.1 

26.6 

29.9 

26.3 

26.6 

27.9 

26.6 

27.09 

1.46 

22 

55 

0.75 

0.001 

3.9 

3.7 

3.4 

3.7 

4.0 

4.1 

3.3 

3.5 

3.6 

3.6 

3.07 

0.25 

23 

55 1 

0.75 

0.100 

35.2 

36.7 

37.8 

36.8 

38.2 

37.8 

36.6 

36.3 

37.9 

37.0 

37.02 

0.92 

24 

10 

0.75 

0.051 

15,4 

15.9 

19.0 

16.5 

17.2 

19.9 

15.7 

14.7 

18.0 

25.9 

17.81 

3.28 

25 

100 ! 

0.75 

0,051 

24.0 

24.8 

24.2 

24.5 

27,0 

25.9 

25.1 

25.2 

26.1 

25.5 

25.24 

0,93 

26 

65 i 

0.50 

0.051 

21.8 

23.0 

24.9 

22.7 1 

23.0 

21.6 

23.8 

20.7 

22.3 

23.1 

22.69 

1.18 

27 

56 

1.00 

0.051 

24.4 

24.6 

25.3 

26.0 i 

24.9 

25.6 

24.0 

23.1 j 

25.6 

25.8 

24.94 

0.93 


Table 4.57; Data obtained from 27 trials on welded beam problem case 1 
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The significant effects of the [)araincters, their interactions and the quadratic effects are 
shown in Table 4.58 for average and in Tabic 4.59 for standard deviation. 


Design parameter 

Effect 

MP 

mutation prob. 

29.352 

PS*MP 

population size * mutation prob. 

4.8 

PS 

population size 

2.74 

MP*MP 

mutation prob. (quad) 

2.11 

CP*MP 

crossover prob.* mutation prob. : 

1.36 

PS*PS 

population size(quad.) 

-1.35 

CP 

crossover prob. 

1.34 


Table 4.58: Effects table for average case 1 


Design parameter 

Effect 

PS 

population size 

-3.33 

MP*MP 

mutation prob. (quad.) 

-1.9 

PS*MP 

population size * mutation prob. 

-1.00 


Table 4.59: Effects table for standard deviation case 1 
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Response surface study 


rrh the average is fitted. Tfie 

Based on 27 trials data in Table 4.57, the best surface throUS* 
resulting equation is : 


on-line performance 


10.522 - 0.09798PS' - 0 
-18.4Mp2-f-69.552CP 


SUCF" + 0.00067PS2 


The best fit surface through standard deviation is: 
log((7) = 1.0272 - 0.01377P5 + 53.68Mp2 - 0.0341P5 



Analysis by direct observation 


• Set PS = 100 and MP = 0.0505 in order to minimize 


• Set CP = 0.5 in order to reduce the on-line perforinarxc®- 


Analysis by Variation transmission 

^ equation 4.34 gives: 

Minimization of equation of standard deviation obtained froff^ 

a = 1.168 at PS = 10, CP = 0.818, MF = 

j. jon 4.35 will be 2.966, at ps 

the standard deviation of on-line performance from eq^J- 

.-„ri found from equation 435 

= 10, CP = 0.818 and MP = 0.064 . The standard 

It analysis of variation 

by direct observation method is 0.68. This means that re^ 
transmission method is not good. 

r\ 0505 and standard deviation 

So with mean at PS = 100, CP = 0.5 and MP = 

, , is expected 

^ps = 5, acp = 0.05, and cmp = 0.005 final capability study 


improvement. 
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4. 4. 2.1 Final Capability Study for Case 1 

deviation is 
the control 

limits hence the process is capable. Now we try for next case. 


I'liti data are given in 'I'able 4.Gt). The control chart for average and standard 
also given in Figure 4.35 and Figure 4.36 . We see that all points are within 


Trial 

number 

Run number | 

average 

std. 

dev. 

I 

5 

3 

4 

5 

g 

7 

g 

S 

10 

1 

23.9 

26.3 

25.6 

25.4 

27.2 

28.1 

26.1 

26.7 

25.0 

25.2 

26.0 

1.19 

2 

23.6 

22.6 

22.1 

23.2 

21.2 

23.5 

22.9 

23.6 

23.4 

23,2 

22.9 

0.75 

3 

21.1 

21.4 

20.5 

21,0 

22.4 

23.1 

21.7 

21.6 

23.2 

21.6 

21.8 

0.89 

4 

23.4 

25.4 

25.7 

24.0 

22.0 

24.2 

24.9 

23.1 

23.2 

23.2 

23.9 

1.17 

5 

21.5 

21.0 

21.1 

20.0 

23.4 

20.3 

21.9 

21.6 

20.9 

23.6 

21.5 

1.18 

6 

24.1 

24.6 

23.6 

23.8 

24.4 

26.8 

24.8 

24.3 

26.7 

24.9 

24.8 

I.IO 

7 

25.0 

22.8 

23.7 

22.9 

24.3 

25.4 

24.3 

25.1 

25.3 

25.3 

24.4 

0.99 

8 

22.4 

21.5 

19.7 

23.7 

21.8 

24.5 

22.1 

20.8 

22.8 

22.4 

22.2 

1.37 


27.0 

28.8 

27.5 

20.7 

27.6 

27.6 

20.2 

27,3 

27.3 

25.1 

27.2 

1.00 

10 

23.5 

23.9 

26.3 

25.3 

24.2 

24.5 

26.9 

23.2 

23.9 

24.6 

24.6 

1.20 

11 

26.9 

26.3 

27.2 

30.0 

27.6 

27.8 

27.9 

30.1 

27.9 

28.4 

28.0 

1.22 

12 

24.4 

25.7 

27.7 

27.0 

28.3 

29.9 

28.3 

26.5 

28.9 

26.5 

27.3 

1.63 

13 

24.0 

26.3 

25.5 

24.9 

25.3 

23.8 

27.6 

24.9 

24.0 

27.5 

25.4 

1.38 

14 

22.4 

24.4 

22.7 

21.8 

24.2 

23.0 

23.1 

23.2 

24.2 

25.5 

23.4 

1.10 

15 

21.6 

22.3 

20.7 

21.2 

20.5 

20.7 

21.1 

20.7 

21.5 

20.6 

21,1 

0.67 

16 

25.5 

24.4 

25.2 

25,7 

24.8 

25.5 

23.7 

27.0 

25.0 

23.8 

25.1 

0.09 

17 

24.6 

22.7 

24.8 

25.2 

25.7 

23.1 

25.7 

25.4 

26.0 

24.7 

24.8 

1.09 

18 

22.5 

22.2 

23.2 

23.4 

21.2 

20.6 

23.6 

21.2 

21.0 

21.4 

22.0 

1.08 

19 

27.0 

25.9 

25.1 

27.4 

27.1 

27.8 

27.9 

28.1 

28.2 

27.9 

27.2 

1.02 

20 

23.6 

24.7 

24.6 

25.9 

23.0 

24.2 

23.9 

23.7 

23.0 

23.1 

24.0 

0.92 


Table 4.60; Final capability study for case 1 welded beam problem 




Figure 4.35: Control chart for average of on- Figure 4.36; Control chart for standard devi- 
line performance after variation reduction ation of on-line performance after variation 
for case 1 on welded beam problem reduction for case 1 on welded beam problem 


4.4.3 CASE 2 (Proportionate Selection and Uniform Crossover) 

27 trial experiment data are shown in Table 4.61. 
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no. 

Input 

variables 

ftun number 



PS 



j— 

2” 

3^ 

If- 


ZZQ 



zzo 

miyiij 

ECEansii 


1 

100 

0.50 

"" O.ioo 

42.1 


4t.4 


41 ft 


40 3 

■ M'a " 

i 

i 

4t.il 


■HHIHIJBI 

U 

too 

1 .00 

O.tMII. 

.1 n 

0 rt 

« a 

ft 0 

8 0 


0 8 

4.0 

fi.3 

8.7 


!,|i2 

a 

10 

0.50 

0 ItM) 

.1«.5 

27. 1 

30,3 

33.8 

27.5 

WBm 

26.8 

28.3 

20.3 

29.7 

31.10 

5.77 

4 

100 

0.50 

O.OOl 

5.7 

5.6 

5.3 

3.7 

4.9 

■xi 

6.4 

4.7 

6.3 

0.3 

5.70 

1.16 

5 

10 

1 .00 

0.001 

11. 1 

18.5 

22.7 

7.2 

5.6 

7.4 

8.2 

7.4 

6.9 

16.6 

11.16 

5.94 

6 

10 

1.00 

0.100 

29.3 

32.1 

28.3 

27.3 

27.2 

51.7 

32.6 

28.5 

28.9 

31.1 

31.70 

7.30 

7 

100 

1.00 

0.100 

43.6 

45.0 

45.6 

43.6 

41.8 

42.6 

45.9 

44.8 

47.7 

43.6 

44.43 

i.73 

8 

10 

0.50 

0.001 

7.5 

8.5 

7.9 

8.1 

8.9 

5.9 

8.1 

9.4 

13.3 

5.3 

8.29 

2.17 

9 

55 

0.75 

0.051 

25.0 

25.8 

22.6 

25.2 

26.3 

26.0 

25.7 

23.3 

22.6 

24.1 


1.40 

20 

55 

0.50 

0.001 

3.8 

3.2 

3.0 

3.1 

3.1 

3.1 

3.4 

3.1 

3.0 

3.0 

3.16 

0.24 

11 

55 

0.50 

0.100 

39.9 

38.2 

39.3 

37.5 

36.9 

37.8 

39.6 

36.1 

38.4 

38.1 

38.19 

1.20 

12 

55 

1.00 

0.001 

5.2 

3.2 

3.0 

3.1 

2.9 

3.3 

3.2 

3.0 

3.3 

3.1 

3.35 

0.86 

13 

55 

1.00 

0.100 

40.4 

39.9 

40.9 

39.7 

40.9 

42.7 

40.6 

39.8 

40.2 

42.0 

40.70 

0.97 

14 

10 

0.75 

0.001 

11.2 

15.8 

12.8 

9.2 

4.6 

9.7 

8.6 

7.8 

7.7 

16.6 

10.40 

3.76 

15 

10 

0.75 

0.100 

86.7 

29.9 

26.9 

26.6 

33.8 

29.4 

34.1 


24.6 

35.1 

35.74 

18.25 

16 

100 

0.75 

0.001 

5.4 

7.2 

8.0 

7.8 

7.0 

7.8 

7.2 

11.4 

10.1 

4.3 

7.61 

2.05 

17 

100 

0.75 

0.100 

43.5 

41.2 

43.9 

44.4 

40.7 

42.1 

44.2 

42.7 

42.3 

42.2 

42.72 

1.25 

18 

10 

0.50 

0.051 

17.1 

23.7 

32.8 

23.6 

21.1 

30.9 

16.4 


26.6 

18.8 

22.72 

5.93 

19 

10 

1.00 

0.051 

24.3 

23.6 i 

22.7 

26.7 

17.5 

24.7 

25.6 


17.2 

17.1 ^ 

21.94 

3.69 

20 

100 

0.50 

0.051 

25.8 

28.7 

25.1 

24.8 

25.6 

26.8 

28.5 

24.1 

26.1 

26.4 

26.19 

1.49 

21 

100 

1.00 

0.051 

27.7 

27.1 

29.9 

27.6 

27.7 

29.8 

30.1 

27-9 

29.4 

28.9 

28.60 

1-12 

22 

55 

0.75 

0.001 

3.0 

3.0 

3.1 

3.3 

5.5 

3.1 

3.1 

10.7 

2.9 

3.2 

4.08 

2.46 

23 

55 

0.75 i 

0.100 1 

38.8 

37.4 

39.9 

40.5 

40.5 

40.1 

37.5 

38.7 

39.8 

39.1 

39.24 

1.13 

24 

10 

0.75 

0.051 

31.4 

21.7 

21.2 

16.1 

17.6 

30.5 

18.8 

17.3 

15.5 

19.6 

20.98 

6.62 

25 

100 

0.75 

0.051 

26.5 

27.3 

28.6 

26.3 

26.6 

26.7 

28.2 

28.3 

26.2 

26.0 

27.06 

0.96 

26 i 

55 

0.50 

0.051 

21.9 

23.3 

24.7 

23.2 

24.1 

25.6 

21.9 

21.7 

24.7 

26.3 

23.74 

1.62 

27 i 

55 

1.00 

0.051 

24.6 

26.0 

25.2 

26.8 

25.5 

27.5 

24.1 

23.9 

26.7 

26.1 

25.64 

1,20 


Tabic 4.61: Data obtained from 27 trials on welded beam problem case 2 


The significant eft'ects of the parameters, their interactions and the quadratic effects are 
shown in Table 4.62 for average and in Table 4.63 for standard deviation. 


Design parameter 

Effect 

MP 

mutation prob. 

31.65 

PS*MP 

population size * mutation prob. 

6.65 

PS 

population vsize 

3.95 

MP*MP 

mutation prob. (quad) 

2.12 

CP 

crossover prob. 

1.52 

PS*PS 

population size (quad) 

-1.004 


Table 4.62: Effects table for average case 2 
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Design parameter 

Effect 

PS 

population size 

-5.06 

PS*MP 

{population size * iniitation prob. 

-3.4 

PS*PS 

poptilaiion sizc(c|ua(I.) 

-2.75 

MP 

mutation prob. 

2.04 

CP*CP 

crossover prob. (quad.) 

1.581 


Table 4.63: Effects table for standard deviation case 2 


Response surface study 

Based on 27 trials data in Table 4.61, the best surface through the average is fitted. The 
resulting equation is : 

on-line performance = 7.179 - 0.083P5 -I- 3.167CP -1- 232.29MP 4- 0.000482P52 

+72.2MP^ -b 1.469P5 x MP (4.36) 


The best fit surface through standard deviation is: 

log(a) = 1.6589 - 0.062P5 + 10.81AfP-b0.000485P52 4- 0.2289C'P2 

-0.1269P5 X MP (4.37) 

Analysis by direct observation 

• Set PS = 62 , CP = 0.5 and MP = 0.001 in order to minimize variation. 

Analysis by Variation transmission 

Minimization of equation of standard deviation obtained from equation 4.36 gives : 
cj = 1.275 at PS = 10, MP = 0.033 

Now putting these values in the equation 4.36 , and minimizing that will give value of 


CP at which we have to target it. 
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minimum on-line performance = 16.2 at CP = 0.5. 

standard dcviat.ion from e(|uation 4.37 will he 4.3 at PS=10 , CP =0.5 and MP = 
0.033. 

From direct observation method standard devaition of on-line performance = 0.77. We 
can see that result from both the method doesn’t match. 

So with mean at PS = 62, CP = 0.5 and MP = 0.001 and standard deviation aps = 
5, acp = 0.05, and cmp — 0.005 final capability study is done to verify the expected 
improvement. 

4. 4. 3.1 Final Capability Study for Case 2 

The data are given in Table 4.64. The control chart for average and standard deviation is 
also given in Figure 4.37 and Figure 4.38 . We see that all points are within the control 
limits hence the process is capable. Now we try for next case. 

Table 4.64; Final capability study for case 2 welded beam problem 


Trial 

Run number j 

average 

mmm 

number 

1 



4 


6 

7 


5 


1 

4.4 

WgQII 


7.5 

6.2 

4.2 

6.3 

WSEM 

10.9 

11.6 

8.2 

MESM 

2 

3.6 



3.4 

3.4 

3.5 

3.3 

ms 

4.2 

3.2 

3.6 

0.70 

3 

3.8 

3.2 

8.8 

3.2 

3.1 

3.0 

4,1 

3.0 

2.9 

2.7 

3.8 

1.81 

4 

4.1 

3.6 

3.6 

4.1 

4.3 

3.9 

2.9 

3.0 

2-7 

2.7 

3.5 

0.62 

5 

4.2 

3.0 

2.7 

2.8 

2.7 

2.9 

2.7 

2.7 

2.8 

2.9 

2.9 

0.45 

6 

4.3 

5.4 

5.4 

4.6 

5.0 

5.1 

5.5 

5.2 

4.7 

5.5 


0-39 

7 

4.4 

4.5 

4.5 

5.0 

4,0 

4.2 

5.0 

3.9 

4.6 

5.2 


0.43 

8 

4.2 

2.7 

6.8 

6.2 

9.7 

2.7 

mm 

5.3 

8.8 

3.2 


2.53 

9 

6.9 

9.2 

7.6 

7.2 

7.5 

7.4 

K la 

7.7 

7.7 

7.7 

7.6 

■iHiif 1 

10 

3.6 

4.1 

4.6 

0.3 

7.7 

7.1 

■H 

9.8 

7.4 

7.3 

6.2 

2.05 

1 1 

7.4 

12.2 

11.8 

8.0 

9.9 

8.5 

BEI 

11.1 

8.7 

8.8 

10.0 

2.15 

12 

6.8 

9.0 

11.0 

7.2 

6.8 

16.9 

10.8 

8.5 

9.4 

10.1 

9.6 

2.99 

13 

4.1 

5.0 

7.3 

5.1 

5.4 

3.8 

5.0 

9.5 

5,0 

8.0 

5.8 

1.83 

14 

7.9 

12.7 

4.1 

6.7 

3.5 

3.4 

3.2 

5.4 

6.4 

3.4 

5.7 

2.96 

16 

4.2 

3.0 

2,8 

2.7 

2.8 

2.8 

2.8 

2.9 

2.8 

2.9 

3,0 


16 

6.0 

6.0 

6.9 

5.0 

5.7 

4.7 

5.8 

6.4 

4.7 

5.6 



17 

3.9 

4.3 

6.4 

4.3 

3.7 

4.2 

3.9 

3.7 

3.6 

4.1 


0.79 

18 

3.6 

3.1 

9.8 

7.7 

9.9 

2.6 

8.5 

9.9 

5.3 

5.0 


2.95 

19 

9.4 

8.0 

8.2 

6.5 

8.9 

9.1 

6.7 

10.2 

8.9 

8.0 

8.4 

1.15 

20 

10.7 

4.1 

4.5 

6.7 

5.0 

6.6 

3.7 

6.2 ; 

6.5 

3.6 

6.8 

2.13 


4.4.4 CASE 3 (Tournament Selection and Single Point Crossover) 


27 trials experiment data are shown in Table 4.65.. 























Average 




welded beam problem case 3 


Table 4.65: Data obtained from 27 trials on 
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The signilicaut cirecis of the parameters, their interactions and tlie quadratic effects are 


shown 


in Table 4. (id for avc'rage 


and in Tables 4.G7 for standard devi atioii. 


Design parameter 

Effect 

MP 

mutation prob. 

-1445.31 

PS 

I)opulati()n size 

1411.67 

PS+MP 

l)opulation size * mutation prob. 

-1231.45 

MP*MP 

mutation prob. (quad) 

-234.26 

PS*PS 

population size (quad) 

-125.246 

PS*CP 

population size * crossover prob. 

36.756 

CP*CP 

crossover prob. (quad.) 

32.63 


Table 4.GG: Effects table for average case 3 



l)(hsij.!;u pa.raiu(^t(n' 

Effect 

MP' 

mutation prob. 

-210.4 

PS*MP 

I)opulation size * mutation prob. 

73.01 

MP*MP 

mutation prob. ((juad) 

-46.488 

PS*CP 

population size * crossover prob. 

34.4 

CP 

crossover prob. 

-32.178 

CP*MP 

cro.s.sover prob. * mutation prob. 

27.4 

CP*CP 

crossover prob. (quad.) 

19.65 

PS*PS 

population size (quad.) 

-16.05 

PS 

population size 

-9.49 


Table 4.G7: Effects table for standard deviation case 3 


Response surface study 

Based on 27 trials data in Table 4.G5, the best surface through the average is fitted. The 
resulting equation is : 

on-line performance = 231.01C + 21.3F5' — 4415.04MP -f 0.0639P5 — 88.35CP 

-f 54G69.87MP2 + 2.129P5 x CP - 282.48PS' x MP (4.38) 

The best fit surface through standard deviation is: 


log(cr) = 8.588 + 0.034P5-8.79GCP-52.G93MP-0.0003PS^-b4.668CP^ 




4.4 'rev'll. Fimc.t.ion 4: W(4ckul I]ea,in Problom 

-dOSoOMP'^ ~ 0.0031P5 X CP + 27.428CP x MP 
+0.2531 P.9 X MP 

Analysis by direct observation 

• Set PS = 10 , CP = 0.75 and MP = 0.1 in order to minimize variation. 

Analysis by Variation transmission 

Minimization of ecination of standard deviation obtained from equation 4.38 gives ; 

a = 8.075 at PS = 19 ,CP = 0.5 and MP = 0.088 

Putting the value PS = 19, CP = 0.5 and MP = 0.088 in the equation 4.39 of standard 
deviation of on-line performance we get a = 8.77. The standard deviation found from 
equation 4.39 l)y direct observation method is 2.85. This means that the result from variation 
transmission method is not good. 

So with mean at PS = 10, CP = 0.75, and MP = 0.1 and standard deviation ctp.s = 
5, (7cp = 0.05, and cmp = 0.005 final capability study is done to verify the expected 
improvement. 

4.4.4. 1 Final Capability Study for Case 3 

The data arc given in Table 4.68. The control chart of average and standard deviation is 
also given in Figure 4.39 and Figure 4.40. We see that all points are within the control limits 
ex(q>t two points whi(4i ar(^ slightly alrove UCL in Figure 4.40 hence the ])rocess is capable. 
Now we t.ry for next case. 

4.4.5 CASE 4 (Tournament Selection and Uniform Crossover) 

27 trials experiment data are shown in Table 4.69. 
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'IVinl 

!iumbcr 


Run number | 

average 

Btd. 

1 

2 

3 

4 

5 

6 

7 

8 


10 


dev. 

1 

149.55 

146.51 

147.84 

151.61 

154.60 

153.50 

157.62 

156.84 

160.63 

162.38 

153.11 

~TTr“ 

2 

112. 14 

106.74 

107.90 

113.61 

112.98 

105.61 

106.03 

104.38 

113.29 

117.14 

109.98 

4.35 

3 

107.47 

108.45 

126.55 

113.90 

112.97 

113.36 

102.17 

105.50 

109.13 

111.86 

111.14 

6.59 

-1 

88.45 

74.97 

68.95 

63.97 

77.03 

87.31 

95.86 

83.53 

81.19 

63.37 

78.46 

10.83 

5 

111.10 

106.94 

119. 66 

106.39 

111, 07 

107.00 

111.60 

104,66 

112.23 

113.18 

110.37 

4.37 

G 

122.48 

122.23 

107.06 

144.25 

112.55 

120.72 

131.84 

136.77 

100.68 

130.05 

122.86 

13.49 

7 

85.25 

112.40 

82.34 

87.59 

85.47 

111.71 

82.82 

99.71 

78.59 

104.52 

93.04 

12.80 

8 

116.09 

104.88 

121.34 

110.65 

117.10 

109.29 

108.68 

103.61 

100.02 

120.65 

111.23 

7.33 

9 

161.98 

139.02 

146.82 

120.02 

142.10 

148.17 

142,71 

129.92 

145.17 

158,04 

144.06 

10.94 

to 

106.80 

104.06 

105.63 

1 10.58 

102.91 

107.. 39 

113.15 

103.80 

104,45 

109.49 

106.82 

.1.35 

t 1 

192.94 

197.07 

188,24 

203.24 

201.66 

179.99 

216.33 

199.93 

208.80 

201.12 

198.83 

10.06 

12 

184.34 

198.28 

187.94 

195.06 

193.09 

188.93 

170.78 

193.42 

188.95 

187.59 

188.84 

7.67 

13 

128.77 

134.42 

134.55 

135.04 

137.35 

133.09 

133.01 

130.11 

135.62 

138.53 

134.05 

2.98 

14 

118.79 

111.76 

102.89 

113.29 

114.55 

112.54 

102.56 

113.67 

103.24 

110.89 

110.42 

6.60 

16 

115.23 

111.44 

105.62 

108.44 

102.89 

121.50 

116.27 

115.49 

102.00 

116.00 

111.49 

6.51 

IG 

127.69 

119.19 

109.91 

139.35 

102.01 

122.60 

139.80 

144.40 

101.99 

123.86 

123.09 

15.28 

17 

84.58 

109.01 

80.50 

85.19 

89.63 

109.02 

81.65 

99.79 

87.26 

104.07 

93.13 

1 1 .33 

18 1 

1 17.86 

U1.05 

105.05 

115.66 

113.94 

112.92 

113.66 

110.79 

107,70 1 

I1T.24 

112.69 

4.06 

19 

144.99 

148.32 

140.74 

131.06 

180.59 

142.64 

125.89 

133.99 

136.59 : 

161.58 

144.64 

16.08 

20 1 

107.64 

115.67 

104.76 

111.84 

114.71 

106,45 

106.55 

102.00 

101.29 , 

106.79 

107.77 

4.90 


d'able 4. 68: bdnal ca.pability study for case 3 welded beam problem 



Trial Number Trial Number 


Figure 4.39: Control chart for average of on- Figure 4.40: Control chart for standard devi- 
line performance after variation reduction ation of on-line performance after variation 
case 3 on welded beam problem reduction case 3 on welded beam problem 
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178.65 

13 

5 5 

1 ,00 

0.1 oo 

313.4 

304.7 

325.3 

303.3 

278.5 

373.8 

287.1 

311.0 

378.9 

360.8 

323.68 

35,54 

14 

10 

0.75 

0.001 

452.4 

440,3 

274.4 

402.7 

441,9 

492.4 

327.7 

321.9 

329.8 

348.3 

372.78 

75.54 

ir. 

10 

0.75 

0.100 

12741 

137.1 

138.4 

140.0 

130.1 

144.4 

134.2 

133.1 

126.3 

136.0 

134.82 

5.76 

K) 

1 00 

0.75 

0.001 

2048,2 

4447.3 

2000.8 

3506.0 

3459.5 

3161.7 

3173.9 

3443.2 

3359.8 

3383.4 

3274.98 

201.41 

17 

100 

0.75 

0.100 

400. 1 

440.7 

405.4 

404.7 

410.7 

378.5 

456.9 

428.0 

412.5 

474.7 

420.82 

28.67 

18 

10 

0.50 

0.051 

202.2 

207,0 

102.3 

200.3 

191.8 

190.1 

213.4 

198.3 

193.7 

177.0 

196.62 

10.09 

10 

10 

1 ,00 

0,051 

184.1 

182.1 

102.8 

199.6 

198.2 

190.3 

180.3 

191.6 

187.2 

194.0 

189.93 

6,64 

20 

100 

0.50 

0.051 

14.50.2 

1447.4 

1471.0 

1260.1 

1347.4 

1404.7 

1416.6 

1332.1 

1414.0 

1358.3 

1372.05 

54.75 

21 

100 

1.00 

0,051 

1242.7 

1448.0 

1420.5 

1244.7 

1285.6 

1235.6 

1237.5 

1505.7 

1301.3 

1307.8 

1322.03 

97.71 

22 

55 

0.75 

0,001 

1501.5 

1400,7 

1580.1 

1565.0 

1342.4 

1386.7 

1595.7 

1597.8 

1309.6 

1289.4 

1474.95 1 

129.07 

23 

55 

0,75 

0,100 

405.5 

440.8 

417.5 

350.4 

294.0 

294.3 

305.5 

305.8 

312.7 

383.0 

320.85 ’ 

28.48 

24 

10 

0.75 

0 4)5 1 

200, 7 

177 0 

204.7 

105.4 

189.3 

207.8 

194.8 

184.5 

192.1 

179.2 

192.55 

10.31 

20 

100 

0.75 

0.051 

1 102.5 

141,5 5 

1414,2 

1425. 6 

1442.3 

1525,6 

1263.2 

‘ 1376.2 

14,35.4 

1251.9 

1351.25 

109 46 

20 

n n 

0.50 

0.051 

1044.5 

041,0 

01 7.0 

670.1 

625.0 

656.9 

613.8 

715.5 

830.5 

767.9 

780.29 

1 42.07 

27 


1 .00 

04)51 

7204) 

074,0 

050.4 

002.4 

847.6 

661.1 

625.8 

647.1 

838.0 

848.9 

711.63 

96.89 


'I'a.hh' ddi!): Did. a oht.aiiu'd (Voin 27 trials on wedded beam problem ease 4 


Thee signilieani, (dleeXs of tJie parauudx'rs, their interaction and the. quadratic effects are shown 
in Table' 4.70 for a,ve'rag'ee anel in 'I’able 4.71 for standard deviatie)n. 


Design parameter 

Effect 

FS 

population size 

1449.16 

MB 

inutc'iiiou prol). 

-1414.66 

P.S*MP 

population size * mutation prob. 

-1311.74 

MP*MP 

mutation prol). (eiuad) 

-235.22 

PS*PS 

pe)[)ulation size ((piad) 

-119.3 

CP 

e:re)ssovor prob. 

-35.37 

CI>*MP 

e;ros.se)ver prob. * mutation prob. 

29.176 

P.SRIP 

popuhition size * crossover prob. 

-15.7 


Tal)le 4.70: Effects table for average case 4 


Response snrface study 


Based e)u 27 trials dal.a in 'Fable 4. GO, 


f,he best surface through the average is fitte'xl. 


The 


re'seilting eeiuation is : 


= 34.343 + 24.98P.9 + 43.33(77" - 2302.19M P + 0.0GP5^ 


on-line' ix'rfonnance 
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Design parameter 

Effect 

MP 

mutation prob. 

-122.4 

PS 

population size 

77.1G4 

PS*MP 

population size * mutation prob. 

-58.82 

PS*PS 

population size (quad) 

31.345 

MP*MP 

mutation prob. (quad) 

-8.01 

PS*CP ^ 

population size * crossover prob. 

6.994 


V 


Table 4.71; Effects table for standard deviation case 4 

+60357.13MP^ - 0A8PS xCP- 1748.4CF x MP 
-301.38P5 X MP (4.40) 


The l)C.st fit surface; througli standard deviation is: 


log(o-) = 3.319 + 0.0576P5 - 19.287MP - 0.000414PP^ - 46.11MP^ 

+0.0039P5 X CP + 0.045PS x MP (4-41) 


Analysis by direct observation 


• Set PS = 10 , CP = 0.5 and MP = 0.1 in order to ininiinizc variation. 


Analysis by Variation transmission 

Minimization of equation of standard deviation obtained from equation 4.40 gives : 


o = 10.65 at PS = 20, CP = 1.0, and MP = 0.088. 

Putting the value PS = 20, CP =1 and MP = 0.088 in the equation 4.41 of standard 
deviation of on-line performance we get o = 11.11. The standard deviation found from 
equation 4.41 by direct observation method is 4.015. This means that re-suit from variation 

transmission method is not good. 

So with mean at PS = 10, CP = 0.5, and MP = 0.1 and standard deviation aps = 
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5,0'cp — 0.05 and (Jmp — 0.005 final capability study is done to verify the expected improve- 
ment. 

4.4. 5.1 Final Capability Study for Case 4 

The data are given in Table 4.72. The control chart for average and standard deviation is 
also given in Figure 4.41 and Figure 4.42. We see that all points are within the control limits 
except some points which are slightly above UCL hence the process is capable. 


4 Vial 

llun number | 

average 

std. 

dov. 

number 

1 

2 

3 


5 

6 

7 

5 

9 

m 

1 

158.19 

145.41 

162.45 

155.46 

157.73 

152.35 

150.40 

146.47 

150.08 

162.75 

153.79 

6.37 

2 

113.76 

105.90 

107.73 

111.10 

116.56 

113.03 

115.22 

110,88 

117.41 

1 19.53 

113.11 

4.30 

3 

111.77 

101.19 

110.50 

114,30 

112.05 

113.48 

110.84 

111.47 

110.65 

103.26 

109.95 

4.28 

4 

85.30 

147.08 

94,85 

146.78 

63.07 

118,74 

130.73 

137.07 

139.48 

130.86 

120.30 

29.26 

6 

115.27 

110.05 

111.16 

107.08 

108.80 

112.49 

119.24 

113.69 

108.39 

115.10 

112.13 

3.76 

6 

132.67 

127.80 

118.77 

123.21 

103.12 

138.40 

145.03 

107.91 

100.56 

141.67 

124.51 

16,19 

7 

101.60 

145.69 

79.36 

113.46 

105.11 

108.31 

108.81 

95.99 

70.92 

144.10 

107.93 

22.93 

8 

113.22 

121.27 

110.84 

102.03 

108.48 

115.17 

100.74 

110.09 

109.33 

120.70 

112.09 

6.79 

0 

156.25 

140.20 

148.91 

138.40 

142,47 

148.63 

129.92 

120.48 

144.51 

101.05 

144.58 

10,80 

10 

109.33 

106.72 

112.62 

111.39 

113.80 

U8.6B 

108.01 

100.82 

103.17 

105.89 

109.99 

4.41 

U 

205.14 

193.30 

209.74 

193.22 

193.99 

207.46 

205.24 

187.45 

200.05 

207.35 

200.89 

7.97 

12 

186.82 

184.05 

170.08 

177.73 

190.13 

188,97 

187.81 

179.28 

177.61 

173.82 

182.53 

5.71 

13 

129.89 

139.65 

140.04 

132.79 

139.62 

137.72 

130.10 

131.66 

137.36 

133.64 

135.85 

3.64 

14 

110.46 

117.79 

107.70 

113.99 

112.76 

116.41 

114.23 

114.90 

114.66 

111.47 

113.44 

2.94 

15 

112.66 

112.12 

102.32 

105.20 

113.17 

121.60 

115.43 

108.55 

113.31 

114.72 

111.91 

6.44 

16 

132.58 

132.98 

120.91 

137.76 

105,73 

138.03 

146.34 

172.91 

191.33 

119.45 

139.80 

26.46 

17 

114.80 

11.4.87 

70,71 

174.01 

118.12 

160.50 

102.04 

107.48 

79.93 

91.29 

119.98 

34.72 

18 

115.64 

116,37 

112.91 

115.69 

111.03 

111.10 

110.85 

105.68 

113.10 

123.83 

114.29 

4.72 

19 

148.79 

148.07 

152.06 

131.07 

144.41 

141.07 

134.71 

132.78 

101.13 

165.36 

140.12 

11.49 

20 

109.58 

103.47 

112.10 

111.50 

105.53 

104.70 

114.35 

113.64 

106.81 

117.86 

109.96 

4.74 


Table 4.72: Final capability study for case 4 welded beam problem 




Figure 4.41: Control chart for average of on- Figure 4.42: Control chart for standard devi- 
line performance after variation reduction ation of on-line performance alter variation 
case 4 on welded beam problem reduction case 4 on welded beam problem 
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4.4.6 Discussion of Result of Welded Beam Problem 

The Table 4.73 represents the average and standard deviation of the capability study after 
optimization and variation reduction. The value of PS , CP and MP are given at which we 
have to fix it with standard deviation of 5 , 0.05 and 0.005 respectively. 

Table 4.73: Result of welded beam problem 



PS 

CP 

MP 

average 

std. deviation 

CASE 1 

Proportionate selection 
single point crossover 

100 

0.5 

0.0505 

24.381 

1.114 

CASE 2 

proportionate selection 
uniform crossover 

G2 

0.5 



0.001 

5.753 

1.956 

CASE 3 

tournament selection 
single point crossover 

10 

0.75 

0.1 

123.29 

9.086 

CASE 4 

tournament selection 
uniform crossover 

10 

0.5 

0.1 

129.16 

14.211 


We see that if our aim is to minimize only on-line performance then we will choose pro- 
portionate selection as selection operator and uniform crossover as crossover operator (case 
2). But if our aim is to minimize only variation in on-line performance then proportionate 
selection with single point crossover (case 1) will be chosen. But in this case on-line perfor- 
mance is very high. So we will choose the case in which variation in on-line performance is 
less and also on-line performance is less. So for welded beam problem proportionate selection 
with uniform crossover will give better result. 

In every case PS, CP and MP which are to be targeted, vary. Since we found that 
proportionate selection with uniform crossover is better, we will set PS at 62, CP at 0.5 and 
MP at 0.001 as mean with standard deviation of 5 for PS, 0.05 for CP and 0.005 for MP. 











Chapter 5 
Conclusion 


In this thesis, a methodology to obtain GA parameters and operators which will produce 
consistent performance of GAs on a function is designed. The methodology is designed 
bixsed on the robust design technique commonly used in product design. Five different GA 
parameters and operators have been chosen for this study. The result on four functions are 
summarized in the following table. The average column in the table is the average of 20 
trials of mean of on-line performance in 10 runs. The standard deviation is within subgroup 
standard deviation of final capability study. 






average 

std; deviation 

Function FI 

Proportionate selection 
uniform crossover 

105 

0.54 

0.0818 

24.83 

0.414 

Function F2 

proportionate selection 
uniform crossover 

55 

0.5 

0.1 

295.88 

15.23 

Rastrigin’s 

function 

proportionate selection 
uniform crossover 

250 

1.00 

0.1 

379.61 

5.89 

Welded beam 
problem 

proportionate selection 
uniform crossover 

62 

0.5 





We see that in all cases proportionate selection with uniform crossover makes the 
process capable (producing on-line performance within lower and upper bounds). So for 
GA to be capable, proportionate selection with uniform crossover is recommended. Since 
none of the previous approaches mentioned in chapter 1 considered the selection operator 
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and crossover oirerator together as tlie parameter affecting performance of GA, we can not 
compare our results with any past studies. However, there are a number of aspects of our 
studies, which agree with that of previous studies. 

In all the test functions, we see that higher population size is recommended, which 
is also suggested by De Jong (1975) that increasing population size reduces the stochastic 
efiect ol random sampling. Goldberg, Deb and Clark (1992) suggested that the population 
size should depend on the problem being solved (mainly on the function’s nonlinearity). Our 
results also agree that the most complicated problem (Rastrigin’s function) require a large 
population. The value of CP is recommended higher for test function 3. For other test 
function a low value of CP is recommended. The value of MP is suggested higher for test 
function 1, 2 and 3. For test function 4, low value of MP is recommended. Since the optimal 
values of PS , CP and MP are different for different problems, hence it is not recommended 
to use a fixed value of these parameters for every problem. However a variation reduction 
technique as suggested here can be used before applying GA to any arbitrary problem. It 
is observed that for unimodal, well behaved function (such as De Jong function FI), a large 
population size with a large mutation probability is recommended, if the minimization of 
the on-line performance is the goal. For a unimodal function but with flat optimal basin, 
a smaller population size gives consistent performance of GAs. For a multimodal function 
(such as F3) , a large population size with a large crossover probability is recommended. For 
any other type of function, a similar study can be performed and best GA parameters and 
operators can be obtained. > 

To show how this method is applicable to an engineering design problem, the technique 
has been applied to a welded beam design problem. When GAs have been run with suggested 
GA i)arameters and operators consistent GA performance has been observed. 

Since statistical tools such as design of experiment are used to find best parameters for 
a consistent GA performance, the technique is expected to give best GA parameters. We 
propose to develop a software implementing robust design technique suggested in this thesis. 
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Such a software can be used by designers and GA practicers to find the best setting of GA 
paraineters and operators. 
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